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— =?  The  analysis  presented  in  this  report  focuses  on  three  attributes  of  the 
inlet  environment:  (a)  the  inlet  hydraulics,  (b)  the  longshore  currents 
adjacent  to  the  inlet,  and  (c)  the  seasonal  morphologic  change  of  the  North 
Inlet  tidal  deltas  and  adjacent  beaches. 

North  Inlet  is  hydraulically  ebb  dominated.  For  the  throat  section,  the 
peak  ebb  velocity  exceeded  the  peak  flood  by  a factor  of  1.22.  The  model 
presented  to  account  for  this  difference  explains  the  ebb  dominance  as  a result 
of  the  different  efficiency  of  water  exchange  between  the  ocean  and  the  bay  at 
high  and  low  tide. 

The  longshore  currents  off  Debidue  Island  were  found  to  be  significantly 
controlled  by  the  wind  stress.  In  a multiple  stepwise  regression  procedure, 
the  longshore  component  .of  wind  velocity  was  found  to  explain  more  o^the 
variance  in  the  observed  longshore  current  velocity  than  any  other  measured 
environmental  parameter. 

Topographic  mapping  of  inlet  shoals  and  adjacent  beaches,  combined  with 
bathymetric  profiling  of  the  throat  and  the  major  channels,  suggest  that  there 
is  a sediment  exchange  between  the  channels  and  the  beaches.  During  periods  of 
•fair  weather,  the  beaches  accrete  and  the  channels  appear  to  scour.  During 
high-energy  conditions,  the  reverse  seems  to  occur. 
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PREFACE 


1.  The  Corps  of  Engineers,  through  its  Civil  Works  program,  has  spon- 
sored, over  the  past  23  years,  research  ii  to  the  behavior  and  characteris- 
tics of  tidal  inlets.  The  Corps'  interest  in  tidal  inlet  research  stems 
from  its  responsibilities  for  navigation,  beach  erosion  prevention  and 
control,  and  flood  control.  Tasked  with  the  creation  and  maintenance  of 
navigable  U.S.  waterways,  the  Corps  dredges  millions  of  cubic  yards  of 
material  each  year  from  tidal  inlets  that  connect  the  ocean  with  bays, 
estuaries,  and  lagoons.  Design  and  construction  of  navigation  improve- 
ments to  existing  tidal  inlets  are  an  important  part  of  the  work  of  many 
Corps'  offices.  In  some  cases,  design  and  construction  of  new  inlets  are 
required.  Development  of  information  concerning  the  hydraulic  character- 
istics of  inlets  is  important  not  only  for  navigation  and  inlet  stability, 
but  also  because  inlets,  by  allowing  for  the  ingress  of  storm  surges  and 
egress  of  flood  waters,  play  an  important  role  in  the  flushing  of  bays 
and  lagoons. 

2.  A research  program,  the  General  Investigation  of  Tidal  Inlets 
(GITI),  was  developed  to  provide  quantitative  data  for  use  in  design 

of  inlets  and  inlet  improvements.  It  is  designed  to  meet  the  following 
objectives : 

To  determine  the  effects  of  wave  action,  tidal  flow,  and  related 
phenomena  on  inlet  stability  and  on  the  hydraulic,  geometric,  and 
sedimentary  characteristics  of  tidal  inlets;  to  develop  the  knowl- 
edge necessary  to  design  effective  navigation  improvements,  new 
inlets,  and  sand  transfer  systems  at  existing  tidal  inlets;  to 
evaluate  the  water  transfer  and  flushing  capability  of  tidal 
inlets;  and  to  define  the  processes  controlling  inlet  stability. 

3.  The  GITI  is  divided  into  three  major  study  areas:  (a)  inlet 
classification,  (b)  inlet  hydraulics,  and  (c)  inlet  dynamics. 

a.  Inlet  Classification.  The  objectives  of  the  inlet  classifi- 
cation study  are  to  classify  inlets  according  to  their  geometry,  hydrau- 
lics, and  stability,  and  to  determine  the  relationships  that  exist  among 
the  geometric  and  dynamic  characteristics  and  the  environmental  factors 
that  control  these  characteristics.  The  classification  study  keeps  the 
general  investigation  closely  related  to  real  inlets  and  produces  an 
important  inlet  data  base  useful  in  documenting  the  characteristics  of 
inlets. 

b.  Inlet  Hydraulics.  The  objectives  of  the  inlet  hydraulics 
study  are  to  define  tide-generated  flow  regime  and  water  level  fluctua- 
tions in  the  vicinity  of  coastal  inlets  and  to  develop  techniques  for 
predicting  these  phenomena.  The  inlet  hydraulics  study  is  divided  into 
three  areas:  (1)  idealised  inlet  model  study,  (2)  evaluation  of  state- 
of-the-art  physical  and  numerical  models,  and  (3)  prototype  inlet 


t 


j 


L 


I 


! 


| 

j 

| 

r 


i 


(1)  The  Idealized  Inlet  Model.  The  objectives  of  this 
model  study  are  to  determine  the  effect  of  inlet  configurations  and 
structures  on  discharge,  head  loss  and  velocity  distribution  for  a num- 
ber of  realistic  inlet  shades  and  tide  conditions.  An  initial  set  of 
tests  in  a trapezoidal  inlet  was  conducted  between  1967  and  1970.  How- 
ever, in  order  that  subsequent  inlat  models  are  more  representative  of 
real  inlets,  a number  of  "idealized"\models  representing  various  inlet 
morphological  classes  are  being  developed  and  tested.  The  effects  of 
jetties ^nd  wave  action  on  the  hydraulics  are  included  in  the  study. 

(2)  Evaluation  of  State-of-the-Art  Modeling  Techniques. 

The  objectives  of  this  part  of  the  inlet  hydraulics  study  are  to  deter- 
mine the  usefulness  and  reliability  of  existing  physical  and  numerical 
modeling  techniques  in  predicting  the  hydraulic  characteristics  of  inlet- 
bay  systems,  and  to  determine  whether  simple  tests,  performed  rapidly  and 
economically,  are  useful  in  the  evaluation  of  proposed  inlet  improvements. 
Masonboro  Inlet,  North  Carolina,  was  selected  as  the  prototype  inlet  which 
would  be  used  along  with  hydraulic  and  numerical  models  in  the  evaluation 
of  existing  techniques.  In  September  1969  a complete  set  of  hydraulic 
and  bathymetric  data  was  collected  at  Masonboro  Inlet.  Construction  of 
the  fixed-bed  physical  model  was  initiated  in  1969,  and  extensive  tests 
have  been  performed  since  then.  In  addition,  three  existing  numerical 
models  were  applied  to  predict  the  inlet's  hydraulics.  Extensive  field 
data  were  collected  at  Masonboro  Inlet  in  August  1974  for  use  in  evaluat- 
ing the  capabilities  of  the  physical  and  numerical  models. 

(3)  Prototype  Inlet  Hydraulics.  Field  studies  at  a number 
of  inlets  are  providing  information  on  prototype  inlet-bay  tidal  hydraulic 
relationships  and  the  effects  of  friction,  waves,  tides,  and  inlet  mor- 
phology on  these  relationships. 

a.  Inlet  Dynamics.  TTie  basic  objective  of  the  inlet  dynamics 
study  is  to  investigate  the  interactions  of  tidal  flow,  inlet  configura- 
tion, and  wave  action  at  tidal  inlets  as  a guide  to  improvement  of  inlet 
channels  and  nearby  shore  protection  works.  The  study  is  subdivided 
into  four  specific  areas:  (1)  model  materials  evaluation,  (2)  movable- 
bed  modeling  evaluation,  (3)  reanalysis  of  a previous  inlet  model  study, 
and  (4)  prototype  inlet  studies. 


(1)  Model  Materials  Evaluation.  This  evaluation  was  initi- 
ated in  1969  to  provide  data  on  the  response  of  movable-bed  model  mate- 
rials to  waves  and  flow  to  allow  selection  of  the  optimum  bed  materials 
for  inlet  models. 


(2)  Movable-Bed  Model  Evaluation.  The  objective  of  this 
study  is  to  evaluate  the  state-of-the-art  of  modeling  techniques,  in 
this  case  movable-bed  inlet  modeling.  Since,  in  many  cases,  movable-bed 
modeling  is  the  only  tool  available  for  predicting  the  response  of  an 
inlet  to  improvements,  the  capabilities  and  limitations  of  these  models 
must  be  established. 
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(3)  Reanalysis  of  an  Earlier  Inlet  Model  Study.  In  1975, 
a report  entitled,  "Preliminary  Report:  Laboratory  Study  of  the  Effect 
of  an  Uncontrolled  Inlet  on  the  Adjacent  Beach,"  was  published  by  the 
Beach  Erosion  Board  (now  CERC) . A reanalysis  of  the  original  data  is 
being  performed  to  aid  in  planning  of  additional  GITI  efforts. 

(4)  Prototype  Dynamics.  Field  and  office  studies  of  a 
number  of  inlets  are  providing  information  on  the  effects  of  physical 
forces  and  artificial  improvements  on  inlet  morphology.  Of  particular 
importance  are  studies  to  define  the  mechanisms  of  natural  sand  bypassing 
at  inlets,  the  response  of  inlet  navigation  channels  to  dredging  and 
natural  forces,  and  the  effects  of  inlets  on  adjacent  beaches. 

4.  This  report  presents  results  of  the  second  phase  of  a field  study 
to  define  the  hydraulics  and  dynamics  of  North  Inlet,  South  Carolina,  a 
natural  Atlantic  coast  tidal  inlet.  Detailed  bathymetric  mapping  of  the 
intertidal  and  shallow  subtidal  zones  is  used  to  define  the  seasonal 
morphologic  variability  of  the  inlet  and  adjacent  beaches.  Wave  and 
tidal  data  are  used  to  correlate  observed  bathymetric  changes  with  the 
causative  processes,  and  to  provide  basic  information  on  wave  conditions 
and  inlet  hydraulic  characteristics.  Results  for  the  1974-75  field  year 
were  included  in  Finley  (1976). 
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CONVERSION  FACTORS,  U.S.  CUSTOMARY  TO  METRIC  (SI) 
UNITS  OF  MEASUREMENT 


U.S.  customary  units  of  measurement  used  in  tb’  * report  can  be  converted 
to  metric  (SI)  units  as  follows: 


Multiply 

by 

To  obtain 

inches 

25.4 

millimeters 

2.54 

centimeters 

square  inches 

6.452 

square  centimeters 

cubic  inches 

16.39 

cubic  centimeters 

feet 

30.48 

centimeters 

0.3048 

meters 

square  feet 

0.0929 

square  meters 

cubic  feet 

0.0283 

cubic  meters 

yards 

0.9144 

meters 

square  yards 

0.836 

square  meters 

cubic  yards 

0.7646 

cubic  meters 

miles 

1.6093 

kilometers 

square  miles 

259.0 

hectares 

knots 

1.852 

kilometers  per  hour 

acres 

0.4047 

hectares 

foot-pounds 

1.3558 

newton  meters 

millibars 

1.0197  x 10- 3 

kilograms  per  square  centimeter 

ounces 

28.35 

grams 

pounds 

453.6 

grams 

0.4536 

kilograms 

ton,  long 

1.0160 

metric  tons 

ton,  short 

0.9072 

metric  tons 

degrees  (angle) 

0.01745 

radians 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  Kelvins1 

*To  obtain  Celsius  (C)  temperature  readings  from  Fahrenheit  (F)  readings, 
use  formula:  C * (5/9)  (F  -32), 

To  obtain  Kelvin  (K)  readings,  use  formula:  K * (5/9)  (F  -32)  ♦ 273.15. 
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SYMBOLS  AND  DEFINITIONS 


Ajj  bay  water  surface  area 
Ac  inlet  cross-sectional  area 
a^,  bay  tidal  amplitude 
a0  ocean  tidal  amplitude 

d water  depth 

F tidal  form  number 

f frequency 

g gravitational  acceleration 
Hjj  breaker  height 
H0  deepwater  wave  height 
H1  water  surface  elevation  in  the  bay 
water  surface  elevation  in  the  ocean 

2 

K Kuelegan's  repletion  coefficient 

k breaker  form  parameter 

Lc  length  of  inlet  channel 

Mj  ith  tidal  harmonic  component,  where  Mj  designates  the  diurnal 
component 

m beach  slope 

P tidal  prism 

Pf  wave  energy  flux 

Pjs  longshore  component  of  the  wave  energy  flux 
Qs  annual  longshore  sediment  transport  rate 
Rq  ocean  tidal  range 

r product-moment  correlation  coefficient 
r'  friction  coefficient 

T semidiurnal  tidal  period  (12.42  hours) 

Tq  deepwater  wave  period 

u mean  current  velocity  in  inlet 
V longshore  current  velocity 

W discharge  into  the  bay  from  sources  other  than  inlet 
breaker  angle 
p density  of  water 

o bed  shear  stress 


HYDRAULICS  AND  DYNAMICS  OF  NORTH  INLET,  SOUTH  CAROLINA,  1975-76 

by 

Dag  Nummedal  and  Stanley  M.  Humphries 
I . INTRODUCTION 


1.  General. 

Tidal  hydraulics,  wave  processes,  and  morphologic  changes  at  North 
Inlet,  South  Carolina,  were  monitored  from  July  1974  to  May  1976,  to  de- 
fine the  physical  environmental  parameters  characterizing  this  stable, 
unmodified,  tidal  inlet  on  the  mesotidal  coast  of  the  southeastern  United 
States.  Results  from  the  first  year  of  study,  including  the  geological 
setting,  the  littoral  process  variability,  the  hydraulic  parameters  of 
the  inlet,  and  detailed  analysis  of  longshore  wave  energy  flux,  sediment 
transport  characteristics,  and  the  development  of  the  inlet  and  associ- 
ated shoals  over  historic  time,  are  presented  in  Finley  (1976). 

The  present  report  completes  the  presentation  of  hydraulic,  litto- 
ral, and  morphologic  data  for  the  2-year  study  period.  To  complement  the 
work  of  Finley  (1976) , this  report  focuses  on  a more  detailed  analysis 
of  the  travel  of  the  tidal  wave  through  the  inlet,  the  relationship  be- 
tween seasonal  changes  in  mean  sea  level  (MSL)  and  the  morphologic  re- 
sponse, and  the  question  of  the  mechanics  of  longshore  current  genera- 
tion. 

North  Inlet  is  located  in  Georgetown  County,  South  Carolina,  and  pro- 
vides drainage  to  the  Atlantic  Ocean  for  about  24  by  106  square  meters 
of  marsh  behind  Debidue  and  North  Islands  (Fig.  1).  The  marsh  is  bound- 
ed to  the  south  by  Winyah  Bay,  a large  estuary  formed  at  the  confluence 
of  the  Pee  Dee  and  Waccamaw  Rivers.  To  the  west  or  north,  the  marsh 
borders  directly  on  Pleistocene  units  of  the  lower  Coastal  Plain  phys- 
iographic province.  Further  details  on  local  stratigraphy  are  provided 
in  Finley  ( 1976) . 

2.  Processes . 

The  hydraulic  behavior  of  a mesotidal  inlet  like  North  Inlet  is  a 
function  of  the  adjacent  ocean  tidal  range,  the  wave  climate,  the  size 
of  the  bay,  and  the  degree  of  development  of  marsh  within  the  drainage 
area.  To  facilitate  an  evaluation  of  how  representative  the  findings 
at  North  Inlet  would  be  for  other  southeast  inlets  of  similar  size,  the 
regional  process  variability  between  Cape  Hatteras  and  south  Florida 
was  reviewed.  North  Inlet  was  found  to  occupy  a transitional  position 
between  wave-dominated  and  tide-dominated  coastal  segments.  Therefore, 
process  characterization  of  the  inlet  system  entailed  documenting  the 
mean  values  and  the  ranges  of  all  pertinent  variables  describing  the 
nearshore  wave  climate  and  the  inlet  currents. 
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The  hydraulics  questions  asked  in  this  study  address  the  problems 
■of  time  asymmetry  and  velocity  asymmetry  in  the  inlet  currents,  the  ef- 
fects of  the  propagation  of  the  tidal  wave  through  shallow  water  on  its 
harmonic  constituents,  the  relation  between  water  surface  differentials 
and  tidal  creek  current  velocities,  inlet  throat  equilibrium,  and  the 
annual  variations  in  MSL  as  recorded  by  the  tide  gage  network.  Because 
the  wave  process  measurements  were  designed  primarily  to  document  en- 
vironmental conditions  rather  than  performing  specific  studies  of  the 
surf  zone,  the  analysis  of  this  data  set  was  limited  to  one  question  of 
fundamental  importance-the  mechanics  of  the  generation  of  longshore 
currents.  The  existence  of  an  extensive  data  set  at  one  location  for 
a wide  range  of  conditions  provided  the  basis  for  a statistical  evalua- 
tion of  the  observed  current  velocities.  Stepwise  regression  procedures 
were  used  to  test  how  the  dependent  variable  (the  observed  current 
velocity)  correlates  with  a series  of  simultaneously  measured  indepen- 
dent variables. 

In  addition  to  characterizing  the  dynamics  of  the  North  Inlet  sys- 
tem, the  parameters  provide  a data  base  for  testing  a number  of  basic 
scientific  hypotheses. 

3.  Response. 

The  monitoring  of  morphologic  response  over  the  study  period  in- 
cluded recording  channel  cross-sectional  profiles,  beach  and  nearshore 
profiles,  topographic  mapping  of  selected  intertidal  features,  and 
quarterly  oblique  aerial  photography.  Basic  documentation  of  the  vari- 
ability was  the  primary  objective,  and  all  pertinent  morphologic  data 
have  therefore  been  included,  either  as  text  or  in  data  Appendixes  A 
to  I. 

The  basic  scientific  question  addressed  in  the  analysis  of  the  mor- 
phologic data  pertained  to  the  seasonal  variability  in  selected  inter- 
tidal features  and  their  significant,  and  apparently  cyclic,  pattern 
of  change  in  response  to  the  fluctuations  in  MSL.  Seasonal  sea  level 
fluctuation  apparently  has  a much  more  profound  effect  on  intertidal 
sediment  distribution  than  has  been  recognized.  The  monthly  surveys 
of  beach  profiles  and  channel  cross-sectional  profiles  provide  a data 
base  for  testing  Escoffier's  (1940)  basic  stability  concept:  If  inlet 
stability  is  maintained  by  increased  tidal  scour  subsequent  to  an  in- 
filling of  the  main  channel  by  increased  wave  action,  then  channel  pro- 
files and  beach  profiles  adjacent  to  the  inlet  should  be  out  of  phase 
with  respect  to  erosion  and  accretion.  The  North  Inlet  morphologic  data 
appear  to  support  this  stability  concept. 

II.  REGIONAL  VARIABILITY 


1.  Introduction. 

Processes  and  morphology  at  North  Inlet  reflect  its  location  rela- 
tive to  the  global  atmospheric  circulation  pattern,  the  tides  of  the 
Atlantic  Ocean  and  their  transformation  by  the  southeast  Atlantic 
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Continental  Shelf,  and  the  physiography  of  the  Coastal  Plain.  On  this  re- 
gional scale.  North  Inlet  occupies  an  interesting  transitional  location: 
to  the  north,  the  Outer  Banks  are  clearly  wave  dominated;  to  the  south, 
the  shores  of  southern  South  Carolina  and  Georgia  are  clearly  tide  domi- 
nated. Physiographically,  the  coastal  zone  north  of  North  Inlet  is 
characterized  by  either  straight  beaches  eroding  the  Pleistocene  main- 
land (e.g..  Myrtle  Beach;  Brown,  1975)  or  narrow  Holocene  barriers 
fringing  a system  of  wide  shallow  lagoons  (e.g.,  the  Outer  Banks  and 
Pamlico  Sound).  To  the  south,  the  South  Carolina  and  Georgia  coasts 
have  stubby  Holocene  and  Pleistocene  barrier  islands  backed  by  exten- 
sive intertidal  marsh  and  tidal  creeks  (Hoyt,  1967;  Hubbard  and  Barwis, 
1976) . These  marsh-filled  bays  create  a hydraulic  geometry  similar  to 
that  studied  by  Oliveira  (1970)  for  bays  with  sloping  banks,  and  tend 
to  develop  ebb-dominated  tidal  inlets. 

2.  Tidal  Range. 

Based  on  the  morphology  of  sandy  shorelines  of  the  world,  Davies 
(1964)  recognized  that  tidal  range  is  the  single  most  diagnostic 
variable.  He  proposed  the  following  class  boundaries:  microtidal 
(spring  tidal  range  0 to  2 meters) , mesotidal  (tidal  range  2 to  4 
meters) , and  macrotidal  (tidal  range  greater  than  4 meters) . Accord- 
ingly, the  coast  of  the  southeastern  United  States  would  be  mesotidal 
only  in  southern  South  Carolina  and  Georgia,  and  microtidal  at  North 
Inlet  and  farther  north  (Table  1;  Fig.  2).  By  defining  mesotidal  es- 
tuaries as  those  where  sedimentary  forms  due  to  deposition  by  tidal 
currents  begin  to  dominate  those  built  by  waves  (Hayes,  1975),  it  is 
found  that  the  lower  boundary  of  the  mesotidal  range  in  the  southeast- 
ern United  States  should  be  adjusted  downward,  perhaps  as  far  as  to 
1.6  meters.  The  North  Inlet  morphology  (spring  tidal  range  1.8  meters) 
clearly  fits  the  mesotidal  models  proposed  by  Oertel  (1972)  and  Hayes, 
et  al.  (1973).  The  low  wave  energy  of  the  southeast  United  States 
coast,  compared  to  the  world  average  (Davies,  1973),  is  probably  the 
reason  behind  the  development  of  mesotidal  morphology  at  a lower 
tidal  range. 

The  observed  increase  in  tidal  range  southward  from  the  Outer 
Banks  to  Georgia  reflects  the  effect  of  shoaling  of  the  tidal  wave 
across  the  increasingly  wider  Continental  Shelf.  Silvester  (1974) 
reported  a fair  correlation  between  shelf  width  and  tidal  range  for  the 
margins  of  all  major  world  oceans.  Figure  3 is  a scatter  plot  illus- 
trating the  offshore  distance  of  the  200-meter  depth  contour  versus 
the  corresponding  open- coast  tidal  range.  As  the  tidal  wave  travels 
onto  the  United  States  southeast  shelf,  the  cotidal  lines  intersect 
the  bathymetric  contours  at  a relatively  steep  angle  (Redfield,  1958). 
Therefore,  the  "effective"  shelf  width  of  any  coastal  tide  station  may 
be  larger  than  the  values  given  in  Figure  3,  accounting  for  some  of 
the  scatter.  The  coastal  tidal  range  depends  on  the  amplitude  of  the 
deepwater  ocean  tide,  the  actual  shelf  profile,  and  the  influence  of 
local  coastal  embayments.  Therefore,  the  relationship  derived  from 
Figure  3 does  not  necessarily  pertain  to  all  shelves. 
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Table  1.  Spring  and  mean  tidal  ranges  along  the  southeastern 


coast  of  the  United  States  (from  National  Oceanic 
and  Atmospheric  Administration,  1977) . 


Ocean  station 

Spring 

Mean 

(m) 

(ft) 

(m) 

(ft? 

Kitty  Hawk 

1.16 

3.8 

0.93 

3.2 

Cape  Hatteras 

1.25 

4.1 

1.04 

3.4 

Cape  Lookout 

1.34 

4.4 

1.13 

3.7 

Atlantic  Beach 

1.31 

4.3 

1.10 

3.6 

Myrtle  Beach 

1.83 

6.0 

1.55 

5.1 

Winyah  Bay  entrance 

1.65 

5.4 

1.40 

4.6 

Cape  Romain 

1.68 

5.5 

1.43 

4.7 

Isle  of  Palms 

1.86 

6.1 

1.58 

5.2 

Charleston  Harbor 
entrance 

1.86 

6.1 

1.58 

5.2 

Microtidal1 

Edisto  Beach 

2.10 

6.9 

1.80 

5.9 

Fripp  Inlet 

2.23 

7.3 

1.89 

6.2 

Mesotidal1 

Hilton  Head 

2.38 

7.8 

2.01 

6.6 

Tybee  Light 

2.44 

8.0 

2.07 

6.8 

According  to  classification  by  Davies  (1964) . 
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Figure  3.  Observed  tidal  range  versus  offshore  distance  to  the 

200-meter  bathymetric  contour  for  selected  stations  of 
the  eastern  United  States.  All  data  refer  to  open-coast 
stations.  Circled  data  points  refer  to  the  identified 
tide  stations  (from  National  Oceanic  and  Atmospheric 
Administration,  1977). 
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3.  Wave  Energy. 


Wave  energies  on  the  southeast  shelf  of  the  United  States  are  moder- 
ate on  a world  scale  (Davies,  1973).  Because  the  intensity  of  breaking 
waves  directly  controls  the  rate  of  longshore  sediment  transportation  and 
significantly  affects  the  pattern  of  tidal  delta  sand  dispersal,  regional 
variations  in  wave  climate  influence  inlet  morphological  development. 
Since  wave  energy  flux  generally  predicts  the  rate  of  coarse-grained 
sediment  transport,  deepwater  wave  energy  flux  values  were  calculated 
for  the  southeast  shelf  (Table  2);  data  were  obtained  from  Summary  of 
Synoptic  Meteorological  Observations  (SSMO)  (U.S.  Naval  Weather  Service 
Command,  1975)  for  the  time  period  1963  to  1971.  The  SSMO  data  squares 
are  outlined  in  Figure  4 together  with  the  wave  energy  flux  diagrams. 


Table  2.  Deepwater  wave  energy  flux  values  (in  units  of 
1010  ergs  per  meter  per  second)  for  the  south- 
eastern United  States. 


N. 

NE. 

E. 

SE. 

S. 

SW. 

w. 

NW. 

Cape  Hatteras 

4.8 

3.3 

1.3 

1.2 

1.9 

2.8 

3.0 

2.8 

Charleston 

3.2 

3.5 

1.9 

1.3 

2.2 

2.5 

3.2 

2.6 

Jacksonville 

4.2 

2.2 

1.3 

1.1 

1.8 

1.6 

1.9 

3.2 

Miami 

1.9 

2.2 

2.4 

1.2 

00 

o 

00 

O 

0.7 

1.3 

These  calculations  show  a significant  decrease  in  total  wave  energy 
and  a general  shift  in  the  dominant  flux  directions  toward  the  south. 

The  energy  flux  directions  are  given  as  those  from  which  the  waves  are 
coming  and  are,  therefore,  consistent  with  the  common  format  for  wind 
diagrams.  Wave  energy  flux  is  calculated  from  the  equation: 

Pf  = 9.58  x 109  x Hq  x T0  (1) 

where  H0  is  deepwater  wave  height  (meters) , To  is  wave  period  (seconds) , 
and  Pf  is  the  wave  energy  flux  expressed  in  ergs  per  meter  per  second. 

An  in-depth  discussion  of  the  procedures,  assumptions,  and  errors  in- 
volved in  the  calculations  of  wave  energy  flux  values  from  SSMO  data 
is  presented  in  Walton  (1973)  and  Nummedal  and  Stephen  (1976) . 

4.  Inlet  Size. 

Tidal  range,  coastal  plain  geomorphology,  and  wave  energy  variations 
have  caused  a distinct  regional  trend  in  inlet  size.  This  size  varia- 
tion is  expressed  in  Figure  5 in  terms  of  tidal  prisms,  cross  sections, 
and  hydraulic  radii  of  the  inlet  throats.  The  data  are  taken  from 
Jarrett  (1976),  Finley  (1976)  for  North  Inlet,  and  FitzGerald,  Num- 
medal and  Attaway  (1977)  for  Price  Inlet.  The  significance  of  tidal 
range  as  an  overriding  factor  controlling  inlet  size  is  evident  when 
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WAVE  ENERGY  FLUX 


NORFOLK 


Figure  4.  Wave  energy  flux  diagrams  for  the  southeastern  United 
States  (energy  flux  calculated  from  data  in  U.S.  Naval 
Weather  Service  Command,  1975) . 
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Figure  5.  Size  variation  in  southeastern  tidal  inlets.  Inlets 
with  hydraulic  radii  larger  than  that  at  Price  Inlet 
have  only  ebb  tidal  deltas  (solid  bars);  those  with 
smaller  radii  have  both  ebb  and  flood  tidal  deltas 
(open  bars) . 
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comparing  the  distributions  in  Figures  2 and  5.  In  all  measures  of 
size,  the  inlets  reach  a maximum  in  southern  South  Carolina  and  Georgia 
and  then  decrease  rapidly  into  Florida  and,  to  a lesser  extent,  into 
North  Carolina.  Inlet  shoal  development  shows  an  interesting  relation- 
ship to  the  depth  of  the  main  channel.  All  inlets  listed  in  Figure  5, 
which  are  at  least  as  deep  as  Price  Inlet  (hydraulic  radius  equal  to 
5.2  meters),  appear  only  to  have  developed  ebb  tidal  deltas  (outer 
shoals) ; the  inlets  shallower  than  Price  have  developed  both  ebb  tidal 
and  flood  tidal  deltas  (inner  shoals) . 

5.  Storm  Frequencies. 

Since  the  initiation  of  geological  studies  at  North  Inlet  in  1972, 
only  two  storm  events  have  caused  significant  morphological  changes  at, 
and  adjacent  to,  the  inlet.  The  storm  of  10  and  11  February  1973,  which 
was  classified  as  a moderate  to  severe  storm  (U.S.  Army,  Corps  of  Engi- 
neers, Coastal  Engineering  Research  Center,  1977),  caused  erosion  up  to 
31.3  cubic  meters  per  linear  meter  of  beach  (Finley,  1976).  Another 
storm  on  23  and  24  September  1974,  classified  as  light  to  moderate, 
caused  up  to  13.0  cubic  meters  of  erosion  per  linear  meter  (Kana,  1977). 
Hurricane  Hazel , by  comparison,  eroded  up  to  52  cubic  meters  per  meter 
at  Myrtle  Beach  in  1954  (U.S.  Congress,  1966). 

These  observations  indicate  that  despite  intensive  sediment  trans- 
portation and  redeposition  during  storm  events,  such  events  may  be  too 
rare  to  significantly  affect  the  long-term  morphologic  development  of 
a system  like  North  Inlet.  As  pointed  out  by  Davies  (1973),  the  South 
Carolina  coast  is  south  of  the  American  east  coast  storm  environment. 

Fox  and  Davis  (1976)  also  reported  that  most  extratropical  cyclones 
("northeast"  storms)  pass  to  the  north  of  the  South  Carolina  coast. 

Tropical  storms  and  hurricanes  do  occasionally  hit  the  South  Caro- 
lina coast;  however,  the  frequency  is  substantially  less  than  that  at 
the  Outer  Banks  (Fig.  6).  According  to  Meyers  (1975),  the  annual  prob- 
ability of  hurricane  force  winds  (windspeed  in  excess  of  120  kilometers 
per  hour)  striking  within  an  80-kilometer  segment  of  coast  in  northern 
South  Carolina  is  about  0.050.  This  corresponds  to  the  one  hurricane 
landfall  every  20  years. 

It  was  assumed  at  the  initiation  of  this  study  that  2 weeks  of 
measurements  every  season  would  give  representative  samples  of  the 
annual  population  of  process  parameters  controlling  inlet  morphologic 
response.  The  low  frequency  of  cyclones  (tropical  and  extratropical) 
crossing  the  shoreline  near  North  Inlet  appears  to  justify  that  as- 
sumption. 

III.  TIDAL  HYDRAULICS 

1.  Introduction. 

One  fundamental  objective  of  the  North  Inlet  project  was  to  document 
the  hydraulic  conditions  of  this  unmodified  tidal  inlet.  Three  tide 


Figure  6.  Annual  probability  of  hurricane  landfalls  within 
80-kilometer  segments  on  the  southeast  coast  of 
the  United  States  (from  Meyers,  1975). 
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gages  were  operated  throughout  the  study  period  from  July  1974  to  May 
1976.  Leupold  and  Stevens  float  gages  in  30-centimeter  culverts  sup- 
ported by  5-centimeter  water  pipes  jetted  into  the  channel  bed  were 
installed  at  the  intersection  of  Town  and  Debidue  Creeks  and  at  Jones 
Creek  (Fig.  7).  A Bristol  bubbler  gage  with  the  orifice  just  seaward 
of  North  Island  (Fig.  7)  was  used  to  record  tidal  fluctuations  in  the 
ocean. 

To  complement  the  tide  gaging  program,  current  velocities  were 
monitored  with  an  ENDECO  directional  current  meter  for  a number  of 
13-hour  periods  at  stations  located  at  the  inlet  throat  and  at  cross 
sections  of  Town  and  Jones  Creeks.  Changes  in  the  hydraulic  geometry 
of  the  inlet  system  over  time  were  measured  by  repetitive  longitudinal 
and  transverse  channel  surveys  with  a Bludworth  fathometer. 

Finley  (1976)  presented  a comprehensive  analysis  of  the  first 
year's  tidal  hydraulics  data  in  terms  of  (a)  the  relationship  between 
inlet  throat  cross-sectional  area  and  the  tidal  prism,  (b)  the  varia- 
tions in  tidal  prism  with  tidal  phase,  (c)  the  relation  between  ocean- 
bay  water  surface  differentials  and  current  velocities  for  hydro- 
graphic  measurement  periods,  (d)  variation  in  channel  friction  co- 
efficient (Manning's  n)  over  the  same  tidal  cycles,  and  (e)  the 
efficiency  of  tidal  exchange  through  the  inlet  expressed  by  Keulegan's 
(1967)  repletion  coefficient. 

The  analysis  of  the  second  year's  tidal  hydraulics  data  emphasizes 
other  aspects  of  the  hydraulics  as  well  as  some  of  Finley's  points  in 
order  to  achieve  as  complete  a characterization  of  inlet  behavior  as 
possible.  Computer  printouts  of  a large  number  of  both  raw  and  reduced 
hydraulic  data  are  given  in  Appendixes  A to  G. 


Hydraulic  Geometry 


North  Inlet  drains  about  23.7  x 106  square  meters  of  marsh,  chan- 
nels, and  shallow  bays  (based  on  planimetry  of  U.S.  Department  of  Agri- 
culture vertical  air  photos).  The  drainage  boundary  between  the  southern 
part  of  the  North  Inlet  system  and  the  tributaries  to  Winyah  Bay  is 
poorly  defined  and  fluctuates  somewhat,  depending  on  the  stage  of  the 
rivers  draining  into  the  bay,  the  tidal  phase,  and  seasonal  fluctuations 
in  sea  level.  These  factors  may  introduce  about  5-percent  variability 
into  the  above  figure. 

The  ocean  tidal  wave  is  distorted  by  passage  through  the  tidal  inlet 
(Boon,  1973;  Byrne  and  Boon,  1976;  Nuamedal,  FitzGerald,  and  Humphries, 
1976)  due  to  the  presence  of  water  level  dependent  impedance  factors 
like  changing  channel  cross  section,  variations  in  bay  surface  area, 
and  bottom  friction.  As  demonstrated  by  Brown  (1928)  and  Keulegan 
(1967),  the  ratio  between  the  bay  surface  area,  Aj,,  and  the  inlet 
throat  cross-sectional  area,  significantly  affects  the  efficiency 
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Figure  7.  Location  of  the  three  tide  gages  relative  to  the  main 
channels  and  barrier  islands.  Mn  refers  to  the  mean 
annual  tidal  ranges  determined  at  the  gage  stations. 


of  water  exchange  between  the  ocean  and  the  marsh  (bay) . Values  for  A5 
were  obtained  by  planimetry  of  the  creek  surface  area  on  one  vertical 
air  photo  (CDW  2K)  taken  near  mean  low  water  (MLW)  and  another  (USDA 
A40)  obtained  near  MSL.  Unfortunately,  A^  could  not  be  determined  for 
spring  low  water.  For  any  water  stage  above  MSL,  significant  parts 
of  the  marsh  are  inundated,  and  A^  cannot  be  estimated  except  at  spring 
high  tide  when  the  vegetation  change  at  the  high  water  line  can  be  used 
to  trace  the  entire  marsh  outline.  Ac  values  were  directly  obtained 
from  fathometer  records  at  the  inlet  throat,  calibrated  by  the  surveyed 
distance  between  two  transect  end  markers.  To  relate  Aj,  and  Ac  to  the 
same  tide  stages,  spring  high  was  taken  as  90  centimeters  above  MSL, 
and  "MLW  was  set  at  70  centimeters  below  MSL.  Table  3 summarizes  the 
results. 


Table  3.  Bay  water  surface  area,  Aj,,  and 
inlet  throat  flow  area,  Ac,  as 
functions  of  water  level. 


Water  level 

Ab 

(106  m2) 

(m2) 

Spring  high  water 

23.66 

1,680 

Mean  sea  level 

7.33 

1,475 

Mean  low  water 

6.32 

1,330 

Ab  varies  by  a factor  of  3.74  through  the  tidal  cycle;  Ac  only  changes 
by  about  1.26  for  the  same  water  level  range.  Therefore,  the  ratio 
Ab/Ao  which  controls  inlet  impedance  according  to  Keulegan's  (1967) 
analysis,  is  2.9  times  larger  at  spring  high  than  at  MLW.  These  varia- 
tions appear  to  be  typical  for  South  Carolina  inlets.  FitzGerald, 
Nummedal,  and  Kana  (1977)  documented  the  Ajj/Ac  ratio  for  Price  Inlet 
as  about  2.7  times  larger  at  spring  high  than  at  MLW. 


These  conditions  are  at  variance  with  the  geometry  of  the  Wacha- 
preague  Inlet  system  in  Virginia  (Byrne,  Bullock,  and  Tyler,  1975), 
which  includes  large  open  bays  that  reduce  Ab  variations. 


A complicating  factor  which  is  not  pursued  any  further  in  this  study 
is  the  time  of  travel  of  the  tidal  wave  in  the  shallow  creek  system. 

The  distance  from  the  North  Inlet  throat  section  to  the  most  distant 
head  of  a tidal  creek  is  about  6 kilometers,  with  a "typical"  channel 
depth  at  low  water  of  2 meters.  The  celerity  of  a free  wave  in  such 
a channel  is  4.43  meters  per  second,  giving  a maximum  traveltime  of  22 
minutes  for  the  tidal  wave  to  reach  the  perimeter  of  the  tidal  system. 
Water  level  adjustments  on  the  shallow  intci-channel  marsh  surface  are 
substantially  slower.  Therefore,  at  any  tide  stage,  there  are  water 
surface  slopes  within  the  tidal  creek-marsh  system  which  affect  the  hy- 
draulics at  the  measured  sections  in  an  unknown  way. 
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The  rapid  changes  in  channel  hydraulic  geometry  required  repetitive 
surveys  of  designated  cross  sections  on.  a monthly  basis.  Figure  8 shows 
the  locations  of  the  monitored  sections  as  well  as  the  "typical"  profile 
for  the  1975-76  field  year.  All  fathometer  profiles  have  been  tied  to 
a common  reference  datum  used  for  the  tide  gages  by  means  of  the  gage 
record  obtained  at  the  time  of  each  survey.  Appendix  A contains  pro- 
files obtained  for  (a)  the  inlet  throat,  (b)  Town  Creek,  (c)  Jones 
Creek,  and  (d)  the  flood  tidal  delta,  from  Jurie  1975  to  May  1976. 


The  inlet  throat  profile  shows  a distinct  separation  between  the 
nearly  8-meter-deep  main  ebb  channel  to  the  south  and  a wide  intertidal 
to  shallow  subtidal  platform  made  up  of  the  channel-margin  linear  bar 
and  marginal  flood  channel  to  the  north.  Only  a small  percentage  of  the 
tidal  prism  passes  over  the  intertidal  platform.  However,  large  varia- 
tions in  throat  cross-section  area  are  accounted  for  by  depth  changes  on 
this  platform.  The  total  throat  cross  section  varies  by  a factor  of  1.34 
between  a high  of  1,679  square  meters  in  November  and  an  alltime  low  of 
1,255  square  meters  in  May.  This  variation  is  probably  a reflection  of 
the  seasonal  change  in  N6L.  High  MSL  in  November  increases  the  mean 
daily  tidal  prism  and  scours  the  channel  floor.  However,  the  low  MSL  in 
May  reduces  the  tidal  prism.  The  throat  section,  in  turn,  adjusts  by  a 
decrease  in  size  in  accordance  with  Escoffier’s  (1940)  stability  concept. 


The  cross-sectional  areas  of  Town  and  Jones  Creeks  (Table  4)  appear 
to  undergo  seasonal  variations  in  phase  with  the  throat  section.  De- 
Alteris  and  Byrne  (1975)  recorded  a seasonal  pattern  of  variation  in 
inlet  flow  area  consistent  with  the  one  detected  at  North  Inlet. 


Table  4.  Monthly  channel  cross-sectional  areas  below  MSL 
(in  square  meters) . 


Date 

Town  Creek 

Jones  Creek 

Inlet 

Throat 

June  1975 
July 

1,091 

1,087 

(+4)1 
(4  3) 

490  (-2) 

520  (44) 

1,585 

(+7) 

Aug. 

1,034 

(-2) 

470  (-6) 

Sept . 

1,124 

(+7) 

483  (-4) 

Oct. 

1,238 

(415) 

510  (42) 

1,445 

(-2) 

Nov. 

Dec. 

1,056 

(45) 

500  (0) 

1,679 

(+12) 

862 

(-18) 

(+11) 

407  (-19 

548  (49) 

Jan.  1976 

Feb. 

1,177 

1,298 

(-12) 

Mar. 

863 

(-18) 

515  (43) 

1,676 

(+12) 

Apr. 

1,033 

(-2) 

593  (416) 

1,388 

(-6) 

May 

956 

(-9) 

473  (-6) 

1,255 

(+15) 

Mean 

1,051 

SOI 

1,475 

Percent  change  from  the  mean. 
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3.  Tidal  Currents. 


a.  Theory.  Keulegan  (1967),  King  (1974),  and  Seelig,  Harris,  and 
Herchenroder  (1977)  have  presented  theoretical  models  which  permit  the 
derivation  of  inlet  current  histories  given  the  ocean  tide  and  the  hy- 
draulic geometry  of  the  system  as  inputs.  Keulegan' s model  and  its 
applicability  to  the  North  Inlet  system  was  discussed  by  Finley  (1976) . 
The  models  by  King  and  Shemdin  (1975)  and  Seelig,  Harris,  and  Herchen- 
roder (1977)  include  factors  not  represented  in  Keulegan 's  model:  iner- 
tial effects,  discharge  into  the  bay  from  sources  other  than  the  tidal 
inlet,  and  variations  in  bay  and  throat  section  areas. 

The  equations  governing  flow  from  the  ocean  through  a constricted 
inlet  into  a bay  are  the  equations  of  motion  and  of  continuity.  These 
can  be  formulated  as: 

g(H2  " Hl)  = k ^ + 7 (1  + ^ k)  u|u|  (2) 

and 


u A„ 


dH 

+ w - dT 


(3) 


where  u is  the  mean  velocity  at  the  inlet,  Lc  is  the  entrance  length, 
r'  is  the  square  law  friction  coefficient  given  by  r'  = a/p  x u2.  H2 
and  Hi  are  the  ocean  and  bay  surface  elevations,  W is  discharge  into 
the  bay  by  sources  other  than  the  inlet,  d is  the  water  depth  in  the 
inlet,  p is  the  density  of  water,  and  o is  the  bed  shear  stress. 


The  equations  show  that  (a)  the  water  surface  elevation  differ- 
ential is  proportional  to  the  square  of  the  mean  current  velocity,  (b) 
that  the  flood  duration  is  different  from  the  ebb  duration  when  Ac  and 
Ab  vary  over  the  tidal  cycle,  (c)  that  peak  (and  mean)  ebb  velocity  may 
be  different  from  peak  (and  mean)  flood  velocity,  and  (d)  that  the  bay 
tidal  range  generally  is  different  from  that  of  the  ocean  (Nummedal, 
FitzGerald,  and  Humphries,  1976). 


This  section  discusses  the  North  Inlet  hydraulic  data  in  light 
of  the  hydraulic  principles  outlined  above. 
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b.  Data.  Tidal  current  velocity  data  were  obtained  during  23  tidal 
cycles  between  June  1975  and  May  1976.  Locations  of  the  individual 
current  measurement  stations  are  shown  in  Figure  9.  The  distribution 
of  observations  per  station  is  summarized  in  Table  5.  Each  current 
velocity  was  calculated  as  the  arithmetic  mean  of  individual  readings 
obtained  from  a vertical  array  of  current  meter  positions  spaced  1 
meter  apart.  As  discussed  by  Finley  (1976),  the  vertical  current  velo- 
city variation  was  invariably  low,  because  density  stratification 
generally  is  absent,  and  bottom  drag  only  affects  the  near-bottom  flow 
layer.  For  this  analysis,  water  surface  differentials  were  directly 
measured  off  the  raw  tide  records  and  hand-plotted.  The  gages  had 
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Location  of  hydrography  stations  at  North  Inlet 
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been  surveyed  to  a common  reference  datum  (Finley,  1976).  The  discon- 
tinuous tide  records  made  it  impossible  to  derive  water  surface  differ- 
ential time  series  for  more  than  half  of  the  data  sets  (Table  5) . 

Sample  plots  of  current  velocities  versus  the  corresponding  water 
surface  differentials  are  shown  in  Figure  10;  the  complete  data  set  is 
contained  in  Appendix  B. 

The  water  surface  differential  plotted  in  the  Jones  Creek  figures 
refers  to  that  between  the  ocean  and  Jones  Creek  gages.  The  differen- 
tial tised  in  the  Town  Creek  and  inlet  throat  figures  refers  to  that  be- 
tween the  ocean  and  Town  Creek  gages. 

c.  Analysis . Figure  10  and  Appendix  B demonstrate  that  the  current 
velocity  curves  are  asymmetric  with  higher  ebb  than  flood  velocities. 
They  are  also  asymmetric  in  time,  with  peak  ebb  velocity  occurring  a- 
bout  halfway  between  ocean  high  and  low  water;  peak  flood  velocity  oc- 
curs shortly  before  ocean  high.  The  flood,  consequently,  lasts  longer 
than  the  ebb.  Examples  of  similar  time-velocity  asymmetries  are  known 
from  many  tidal  systems  (Klein,  1970;  Hayes,  et  al.,  1973;  Boothroyd 
and  Hubbard,  1975;  Hine,  1975)  and  reflect  the  presence  of  differential 
impedance  to  inlet  flow  at  high  and  low  water. 


A number  of  environmental  factors  affect  each  set  of  tidal  hydro- 
graphic  observations.  Therefore,  reliable  quantification  of  hydraulic 
process  parameters  must  be  based  on  a statistical  analysis  of  a rep- 
resentative sample. 

For  the  statistical  analysis,  the  following  variables  were  defined 
from  Appendix  C:  VFMAX,  peak  flood  velocity  in  a given  tidal  cycle 
(centimeters  per  second);  VEMAX,  peak  ebb  velocity  (centimeters  per 
second);  DELTF,  time  (in  hours)  between  peak  flood  velocity  and  subse- 
quent high  water  in  the  ocean;  DELTE,  time  between  high  water  in  the 
ocean  and  subsequent  peak  ebb  velocity;  LAGH,  time  between  ocean  high 
tide  and  subsequent  station  slack;  LAGL,  time  between  ocean  low  tide  and 
subsequent  station  slack.  Data  sets  of  100  for  the  1975-76  season  (App. 
C)  and  45  for  the  1974-75  season  were  obtained.  Summary  statistics  are 
presented  in  Table  6.  Both  years,  peak  ebb  velocities  exceeded  peak 
flood.  For  all  velocity  stations  combined,  the  peak  ebb  velocity  in 
1975-76  was  82.8  centimeters  per  second  which  was  compared  to  70.7 
centimeters  per  second  for  the  flood,  an  average  peak  ebb  velocity 
1.17  times  the  average  peak  flood.  In  1974-75,  the  peak  ebb  velocity 
exceeded  peak  flood  by  as  much  as  1.26.  Peak  velocities  during  in- 
dividual tidal  cycles  covered  a wide  range.  Over  the  2 years  of  study, 
peak  ebb  velocities  ranged  from  30  to  156  centimeters  per  second; 
peak  flood  velocities  ranged  from  24  to  128  centimeters  per  second. 

The  ebb  dominance  of  the  inlet  circulation  is  even  more  pronounced 
when  the  individual  velocity  sections  are  considered  separately 
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Table  6.  Summary  statistics  on  North  Inlet  current  parameters 
and  hydraulic  lag  constants,  1974  to  1976. 


Variable1 

N 

Mean 

Median 

Std.  dev. 

Min. 

Max. 

1974-75 

VFMAX 

45 

70.60 

70 

22.94 

26 

128 

VEMAX 

40 

88.43 

84 

18. 10 

59 

156 

DELTF 

42 

1.55 

1.6 

0.65 

0.2 

2.9 

DELTE 

32 

3.45 

3.3 

0.45 

2.5 

4.4 

UGH 

35 

0.86 

0.75 

0.39 

0.1 

2.0 

UGL 

20 

0.87 

0.65 

0.54 

0.3 

2.1 

1975-76 

VFMAX 

92 

70.74 

67 

22.60 

24.0 

125 

VEMAX 

91 

82.80 

80 

25.57 

30.0 

150 

DELTF 

92 

1.33 

1.25 

C .62 

0.0 

3.0 

DELTE 

82 

3.06 

3.05 

0.72 

0.9 

4.8 

UGH 

89 

0.81 

0.8 

0.37 

0.0 

1.9 

UGL 

58 

0.52 

0.45 

0.32 

0.0 

1.6 

Barnes  defined  in  text,  symbols  listing,  and  Appendix  C. 


(Table  7).  For  the  throat  section,  the  peak  ebb  velocity  exceeded  the 
peak  flood  by  a factor  of  1.32  in  1974-75  and  1.22  in  1975-76. 

This  degree  of  ebb  current  dominance  over  flood  currents  through  the 
inlet  throat  is  comparable  to  Price  Inlet  where  FitzGerald,  Nummedal,  and 
Attaway  (1977)  document  an  average  peak  ebb  velocity  16  percent  larger 
than  the  average  flood  over  a data  set  of  12  complete  tidal  cycles.  How- 
ever, New  Corpus  Christi  Pass,  connecting  the  microtidal  Gulf  of  Mexico 
to  the  open  Corpus  Christi  Bay,  is  dominated  by  floodtides.  The  average 
peak  flood  velocity,  over  eight  tidal  cycles  measured  in  1972-73,  exceeds 
the  ebb  by  as  much  as  33  percent  (Behrens,  Watson,  and  Mason,  1977). 

These  current  asymmetries  are  consistent  with  equations  (2)  and  (3) 
when  Aj,  and  Ac  vary  during  the  tidal  cycle  in  accordance  with  local 
hydraulic  geometries. 

Currents  at  the  Town  Creek  and  Jones  Creek  velocity  sections  are 
weaker  than  those  at  the  throat  section  by  about  83  and  76  percent, 
respectively,  for  peak  flood,  and  by  79  and  70  percent,  respectively,  for 
peak  ebb.  Calculations  indicate  that  about  7 percent  of  the  mean  tidal 
prism  (including  flow  up  Sixty  Bass  Creek)  is  stored  between  the  throat 
and  the  two  velocity  sections,  and  that  the  combined  flow  areas  of  Town 
Creek  and  Jones  Creek  are  5 percent  larger  than  the  cross  section  at  the 
throat  (Table  5).  These  two  factors  appear  adequate  to  explain  the 
reduced  current  velocities  in  the  tidal  creeks. 
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Table  7.  Summary  statistics,  by  channel  section,  on  North  Inlet 
current  parameters  and  hydraulic  lag  constants. 


Variable1 

1974 

-75 

1975- 

-76 

N 

Mean 

Std.  dev. 

_JL_j 

Mean 

Std.  dev. 

Inlet  throat  (ITN, 

ITS) 

(FCN,  ITN,  ITS) 

VFMAX 

6 

85.13 

33.36 

25 

83.05 

24.89 

VEMAX 

6 

112.00 

26.37 

24 

101.05 

23.48 

DELTF 

6 

1.58 

0.60 

25 

1.63 

0.71 

DELTE 

6 

3.68 

0.58 

20 

2.94 

0.59 

LAGH 

6 

1.02 

0.39 

22 

0.80 

0.32 

LAGL 

4 

0.68 

0.53 

14 

0.34 

0.15 

Town  Creek  (TCN1,  TCN2,  TCN3) 

(TCN1,  TCN2,  TCN3) 

VFMAX 

24 

72.74 

18.83 

37 

66.79 

20.55 

VEMAX 

21 

87.86 

11.78 

37 

80.57 

25.44 

DELTF 

21 

1.70 

0.62 

37 

1.34 

0.52 

DELTE 

15 

3.42 

0.42 

37 

3.05 

0.69 

LAGH 

18 

0.79 

0.26 

41 

0.85 

0.39 

LAGL 

8 

0.69 

0.30 

28 

0.59 

0.36 

Jones 

Creek  (JCN1,  JCN2) 

(JCN1,  JCN2) 

VFMAX 

15 

61.33 

20.58 

30 

65.57 

19.96 

VEMAX 

13 

78.46 

6.74 

30 

70.17 

0.61 

DELTF 

15 

1.28 

0.70 

28 

1.10 

0.69 

DELTE 

11 

3.35 

0.46 

30 

2.84 

0.35 

UGH 

11 

0.86 

0.76 

25 

0.74 

0.31 

UGL 

8 

1.09 

0.74 

17 

0.53 

^Names  defined  in  text,  symbols  listing,  and  Appendix  C. 
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Frequency  distributions  for  variables  VFMAX  and  VEMAX  (Apps.  C and 
D)  illustrate  a pronounced  unimodal  distribution,  consistent  with  semi- 
diurnal tides  with  only  minor  diurnal  inequality. 

Variables  DELTF  and  DELTE  express  the  current  time  asymmetry  with 
respect  to  the  ocean  tidal  cycle.  This  time  asymmetry  shows  no  statis- 
tically significant  variation  between  stations  or  years  (Table  7).  Gen- 
erally, peak  flood  occurs  1.44  hours  before  ocean  high,  and  peak  ebb 
occurs  3.26  hours  after  ocean  high  (Table  6).  The  frequency  histograms 
for  these  variables  (Apps.  C and  D)  are  all  strongly  skewed.  Random  fac- 
tors like  meteorologic  tides,  measurement  errors,  and  resolution  limits 
produce  a normal  distribution.  The  existence  of  a sharp  cutoff  on  the 
right  side  of  the  distribution  in  three  out  of  four  histograms  reflects 
the  tendency  of  the  inlet  flow  system  to  establish  equilibrium  by  in- 
creasing surface  slope,  and,  consequently,  current  velocity  in  response 
to  a delay  in  adjusting  the  bay  surface  to  the  elevation  of  the  ocean. 

The  tidal  lags  (i.e.,  the  times  between  ocean  extremes  and  station 
slacks)  were  also  unequal  (Table  6).  High  water  lag  averaged  0.81 
hour;  low  water  lag  averaged  only  0.52  hour  (mean  lag  values  for  the 
1974-75  season  are  not  included  here  because  the  standard  deviations 
of  the  distributions  are  as  large  as  the  means) . The  frequency  histo- 
grams (Fig.  11)  illustrate  the  difference  between  high  and  low  water 
lags  even  better  than  the  means.  The  low  water  lag  mode  for  1975-76 
was  at  zero,  illustrating  that  there  is  a larger  number  of  measured 
tidal  cycles  with  no  low  water  lag  than  any  other  class.  The  high 
water  mode,  on  the  contrary,  was  at  1 hour  for  1975-76  and  at  0.86  hour 
for  1974-75.  Equations  (2)  and  (3)  predict  less  flow  impedance  (or 
greater  exchange  efficiency)  at  low  values  of  A(j/Ac  than  at  high;  i.e., 
the  bay  water  level  adjusted  much  more  easily  to  that  of  the  ocean  at 
low  water  than  at  high.  This  agrees  with  the  observations. 

Different  lags  for  high  and  low  water  were  also  observed  at  Govern- 
ment Cut,  Florida  (O'Brien  and  Dean,  1972).  At  high  tide,  the  lag  was 
520  or  1.79  hours.  At  low  tide,  it  was  60°  or  2.07  hours.  The  fact 
that  low  tide  lag  exceeded  high  tide  lag  was  probably  due  to  the  hy- 
draulic geometry  of  the  system.  Government  Cut  links  the  ocean  to  the 
large  Miami  Harbor  and  Biscayne  Bay  rather  than  a marsh-tidal  creek 
system. 

It  is  apparent  at  this  point  that  knowing  the  average  lag  and  the 
average  Keulegan  repletion  coefficient  for  a tidal  cycle  is  not  ade- 
quate. Net  sediment  discharge  and,  consequently,  the  morphological 
development  of  the  entire  inlet  system,  are  controlled  by  the  flow 
asymmetry  which  is  related  to  the  "variation"  in  the  lag  and  repletion 
coefficient  over  time.  Table  8 summarizes  the  high  tide  and  low  tide 
median  lag  values  for  North  Inlet  for  the  1974-75  seasons  (median 
values  represent  the  central  tendency  of  skewed  distributions  like 
those  in  Fig.  11  better  than  the  mean).  K is  consistently  larger 
at  low  water  than  at  high,  demonstrating  a more  efficient  water  ex- 
change. 
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Figure  11.  Histograms  summarizing  the  distribution  of 


high  water  and  low  water  tidal  lag  at  the 
North  Inlet  throat  section. 
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Table  8.  Median  phase  lags  measured  at  the  throat  section  at  North 


Inlet;  predicted  Keulegan  repletion  coefficients  (K)  and 
relative  bay  tidal  amplitude  (aj,/a0)  (K  and  ajj/a0  from 
O'Brien  and  Dean,  1972). 


Season 

Lag 

Tidal  stage 

Hours 

Degrees 

K 

ab/a0 

1974-75 

Low 

0.65 

18.6 

1.32 

0.95 

High 

0.77 

22.3 

1.22 

0.92 

1975-76 

Low 

0.45 

13.1 

1.60 

0.98 

High 

0.78 

22.6 

1.22 

0.92 

The  predicted  tidal  amplitudes  compare  fairly  well  to  those 
actually  observed.  The  average  of  the  Jones  and  Town  Creek  amplitudes 
relative  to  that  in  the  ocean  is  0.95  for  the  1974-75  data  (Finley, 

1976)  and  0.92  for  the  1975-76  data. 

To  summarize  the  discussion  of  tidal  lags  and  time  velocity  asym- 
metry, Figure  12  illustrates  the  "average"  conditions  at  North  Inlet. 

This  tidal  cycle  has  a period  of  12.4  hours,  and  the  ocean  tide  is  at 
MSL  at  time  t = 0 hour.  At  t * 1.66  hours,  a peak  flood  velocity  of  84 
centimeters  per  second  is  attained;  1.44  hours  later,  at  t = 3.1  hours, 
the  ocean  tide  is  at  its  high  with  bay  high  being  attained  0.77  hour 
later  or  at  t = 3.87  hours.  At  3.26  hours  after  ocean  high,  a peak  ebb 
current  of  107  centimeters  per  second  occurs,  and  2.94  hours  after  that, 
at  t = 9.3  hours,  the  ocean  reaches  low  tide.  However,  the  bay  continues 
to  ebb  to  reach  low  water  0.55  hour  later;  i.e.,  at  t = 9.85  hours.  At 
t = 12.4  hours,  the  ocean  surface  returns  to  MSL.  For  this  tidal  cycle, 
the  inlet  currents  flood  for  a total  of  6.42  hours  and  ebb  for  5.98 
hours,  a difference  of  0.44  hour  or  26  minutes.  This  compares  to  an 
average  duration  difference  of  42  minutes  at  Price  Inlet  (FitzGerald, 
Nummedal,  and  Attaway,  1977)  and  26  minutes  at  Wachapreaque  Dock  in 
Virginia  (Boon,  1973). 

To  test  for  the  correlation  between  inlet  current  velocities  and 
water  surface  differentials,  hourly  velocity  readings  at  each  individual 
station  and  the  corresponding  surface  differentials  (App.  E)  were  ana- 
lyzed by  a simple  correlation  routine.  Table  9 is  a simplified  output 
matrix;  a complete  correlation  matrix  and  some  scattergrams  are  included 
in  Appendix  E.  Variable  names  DELTOW  and  DELJ0N  refer  to  water  surface 
differentials  at  Town  Creek  and  Jones  Creek,  respectively.  DIFT  and  DIFJ 
refer  to  the  absolute  values  of  same.  VELT,  VEU,  VELM  and  AVELT,  AVELJ, 
AVELM  refer  to  the  mean  velocities  and  the  absolute  values  of  the  mean 
velocities  at  Town  Creek,  Jones  Creek,  and  the  inlet  throat,  respectively. 

The  linear  correlations  between  the  mean  velocities  of  any  section 
and  the  corresponding  water  surface  differentials  are  fair.  For  example, 
the  correlation  is  0.78  between  the  instantaneous  mean  velocity  at  the 
Town  Creek  section  (VELT)  and  the  water  surface  differential  between 
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Table  9.  Correlation  matrices  between  water  surface 

differential  values  and  current  velocities,  1975-76. 1 


QELTOW 

DEUON 

VELT 

VEU 

VELM 

DELTOW* 

1 

190 

91 

87 

64 

DEUON 

0.81 

1 

93 

89 

69 

VELT 

0.78 

0.75 

1 

66 

11 

VEU 

0.79 

0.83 

0.97 

1 

22 

VELM 

0.86 

0.88 

0.95 

0.98 

1 

DIFT 

DIFJ 

AVELT 

AVEU 

DIPT 

1 

190 

91 

87 

DIFJ 

0.C3 

1 

93 

89 

AVELT 

0.54 

0.38 

1 

66 

A VEU 

0.50 

0.53 

0.90 

1 

*Nuabers  below  the  diagonal  give  the  Pearson  correlation 
coefficients.  Numbers  above  give  the  nuaber  of  data  pairs 
used  in  the  same  calculations. 

2Nwaes  are  defined  in  text,  symbol  listing,  and  Appendix  E. 


the  ocean  and  Town  Creek  gages  (DELTOW) . Similar  correlation  coeffi- 
cients characterize  the  other  linear  regressions.  However,  the  correla- 
tion coefficients  decrease  when  the  absolute  values  of  the  variables  are 
considered.  The  correlations  between  the  water  surface  differentials 
and  the  square  of  the  velocities  are  also  poor.  Equation  (2)  indicates 
that  the  square  of  the  current  velocity  is  proportional  to  the  water 
surface  differential.  However,  the  fact  that  both  variable  ranges  in 
these  correlations  are  reduced  to  about  one-half  by  using  the  absolute 
(or  square)  values  gives  the  data  an  apparent  increase  in  scatter.  The 
good  linear  correlation  does  support  equation  (2);  failure  to  obtain 
improved  correlation  by  using  the  velocity  squared  does  not  contradict 
it. 

4.  Tidal  Curves. 

All  tide  records  were  digitized  and  hourly  water  surface  elevation 
values  were  output  on  cards  and  verified.  Tidal  curves  were  plotted 
from  the  punched  records.  All  output  was  plotted  on  a monthly  basis 
regardless  of  the  duration  of  the  record.  Sample  outputs  are  shown  in 
Figure  13,  and  the  complete  set  of  plotted  curves  is  in  Appendix  F. 

All  tidal  curves  are  referenced  to  a common  datum  of  205  centimeters 
below  the  "INLET"  triangulation  station.  The  tidal  curves  are  not  con- 
tinuous because  of  logistics  problems.  The  availability  of  tide  re- 
cords for  each  gage  is  summarized  in  Appendix  F. 
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The  curves  illustrate  well  the  semidiurnal  tide,  the  variations  in 
diurnal  inequality,  and  the  tidal  range  fluctuations.  For  the  ocean 
gage,  the  diurnal  inequality  at  tropical  tides  can  reach  55  centimeters 
with  most  of  the  inequality  generally  due  to  unequal  high  water  levels. 
At  equatorial  tides,  the  inequality  is  all  but  nonexistent.  The  daily 
mean  tidal  range  varies  by  as  much  as  60  centimeters  from  an  average 
spring  range  of  175  centimeters  to  a neap  range  of  only  117  centimeters. 
The  mean  tidal  range  and  MSL  were  calculated  for  selected  2-week  periods. 
The  mean  annual  tidal  ranges  were  computed,  in  turn,  as  the  arithmetic 
mean  of  the  2-week  ranges.  Tidal  ranges  for  the  1975-76  field  year  and 
1974-75  (from  Finley,  1976)  are: 


1975-76 

1974- 

(cm) 

(cm) 

Ocean  gage 

167 

151 

Town  Creek 

140 

145 

Jones  Creek 

135 

140 

5.  Water  Surface  Differentials. 

Scarcity  of  overlapping  tide  records  limits  the  number  of  differ- 
ential curves.  Computer-plotted  differential  curves  for  the  period  28 
April  to  12  May  1976  (Fig.  14)  demonstrate  good  correspondence  with  the 
hand-plotted  differentials  in  Appendix  B.  Typically,  peak  differentials 
ranged  from  about  10  to  20  centimeters  at  neap  tides  (1  to  4 May)  to 
about  30  centimeters  at  spring  (12  May).  The  jaggedness  of  the  differ- 
ential curves  is  probably  attributed  to  random  errors.  Normally,  most 
irregularities  occur  near  slack  water  (0  differential)  at  a time  when 
both  measurement  error  and  weather-induced  level  fluctuations  are  most 
noticeable.  Figure  14  demonstrates  that  throughout  most  of  the  time 
period,  peak  ebb  differentials  seem  to  equal  those  for  peak  flood.  At 
springtides,  however,  the  seaward  slopes  are  clearly  the  dominant  ones. 
This  indicates  that  the  ebb  flow  dominance  and  the  consequent  net  sea- 
ward sediment  transport  are  to  a large  extent  caused  by  the  spring  tides. 

6.  Spectral  Analysis. 

A spectral  analysis  was  made  of  the  water  level  time-series  data 
for  a number  of  2-week  periods  of  overlapping  tide  records  for  the 
ocean  gage  and  at  least  one  of  the  two  bay  gages  (Table  10). 


The  first  six  time  periods  (Table  10)  of  records  contained  360 
hourly  water  level  observations.  The  absolute  value  of  the  calculated 
spectral  density  function  was  compared  from  one  record  to  another,  and 
an  increase  or  a decrease  in  the  relative  contribution  of  a given  har- 
monic component  was  determined. 

The  procedure  adjusts  the  mean  of  the  sample  to  zero  and  then 
creates  a raw  periodogram  and  a spectral  density  function  as  measures 
of  the  variance  explained  by  a simple  harmonic  oscillation  of  a given 
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Table  10.  Water  level  time-series  data  subjected  to  spectral  analysis.1 


Date 

Ocean  gage 

Town  Creek 

Jones  Creek 

9 

to  23  July 

1975 

X 

- 

X 

26 

Oct.  to  10 

Nov. 

1975 

X 

X 

X 

8 

to  23  Dec. 

1975 

X 

- 

X 

2 

to  17  Jan. 

1976 

- 

X 

X 

14 

to  29  Mar. 

1976 

X 

X 

X 

1 

to  15  May 

1976 

X 

X 

X 

26 

Oct.  to  31 

Dec. 

1975 

- 

- 

X 

Spectra  are  presented  in  Figure  15  and  Appendix  G. 

frequency.  The  cutoff  frequency  is  that  corresponding  to  twice  the 
observation  interval.  In  this  analysis,  T cutoff  equals  2 hours. 
Since 


f cutoff  = ir  . 

Calculations  were  made  for  180  frequency  values  less  than  it;  the  fre- 
quency interval  between  each  successive  spectral  density  estimate 
was  0.0175  cycle  per  hour.  To  facilitate  reading  Figure  15  and  the 
spectra  in  Appendix  G,  Table  11  presents  the  relation  between  fre- 
quency and  period  calculated  according  to  equation  (4) . 

Table  11.  Radian  frequency,  f,  in  cycles  per  hour  versus  period,  T. 


f 3.14 

2.09 

1.57  1.26  1.05  0.89 

0.78 

0.70 

0.63  0.57  0.52 

T 2 

3 

4 5 6 7 

8 

9 

10  11  12 

The  ordinate  of  the  spectral  density  diagrams  is  expressed  in  units 
of  meters  squared  times  hours  or  (2ir/f)  meters  squared. 

• The  distinctly  semidiurnal  tidal  characteristics  at  North  Inlet 

are  expressed  by  the  dominance  of  the  0.5  cycle  per  hour  peak  which 
corresponds  to  a period  of  12.4  hours.  Typically,  the  spectral  den- 
sity value  for  this  semidiurnal  peak  is  about  50  times  larger  than 
that  for  the  diurnal  peak  (Table  12).  The  semidiurnal  solar  constit- 
uent, Sy,  and  the  lunar  constituent,  Mj,  are  too  similar  in  fre- 
quency to  be  resolved  by  the  bandwidths  used;  the  same  applies  to  the 
diurnal  components.  Si  and  Mi.  However,  it  is  not  necessary  for 
these  components  to  be  resolved.  An  approximate  tidal  form  number 
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Figure  15.  Spectral  density  functions  for  2-week  records 
(1  to  15  May  1976)  off  the  ocean.  Town  Creek, 
and  Jones  Creek  gages.  Spectral  density 
(ordinate)  is  expressed  in  log- 10  scales. 
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Table  12.  Spectral  densities  for  tidal  harmonic  components  Mj , M2,  M4, 
Mg,  and  Mg  at  ocean  gage.  Town  Creek,  and  Jones  Creek,  26 
October  to  10  November  1975. 


Station 

Mi 

m2 

m4 

Me 

m8 

Ocean  gage 

0.08 

5.5 

0.004 

0.0005 

0.0003 

Town  Creek 

0.07 

3.6 

0.01 

0.003 

0.0003 

Jones  Creek 

0.08 

3.6 

0.01 

0.003 

0.0004 

M!/M2 

m4/m2 

m6/m2 

m8/m2 

Ocean  gage 

1.5  x 

10‘2 

7.3  x 10"4 

9.1  x 10'5 

5.5  x 10"5 

Town  Creek 

1.9  x 

10'2 

2.8  x 10"3 

8.3  x 10"4 

8.3  x 10'5 

Jones  Creek 

2.2  x 

10'2 

2.8  x 10’3 

8.3  x 10‘4 

1.1  x 10"4 

(Defant,  1958)  can  be  calculated  if  it  is  assumed  that  the  Ml  and  M2 
values  of  Table  12  equal  the  sum  of  the  corresponding  lunar  and  solar 
component  values.  The  mean  values  of  the  form  number  F,  defined  as 


(5) 


are  F (ocean)  equals  0.0145,  F (Jones  Creek)  equals  0.0194,  and  F 
(Town  Creek)  equals  0.0222.  All  demonstrate  a strong  semidiurnal  form. 

Figure  15  and  Appendix  G demonstrate  the  presence  of  shallow- 
water  tidal  components  of  higher  than  semidiurnal  frequency.  These 
include  the  M4,  M^,  and  Mg  components,  the  overtides , labeled  on 
all  spectra,  as  well  as  other  peaks  possibly  due  to  interactions  be- 
tween the  major  tidal  constituents.  Table  12  summarizes  the  magni- 
tude of  the  overtide  spectral  densities  at  the  three  gage  sites  for 
a "typical"  2-week  period,  26  October  to  10  November  1975.  The  high 
frequency  overtides  have  spectral  densities  ranging  from  about  3 x 
10  ^ to  5 x 10  5 times  that  of  the  semidiurnal  component.  Therefore, 
the  significance  of  these  overtides  is  due  more  to  their  amplifica- 
tion with  the  passage  of  the  tidal  wave  through  the  inlet,  than  their 
absolute  magnitudes.  Table  12  demonstrates  that  the  amount  of  vari- 
ance in  the  tide  record,  accounted  for  by  the  M4  overtide,  increases 
by  a factor  of  4 as  the  tidal  wave  progresses  from  the  ocean  to  the 
bay  gages.  The  other  overtides  considered,  M(j  and  Mg,  also  appear  to 
account  for  more  of  the  variance  in  the  bay  spectra  than  in  the  ocean 
spectrum.  Since  the  tidal  wave  had  traveled  across  a wide  shallow 
swash  platform  before  it  was  recorded  by  the  ocean  gage  (Fig.  7),  sig- 
nificant shallow-water  components  were  already  present  in  this  record. 
The  total  amplification  of  the  overtides  from  deep  ocean  water  to  bay 
conditions  was  even  larger  than  indicated  by  these  figures. 
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Low-frequency  harmonics  (periods  more  than  12  hours)  do  not  change 
much  due  to  tidal  wave  progress  through  the  inlet  (Fig.  15  and  Table 
12).  Therefore,  the  tidal  inlet  acts  as  a low-pass  filter  for  tidal 
wave  travel  from  the  ocean  to  the  bay-marsh  system.  The  amplification 
of  the  overtides  suggests  that  inlet  impedance  is  a time-varying  factor 
introducing  higher  order  harmonics,  through  a series  of  nonlinear  terms 
in  the  flow  equations  (eqs.  2 and  3).  The  results  of  the  spectral 
analysis  are  conceptually  consistent  with  the  hydraulic  data  presented 
in  this  section  and  Figure  11.  The  amplification  of  overtide  ampli- 
tudes in  shallow  water  has  also  been  found  at  Wachapreague  (Boon, 

1973)  and  at  Chesterfield  Inlet,  Northwest  Territories  (Budgell,  1976). 

7.  Inlet  Equilibrium. 

Empirical  studies  of  Pacific  coast  inlets  by  O'Brien  (1931)  demon- 
strated that  the  cross-sectional  area  of  the  inlet  throat  is  propor- 
tional to  the  tidal  prism,  according  to  the  equation: 

Ac  = 4.69  x 1CT4  P°-85  # (6) 

where  Ac  is  the  minimum  flow  section  of  the  inlet  (throat)  measured  be- 
low MSL  ( in  square  feet) , and  P is  the  diurnal  or  spring  tidal  range  (in 
cubic  feet).  O'Brien  (1969)  and  Johnson  (1973)  introduced  slightly 
different  equations  expressing  the  same  equilibrium  concept. 

Escoffier  (1940)  and  Inman  and  Frautschy  (1966)  discussed  the  equi- 
librium concept  in  detail,  stating  that  a stable  inlet  cross  section 
is  maintained  by  a balance  between  the  amount  of  sediment  scoured  by 
tidal  currents  and  the  amount  supplied  as  wave- generated  longshore 
drift.  Therefore,  an  inlet  equilibrium  equation  should  include  both 
tidal  prism  and  wave  energy  as  independent  variables.  A development 
in  this  direction  has  been  provided  by  Jarrett  (1976)  who  investigated 
108  tidal  inlets  along  the  coast  of  the  United  States.  Jarrett  con- 
cluded that  natural  or  single-jettied  inlets  on  the  Atlantic  coast 
have  larger  cross  sections  for  a given  tidal  prism  than  on  the  Pacific 
coast,  presumably  because  of  the  lower  wave  energy  on  the  east  coast. 

Using  the  same  symbols  and  units  as  in  equation  (6),  Jarrett's  equa- 
tion for  east  coast  natural  inlets  is: 

A - 5.37  x 10‘6  P1-07  . (7) 

Based  on  channel  hydrography  measurements  and  throat  section  fathometer  » 

profiles,  the  flood,  ebb,  and  mean  daily  tidal  prisms  were  calculated 
(Table  13)  for  each  data  set  between  May  1975  and  May  1976.  Throat 
cross  sections  were  precisely  measured  during  11  hydrography  periods. 

These  data  sets  and  Finley's  (1976)  data  for  the  previous  year  were 
plotted  on  Jarrett's  (1976)  equilibrium  diagram  (Fig.  16). 

All  data  points  plot  well  within  the  95-percent  confidence  limit  on 
Jarrett's  empirical  diagram,  indicating  that  North  Inlet  is  at  a state  of 
long-term  stability.  The  scatter  of  the  data  points  reflects  measurement 
errors,  meteorologic  tides,  tidal  exchange  between  the  North  Inlet  system 
and  Winyah  Bay,  and  significant  annual  variations  in  sea  level. 
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Table  13.  North 

Inlet  i 

ebb  and  flood 

tidal  discharge 

• 

Prism 

Date 

Tide 

Flood 

Ebb 

Mean 

Mean 

Millions  of 

m3 

108ft3 

1975 

24  May 

Spring 

24 . 066 

15.724 

19.895 

7.025 

10  June 

Spring 

29.088 

21.773 

25.430 

8.979 

14  June 

Mean 

19.894 

21.341 

20.617 

7.280 

17  June 

Neap 

18.464 

19.123 

18.794 

6.636 

21  June1 

Mean 

15.631 

12.820 

14.225 

5.023 

29  June 

Mean 

12.942 

18.918 

15.930 

5.625 

26  Aug. 

Mean 

13.658 

15.311 

14.490 

5.116 

26  Sept. 

Mean 

14.630 

16.092 

15.360 

5.424 

29  Sept. 

Neap 

18.446 

18.706 

18.580 

6.561 

2 Oct. 1 

Mean 

24.926 

23.022 

23.980 

8.467 

4 Oct. 

Spring 

19.336 

24.008 

21.672 

7.652 

24  Nov. 

Mean 

18.630 

21.186 

19.908 

7.029 

1976 

7 Jan. 1 

Mean 

11.700 

10.800 

11.250 

3.972 

10  Jan. 1 

Neap 

8.071 

8.586 

8.329 

2.941 

15  Jan. 

Mean 

18.515 

23.274 

20.894 

7.378 

17  Jan. 

Spring 

19.555 

29.530 

24.543 

8.666 

17  Feb.1 

Spring 

20.534 

21.614 

21.070 

7.439 

15  Mar.1 

Spring 

21.395 

18.047 

19.721 

6.964 

22  Mar.1 

Mean 

10.393 

7.844 

9.119 

3.220 

24  Mar. 

Neap 

8.550 

8.233 

8.392 

2.963 

27  Mar.1 

Mean 

11.963 

14.854 

13.410 

4.735 

27  Apr. 1 

Mean 

15.732 

10.145 

12.940 

4.569 

19  May1 

Mean 

13.104 

12.880 

12.992 

4.588 

^ased  on  throat  section  stations. 
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Figure  17  presents  a correlation  diagram  of  10  measurements  of  tid- 
al prism  and  the  ocean  tidal  range,. Rq,  using  both  the  data  in  Table 
13  and  the  recorded  tidal  ranges.  The  regression  equation  is: 

p = 12R0  - 5.10  (8) 

where  P is  measured  in  106  cubic  meters  and  Rq  in  meters.  The  Pearson 
product  moment  correlation  coefficient,  r,  equals  0.94.  This  correla- 
tion clearly  demonstrates  that  tidal  prism  is  a direct  function  of  the 
available  storage  volume  in  the  marsh-creek  system.  A good  linear  cor- 
relation for  Price  Inlet  was  demonstrated  by  FitzGerald,  Nummedal,  and 
Kana  (1977)  by  averaging  hourly  water  level  readings. 

8.  Seasonal  Changes  in  Sea  Level. 

MSL  (relative  to  a reference  datum)  was  calculated  for  2-week  cy- 
cles from  the  available  ocean,  Town  Creek,  and  Jones  Creek  gage  re- 
cords. The  resulting  MSL  curve  is  shown  in  Figure  18.  A seasonal 
pattern  of  high  sea  level  in  the  fall,  an  alltime  low  in  mid-winter, 
and  a secondary  high  in  mid-spring  is  quite  evident.  The  measured 
annual  range  at  the  creek  gages  is  about  34  centimeters,  the  infer- 
red ocean  range  is  about  45  centimeters.  Annual  sea  level  changes 
significantly  affect  the  sedimentation  and  patterns  of  morphologic 
change  of  the  intertidal  zone  (Sec.  V,  4).  To  assess  the  applicability 
of  certain  North  Inlet  morphologic  results  to  inlets  in  general,  the 
geographic  variations  in  and  reasons  for  the  sea  level  fluctuations 
will  be  discussed  in  some  detail. 

Pronounced  seasonal  variations  in  sea  level  are  found  at  nearly 
all  coastal  locations  (Marmer,  1952;  Pattullo,  et  al.,  1955;  Pattullo, 
1966),  although  the  pattern  of  variation  changes  from  place  to  place. 

For  Charleston  Harbor,  Pattullo,  et  al.  (1955)  reported  an  annual 
high  in  October  of  17  centimeters  above  MSL  and  a low  in  January  of 
8.7  centimeters  below  MSL,  based  on  tide  records  from  1921  to  1946. 

There  was  no  indication  of  a secondary  spring  high  in  that  record. 

In  the  Gulf  of  Mexico,  an  interesting  increase  in  the  secondary  spring 
high  toward  the  west  was  observed  (Marmer,  1952;  Pattullo,  et  al. 

1955).  At  Key  West,  no  spring  high  was  observed.  At  Port  Isabel 
(Marmer,  1952),  and  at  the  New  Corpus  Christi  Pass  (Behrens,  Watson, 
and  Mason,  1977),  a pronounced  spring  high  was  observed.  The  sea 
level  curve  in  Behrens,  Watson,  and  Mason  (1977)  corresponds  closely 
to  the  one  for  North  Inlet  (Fig.  18). 

A number  of  factors  can  contribute  to  annual  changes  in  sea  level. 
Pattullo  (1966)  listed  (a)  variations  in  water  temperature,  (b)  varia- 
tions in  salinity,  and  (c)  changes  in  local  atmospheric  pressure; 

Meade  and  Emery  (1971)  listed  river  runoff  in  the  eastern  United 
States.  Each  individual  factor  can  cause  seasonal  deviations  in  sea 
level  of  about  20  centimeters.  Behrens,  Watson,  and  Mason  (1977) 
qualitatively  attribute  the  changes  at  Corpus  Christi  to  seasonal 
variations  in  onshore  and  longshore  winds,  factors  which  Pattullo 
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Mean  Tidal  Prism(106m3) 


Figure  17.  Measured  tidal  prism  versus  ocean  tidal  range 
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(1966)  was  unable  to  assess  quantitatively.  An  additional  set  of  fac- 
tors known  to  have  a significant  influence  on  sea  level  is  the  upwell- 
ing  and  geostrophic  effects  of  ocean  currents  (Wyrtki,  1973;  Komar, 
1976) . 

The  sea  level  curve  observed  at  North  Inlet  for  1975-76  is  clearly 
the  result  of  a multitude  of  factors.  The  steric  effect  (Pattullo, 
et  al.,  1955)  explains  an  annual  range  of  about  25  centimeters,  with  a 
single  maximum  in  October  and  a broad  minimum  in  the  winter  and  spring 
Therefore,  the  steric  effect  is  probably  the  dominant  factor.  North- 
east winds  in  spring  and  fall  could  have  contributed  to  the  secondary 
peak  in  mid-spring  and  the  increase  in  fall  sea  surface  elevation  be- 
yond the  rise  due  to  the  steric  effect.  Offshore  winds,  often  associ- 
ated with  winter  frontal  passage,  may  be  partly  responsible  for  the 
extreme  January  low. 

Figure  19  illustrates  that  tidal  range  variations  also  follow  a 
seasonal  cycle.  However,  this  cycle  is  much  more  pronounced  in  the 
creeks  than  at  the  ocean  gage.  This  implies  that  the  astronomical 
and  meteorological  tidal  forces  have  less  control  on  the  relative  bay 
tidal  range  than  does  the  seasonal  change  in  inlet  efficiency  which 
is  dependent  on  MSL.  The  pronounced  low  in  both  Town  and  Jones  Creek 
tidal  ranges  relative  to  the  ocean  range  in  winter  is  in  qualitative 
agreement  with  the  arguments  in  Section  III,  which  implies  higher  in- 
let impedance,  a smaller  Keulegan  repletion  coefficient,  and  a re- 
duced relative  bay  tidal  range  for  a larger  value  of  the  A^/Ac.  Low- 
ering the  MSL  in  winter  reduces  tidal  prism,  throat  cross  section 
(Table  4),  and  the  Keulegan  repletion  coefficient. 

IV.  LITTORAL  PROCESSES 


1'.  Introduction. 

The  morphology  of  a tidal  inlet  and  associated  shoals  is  controlled 
directly  by  the  waves  through  their  interaction  with  the  tidal  currents 
to  produce  swash  bars  and  channel- margin  linear  bars  and,  indirectly, 
by  the  wave-generated  longshore  drift  which  provides  a source  for  net 
sediment  supply  to  the  system.  In  this  study,  littoral  process  varia- 
tions were  monitored  during  every  2-week  seasonal  study  to  establish 
a basis  for  the  calculation  of  longshore  sediment  transport  rates, 
and  to  evaluate  changes  in  wave  climate  over  the  2-year  study  period 
which  might  be  correlated  to  observed  variations  in  sedimentation  pat- 
terns. 

An  extensive  data  set  of  longshore  current  velocities,  breaker 
parameters,  and  windspeeds  was  collected  during  the  year,  because  of 
the  unresolved  question  of  the  relative  importance  of  the  momentum 
of  the  breaking  wave  and  the  stress  of  the  longshore  wind  in  generat- 
ing the  observed  nearshore  current  (Shepard  and  Inman,  1950;  Harrison 
and  Krumbein,  1964). 
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Figure  19.  Annual  variation  in  tidal  ranges  at  the  three  stations. 


2.  Wind  Patterns. 


Annual  offshore  wind  data  (Finley,  1976)  demonstrate  that  South 
Carolina  coastal  winds  dominantly  come  from  the  northeast  and  the 
southwest.  The  former  component  is  largely  due  to  extratropical  cy- 
clones, most  of  which  are  far  to  the  northeast  of  South  Carolina  but 
control  the  wind  pattern  along  the  entire  east  coast  of  the  United 
States.  South  and  southwesterly  winds  are  generally  anticyclonic  sum 
mer  winds  associated  with  a strong  Bermuda  high. 

Wind  velocities  and  directions  were  measured  by  a Weather  Measure 
W121  wind  indicator  on  an  hourly  basis  at  station  DB1-10.  Wind  velo- 
cities for  1975-76  agree  with  the  1974-75  values  (Finley,  1976).  The 
directional  data  for  the  2 years  were  summarised  in  Figure  20.  The 
seasonal  pattern  shows  a fair  degree  of  consistency.  In  winter  (Janu 
ary),  winds  from  the  north  and  northeast  dominate,  reflecting  frontal 
passage  over  the  study  area.  Strong  northerly  winds,  bringing  cold 
and  generally  dry  air  into  the  southeast,  typically  succeed  such 
fronts.  March  has  northerly  winds  due  to  frontal  passage  combined 
with  an  equal  frequency  of  southerly  winds  associated  with  strong  and 
generally  migrating  subtropical  high-pressure  systems.  By  mid-summer 
the  subtropical  high  assumes  a stationary  position  over  Bermuda,  and 
persistent  southerly  winds  bring  warm  moist  air  into  the  southeast, 
resulting  in  unstable  atmospheric  conditions;  i.e.,  thunderstorms 
with  strong  winds  and  a dispersed  directional  spectrum.  The  fall 
months  have  the  highest  frequency  of  extratropical  cyclones  off  the 
northeast  coast  of  the  United  States,  producing  almost  exclusively 
winds  from  the  northeast  quadrant. 

Resultant  wind  velocity  vectors  for  each  seasonal  study  period 
from  July  1974  to  March  1976  were  computed  for  station  DBI-10  (Fig. 
21).  The  directions  follow  the  scheme  described  above.  The  result- 
ant wind  vectors  for  the  suirmer  and  fall  were  rather  consistent  for 
the  2 years.  Due  to  an  onshore  (azimuth  71°)  wind  vector  in  January 

1975,  as  compared  to  the  offshore  (azimuth  13°)  vector  in  January 

1976,  the  annual  resultant  wind  vector  for  1974-75  was  much  larger 
than  the  one  for  1975-76  and,  more  significantly,  oriented  nearly 
perpendicular  to  the  local  shoreline.  The  magnitude  of  the  annual 
resultant  wind  vector  is  a measure  of  the  directional  variation 
throughout  the  year,  which  differs  from  the  scalar  quantity  of  wind- 
speed.  The  mean  windspeed  was  not  significantly  different  the  2 
years.  However,  the  smaller  contribution  of  northerly  winds  in 
January  1975,  as  compared  to  January  1976,  caused  the  1974-75  result- 
ant vector  to  rotate  24°  clockwise  from  the  1975-76  vector.  Based 
on  this  wind  pattern,  the  total  wave  energy  flux  was  comparable  for 
the  2 years.  However,  the  net  longshore  energy  flux,  oriented  to 
the  south  both  years,  was  significantly  higher  in  1975-76  than  in 
1974-75. 
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WIND  OBSERVATIONS 


A.  2-15  JANUARY  1975 
8-18  JANUARY  1976 


Mini  FIRST  ANNUAL 
= SECOND  ANNUAL 


B.  14-26  MARCH  1975 
15-27  MARCH  1976 


SHORELINE  TREND  i 


C.  26  JULY  • i 
2 AUGUST  1974 

10-21  JUNE  1975 1 


D.  15-28  SEPTEMBER  1974 
22  SEPTEMBER  - 
5 OCTOBER  1975 


Figure  20.  Seasonal  distributions  of  wind  directions  observed  at 
North  Inlet. 


3.  Breaker  Parameters. 


a.  General.  Recently  developed  instrumentation  has  provided  a 
means  of  recording  directional  wave  spectra  on  a continuous  basis 
(Seymour,  et  al.,  1976).  More  commonly,  however,  wave  recorders  are  all 
nondirectional  (Harris,  1972)  and  estimates  of  directionally  dependent 
parameters  like  energy  flux  are  based  on  visual  observations  of  the  near- 
shore wave  field.  Littoral  processes  at  North  Inlet  were  measured  by 
observers  standing  in  the  breaker  zone.  Every  2 daylight  hours  during 
the  2-week  seasonal  observation  period,  the  observers  recorded  breaker 
heights  and  depths  by  holding  a graduated  rod  at  the  primary  breaker. 
Breaker  angle  was  measured  by  a hand-held  protractor  or  Brunton  com- 
pass, and  longshore  currents  by  tracing  the  movement  of  a neutrally 
buoyant  rubber  float  in  the  surf.  Periods  were  measured  by  counting 

11  successive  wave  crests  passing  a fixed  point  in  the  surf.  To  re- 
lease personnel  for  tidal  hydraulics  measurements  and  morphological 
mapping,  littoral  processes  in  1975-76  were  only  recorded  at  station 
DBI-25  which  was  least  influenced  by  shoaling  and  refraction  around 
the  ebb  tidal  delta.  Therefore,  the  parameter  statistics  presented 
in  this  section  should  reflect  conditions  at  a straight  segment  of 
beach  without  any  complicating  nearshore  bathymetry. 

b.  Breaker  Type.  Breakers  were  classified  as  spilling  and  plung- 
ing (Galvin,  1972) . Surging  and  collapsing  breakers  were  not  ob- 
served. The  beach  slope  at  station  DBI-25  was  gentle;  most  varia- 
tions in  the  primary  breaker  form  were  due  to  the  changing  deepwater 
steepness  of  the  waves.  Long-period  swell  arriving  from  offshore 
areas  of  generation  more  commonly  formed  plunging  breakers.  Locally 
generated  steeper  waves  commonly  break  by  spilling.  Summary  statis- 
tics (Tables  14  and  15)  for  2 years  of  observation  (July  1974  to  March 
1976)  indicate  that  the  breakers  were  evenly  distributed,  half  were 
spilling  and  half  were  plunging.  The  standard  error  of  the  mean  is 
quite  large,  so  the  percentage  differences  with  respect  to  wave  ap- 
proach direction  are  not  statistically  significant.  There  are  no 
statistically  significant  differences  between  the  seasonal  percent- 
ages either,  which  indicates  that  locally  generated  wind  waves  and 
distant  swell  account  for  about  the  same  relative  amount  of  the  ob- 
served breakers  at  all  seasons. 

c.  Breaker  Height.  The  hreaker  height  statistics  were  better 
than  those  for  breaker  type,  and  significant  seasonal  variations 
were  detected  (Table  16).  The  2-year  adjusted  average  was  59  centi- 
meters with  seasonal  variations  ranging  from  a low  of  54  centimeters 
in  March  to  a high  of  70  centimeters  in  September.  This  is  consis- 
tent with  the  wind  climate.  Weak,  variable,  and  commonly  offshore 
winds  in  spring  tend  to  generate  small  waves.  A high  frequency  of 
northeast  winds  in  fall  tends  to  generate  larger  waves  reaching 
shore  from  the  northeast.  These  results  vary  somewhat  with  those 

of  U.S.  Army,  Corps  of  Engineers,  Coastal  Engineering  Research 
Center  (1977)  which  indicated  that  the  south  Atlantic  coast  has  a 


Table  14.  Percent  spilling  waves  at  station  DBI-2S. 


Direction1 

NJ 

Mean 

Std.  dev. 

Median 

July  1974  - March  1976 

s 

141 

46.6 

33.0 

49.3 

N 

210 

54. 2 

29.4 

50.5 

u 

112 

40.2 

34.5 

30.1 

Average 

48.9 

32.2 

50 

Analysis  by  season 


July 

1974 

June 

1975 

Sept. 

1974  1975 

Jan. 

1975 

1976 

Mar. 

1975 

1976 

Mean 

3 

48.9 

56.4 

55.0 

44.5 

46.0 

36.0 

si. 5 

Std.  Dev. 

— 

32.6 

28.4 

31.8 

32.7 

32.4 

27.8 

35.2 

* N * Wave  approach  direction  between  east  and  north; 
S • Wave  approach  direction  between  east  and  south; 
//  - Wave  approach  direction  parallel  to  shoreline. 

2 Number  of  observations. 

3 No  reading. 


Table  15.  Percent  plunging  waves  at  station  DBI-25. 


Direction* 

N2 

Mean 

Std.  dev. 

Median 

July  1974  - March  1976  S 

141 

53.4 

33.0 

50.7 

N 

210 

45.5 

29.5 

40.5 

// 

112 

59.8 

34.5 

69.9 

Average 

51.1 

32.2 

50.0 

Analysis  by  season 


July 

1974 

June 

1975 

Sept, 

1974 

1975 

Jan. 

1975  1976 

Mar. 

1975 

1976 

Mean 

3 

51.0 

43.6 

45.0 

55.5 

54.0 

64.0 

47. S 

Std.  Dev. 

— 

32.7 

28.4 

31.8 

32.7 

32.4 

27.8 

35.2 

lN  ■ Wave  approach  direction  between  east  and  north; 
S • Wave  approach  direction  between  cast  and  south; 
//  ■ Wave  approach  direction  parallel  to  shoreline. 

2 Number  of  observations. 

'No  reading.. 


Table  16.  Breaker  heights  (in  meters)  at  station  DBI-25. 


Direction1 

N2 

Mean  i 

Std.  dev. 

Median 

July  1974  - March  1976 

S 

141 

0.17 

0.52 

N 

209 

m 

0.26 

0.60 

LL 

112 

0.53 

0.19 

0.50 

Average 

1 0.59 

0.23 

0.56 

Analysis  by  season 


July  June 

1974  1975 

Sept. 

1974  1975 

Jan. 

1975  1976 

Mar. 

1975  1976 

Mean 

3 0.57 

0.69  0.71 

0.59  0.55 

0.52  0.53 

*N  • Wave  approach  direction  between  east  and  north; 
S • Nave  approach  direction  between  east  and  south; 
//  • Nave  approach  direction  parallel  to  shoreline. 
*Nuaber  of  observations. 

*No  reading. 
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minimum  wave  height  of  about  43  centimeters  in  July  and  a maximum  of 
about  56  centimeters  in  September.  Since  these  latter  figures  are 
averages  from  a number  of  U.S.  Coast  Guard  stations  over  a large  region 
with  variable  deepwater  waves  and  refraction  and  shoaling  characteris- 
tics, this  level  of  difference  is  expected.  Breaker  depth,  angle  be- 
tween the  crest  of  the  breaker  and  the  shoreline,  and  wave  period  were 
measured,  but  not  tabulated  in  this  report.  However,  they  were  used  in 
following  calculations. 

4.  Longshore  Sediment  Transportation. 

The  best  established  predictive  equation  for  the  transportation  of 
sediment  on  beaches  is  based  on  the  longshore  component  of  wave  energy 
flux  (Komar,  1976;  Galvin  and  Vitale,  1977).  Since  wave  observations 
in  this  study  were  made  at  the  surf,  and  the  significant  wave  heights 
were  recorded,  the  correct  expression  is: 

Pls  = 2.78  x 10'2  x H^/2  x sin  2oib  (9) 

where  P^s  is  the  longshore  energy  flux  factor  measured  in  joules  per 
second  per  meter  of  beach  front  (U.S.  Army,  Corps  of  Engineers,  Coast- 
al Engineering  Research  Center,  1977) , Hb  is  the  breaker  height  (in 
meters),  and  ab  is  the  breaker  angle. 

Longshore  energy  flux  factors  at  station  DBI-25  for  the  eight  sea- 
sonal field  studies  are  summarized  in  Table  17.  Since  Pjs  does  not 
scale  linearly  with  Hb  or  ab,  no  mean  values  were  used  in  the  calcula- 
tions. Pjs  was  calculated  for  each  observation  to  include  all  observa- 
tion variability;  the  Pjs  values,  in  turn,  were  averaged  for  specified 
time  periods.  Mean  values  in  Table  17  refer  to  the  mean  of  all  readings 
giving  either  a northerly  or  southerly  flux  component.  The  resultant 
vector,  on  the  other  hand,  is  the  arithmetic  mean  of  all  calculated 
Pis  values,  with  (♦)  assigned  to  all  flux  vectors  oriented  to  the  right 
(south)  and  (-)  assigned  to  all  flux  vectors  oriented  to  the  left  (north). 

Using  only  field  studies  where  a precise  relationship  between  P^s 
and  the  volumetric  rate  of  longshore  sand  transport  had  been  established, 
the  U.S.  Army,  Corps  of  Engineers,  Coastal  Engineering  Research  Center 
(1977)  arrived  at  the  following  rating  equation: 

Qs  = 1.28  x 103  x PJs  (10) 

where  Qs  is  the  annual  longshore  transport  rate  (in  cubic  meters) . Ap- 
plying equation  (10)  to  the  data  in  Table  17,  sediment  transport  rates 
at  station  DBI-25  were  calculated  (Table  18).  Transport  to  the  north 
(Q  north)  and  south  (Q  South),  respectively,  was  computed  from  the  mean 
energy  flux  factors  to  the  north  and  south  for  each  2-week  period.  The 
net  and  gross  transport  rates  were  obtained  as  the  difference  and  sum, 
respectively,  of  Q north  and  Q south.  To  estimate  the  annual  transport 
rates,  the  calculated  values  for  each  2-week  period  were  assumed  to  be 
representative  for  an  entire  3-month  season.  The  validity  of  this  as- 
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Table  17.  Longshore  energy  flux  factors  (in  Joules  per  second 
per  meter)  at  station  DBI-25. 


Date 


July  1974  - Mar.  1976 


July  1974 
June  1975 


Sept.  1974 
1975 


Direction1 


Analysis  by  season 


N 


'Direction  along  shoreline. 

2Number  of  observations. 

3(+)  values  indicate  flux  toward  the  south. 
"*(-)  values  indicate  flux  toward  the  north. 


Resultant  vector 


+ 182 


Jan.  1975 

N 

22 

449 

S 

27 

289 

-42 

1976 

N 

26 

87 

S 

33 

355 

+ 160 

Mar.  1975 

N 

16 

S 

18 

+ 184 

1976 

N 

32 

S 

27 

144 

-32 

Table  18.  Computed  longshore  sediment  transport  rates  at  station  DBI-25. 


sumption  should  be  judged  both  in  light  of  the  breaker  height  data 
(Table  16) , which  indicate  a distinct  and  clearly  repetitive  seasonal 
pattern,  and  the  low  frequency  of  shoreline  crossing  extratropical  or 
tropical  cyclones. 

The  computed  total  rate  for  each  year  indicates  a consistent  gross 
transport  of  about  8 x 105  cubic  meters.  In  the  1974-75  project  year, 
the  net  transport  to  the  south  was  only  8.7  x 104  cubic  meters  or 
about  10  percent  of  the  gross.  This  result  is  consistent  with  the 
wind  climate  studies  which  indicated  that  the  resultant  wind  vector  fo 
1974-75  was  oriented  nearly  perpendicular  to  the  local  shoreline  with 
only  a small  southward  component  (Fig.  21).  However,  the  wind  vector 
for  1975-76  is  more  oblique  to  the  shoreline,  with  a much  larger  south 
ward  component,  consistent  with  a net  transport  to  the  south  of  about 
50  percent  of  the  gross  for  that  year,  or  about  3.9  x 105  cubic  meters 
Transport  rates  on  the  United  States  east  coast  are  about  2.0  x 105 
cubic  meters  per  year,  generally  to  the  south  (Wiegel , 1964).  The 
average  net  transport  rate  for  2 years  at  North  Inlet  is  2.4  x 105 
cubic  meters.  Studies  of  the  growth  of  Debidue  spit  and  the  North 
Inlet  ebb  tidal  delta  from  1925  to  1964  (Finley,  1976)  indicate  an  an- 
nual rate  of  increase  of  4.3  x 105cubic  meters  per  year. 

5 . Longshore  Currents. 

a.  Theory.  Longshore  current  velocities  measured  in  the  surf 
zone  at  station  DBI-25  are  the  vector  resultant  of  velocity  components 
due  to  the  oblique  approach  of  the  breaking  wave,  the  nearshore  cell 
circulation  system,  currents  generated  directly  in  the  nearshore  zone 
by  wind  stress  and,  perhaps,  components  of  a regional  circulation  sys- 
tem. It  is  generally  assumed  that  the  dominant  factor  in  the  genera- 
tion of  longshore  currents  is  the  oblique  approach  of  the  breaking 
wave.  The  two  most  important  parameters  determining  the  longshore  cur 
rent  velocity  are  the  wave  height  and  the  angle  between  the  wave  crest 
and  the  shoreline  (U.S.  Army,  Corps  of  Engineers,  Coastal  Engineering 
Research  Center,  1977).  A number  of  theories  explain  the  generation 
of  such  currents  and  predict  their  magnitude.  These  theories  are 
grouped  according  to  their  considerations  of  (a)  conservation  of  mass, 
(b)  conservation  of  energy  flux,  or  (c)  consideration  of  momentum 
flux  (Galvin,  1967).  Empirical  equations,  based  on  statistical  analy- 
sis of  parameters  describing  the  breaking  wave,  have  also  been  con- 
structed (Harrison,  1968;  Fox  and  Davis,  1972),  though  such  equations 
are  generally  of  limited  value  outside  the  specific  area  of  formula- 
tion. 

The  objective  of  this  section  was  to  evaluate,  through  the  use  of 
a stepwise  multiple- linear  regression  procedure,  whether  parameters 
descriptive  of  the  surf  zone  field  adequately  explain  the  observed 
variability  in  longshore  currents,  or  if  the  inclusion  of  additional 
environmental  parameters  could  significantly  improve  the  ability  to 
predict  such  current  velocity.  Linear  combinations  of  breaker  param- 


eters  and  three  proposed  equations  for  prediction  of  longshore  current 
velocity  were  used  in  the  evaluation.  Galvin  (1963),  basing  his  model 
on  the  continuity  of  water  mass,  arrived  at: 

V=kxgxmxTbx  sin  2ab  (11) 

where  V is  longshore  current  velocity,  g is  the  acceleration  of  gravity, 
m is  beach  slope,  ab  is  breaker  angle,  and  k is  a parameter  of  the 
breaker  form,  here  taken  as  1.0  (Galvin  and  Eagleson,  1965). 

The  equation  by  Longuet-Higgins  (1970)  was  derived  by  considerations 
of  the  conservation  of  momentum  flux  for  breaking  waves.  The  U.S.  Army, 
Corps  of  Engineers,  Coastal  Engineering  Research  Center  (1977)  empiri- 
cally determined  the  proportionality  constant  in  Longuet-Higgins'  equa- 
tion by  fitting  it  to  laboratory  data  by  Galvin  and  Eagleson  (1965)  and 
field  data  from  Putnam,  Munk,  and  Traylor  (1949).  Thus,  the  "modified" 
Longuet-Higgins  equation  reads: 

1/2 

V = 20.7  m(g  Hb)  ' sin  2ab  . (12) 

Fox  and  Davis'  (1972)  empirical  equation  for  their  Lake  Michigan 
data  set, 

Hb  (13) 

V = 2.98  sin  4ab 

was  also  included  in  the  testing.  This  equation  accounted  for  77.8  per- 
cent of  the  variance  in  the  current  velocity  data  at  Sheboygan,  Lake 
Michigan. 

b.  Data.  Descriptive  summary  statistics  for  all  longshore  current 
velocity  readings  obtained  between  July  1974  and  March  1976  are  presented 
in  Table  19.  The  mean  longshore  current  is  stronger  to  the  south  (35.8 
centimeters  per  second  versus  23.8  centimeters  to  the  north).  Extreme 
variability  in  current  velocities  is  demonstrated  by  the  fact  that  the 
standard  deviation  almost  equals  the  mean.  Inman  and  Quinn  (1952)  also 
found  in  their  study  of  longshore  current  variability  on  the  Pacific 
coast  that  the  standard  deviation  often  equaled  or  exceeded  the  mean. 
There  is  an  indication  in  the  data  that  fall  current  velocities  are 
slightly  higher  than  those  at  other  seasons. 

Appendix  H shows  that  the  1975-76  current  velocity  data  ranged  from 
73  centimeters  per  second  to  the  left  (north)  to  151  centimeters  per 
second  to  the  right  (south) , with  5 percent  of  currents  set  to  the  north 
stronger  than  51  centimeters  per  second  and  5 percent  of  those  set  to 
the  south  stronger  than  85  centimeters  per  second. 

The  1974-75  current  data  were  not  included  in  the  following  analy- 
sis because  of  the  large  number  of  missing  wind  observations.  The  1975- 
76  data  set  included  only  those  velocity  readings  which  occurred  at 
least  3 hours  after  sudden  reversals  in  wind  direction;  there  was  a lag 
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Table  19.  Longshore  current  velocities  (in  centimeters  per  second)  at  station  DBI-25. 

Date  Direction1 Mean  Std.  dev.  Median  Resultant  vectors 
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of  some  hours  between  reversals  in  the  longshore  wind  component  and 
corresponding  reversals  in  the  current.  A total  sample  of  250  data 
sets  remained. 

c.  Analysis.  Three  regression  procedures  were  utilized  in  this 
data  analysis : 

(a)  Simple  correlation  was  used  to  test  for  linear  rela- 
tionships between  any  pair  of  variables.  The  Pearson  correla- 
tion coefficient  is  a measure  of  the  degree  of  proportionality 
between  two  variables. 

(b)  Stepwise  regression  enters  one  independent  variable  at 
a time  until  all  are  entered  simultaneously.  Their  order  of  in- 
clusion is  determined  by  the  computer — the  independent  variable 
which  explains  the  largest  amount  of  variance  in  the  dependent 
variable  is  entered  first.  The  others  are  then  entered  in 
order  of  decreasing  variance.  The  proportion  of  the  total 
variance  in  the  dependent  variable  explained  by  an  independent 
variable,  or  a combination  of  variables,  is  expressed  by  the 
multiple  correlation  coefficient,  r2. 

(c)  Multiple  regression  enters  one  independent  variable  at 
a time  in  any  order  specified  by  the  investigator.  Thus,  the 
amount  of  variance  explained  by  any  independent  variable  of 
particular  interest  can  be  assessed. 

The  variable  names  used  in  the  computations  and  the  following 
analysis  are  defined  in  Table  20. 


Table  20.  Variable  names  used  in  multiple  regression  analysis 
of  littoral  processes. 

Name  Definition 

VEL  Observed  longshore  current  velocity  (in  centimeters  per  second) . 

(♦)  indicates  current  to  the  right  (south);  (-)  indicates  current 
to  the  left  (north). 

HGT  Breaker  height  (in  centimeters). 

PER  Breaker  period  (in  seconds). 

WIND  Wind  velocity  (miles  per  hour). 


WINDL  Longshore  component  of  wind  velocity  (miles  per  hour). 

ANGL  Angle  of  wave  orthogonal  relative  to  shoreline. 

WAVL  Sine  of  the  breaker  angle. 

VGAL  Velocity  calculated  by  Galvin's  (1963)  formula  (eq.  11). 

VELH  Velocity  calculated  by  Louguet-Higgins ' (1970)  formula  (eq.  12). 

VEFD  Velocity  calculated  by  Fox  and  Davis'  (1972)  formula  (eq.  13). 
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Table  21  is  a toatrix  of  the  Pearson  product-moment  correlation  coef- 
ficients for  these  10  variables  for  the  1975-76  annual  data  set  of  250 
observations. 


Table  21.  Pearson  correlation  coefficients  between  littoral 
process  variables  defined  in  Table  20. 


VEL 

HGT 

PER 

WINDL 

WAVL 

VGAL 

VELH 

VEL 

HGT 

0.37 

PER 

-0.11 

-0.06 

WINDL 

0.83 

0.36 

-0.07 

— 

WAVL 

0.68 

0.27 

-0.09 

0.68 

— 

VGAL 

0.64 

0.27 

-0.01 

0.66 

0.97 

— 

VELH 

0.69 

0.38 

-0.07 

0.70 

0.97 

0.95 

— 

VEFD 

0.70 

0.46 

-0.12 

0.69 

0.88 

0.81 

0.93 

The  measured  surf  zone  current  velocity  correlates  better  with  the 
longshore  component  of  the  wind  velocity  than  with  any  other  variable 
(r  = 0.83).  Wave  height,  period,  and  the  sine  of  the  breaker  angle  all 
show  a poor  correlation  with  the  current  velocity.  When  these  latter 
parameters  are  combined  as  suggested  in  equations  (10),  (11),  and  (12), 
the  predicted  and  observed  longshore  velocities  show  correlations 
ranging  from  r equals  0.64  for  Galvin's  (1963)  formula  to  r equals  0.70 
for  Fox  and  Davis'  (1972)  formula. 

Bivariate  data  plots  produced  by  the  computer  "SCATTERGRAM"  procedure 
(App.  H)  give  further  information  on  the  nature  of  these  correlations. 

The  plots  of  wave  height  and  period  versus  current  velocity  show  hardly 
any  discernible  trend.  Clearly,  these  parameters  do  not  exert  a dominant 
control  on  the  current  velocity.  Wind  velocity  correlates  with  the  long- 
shore current  velocity  at  r equals  0.32.  A large  amount  of  scatter  is 
observed  at  low  wind  velocities.  As  expected,  wave  angle  shows  a nega- 
tive correlation  ( r = 0.45)  with  the  velocity  because  of  the  direction 
convention.  The  variable  ANGL  (Table  20)  refers  to  the  orientation  of 
the  incoming  wave  orthogonal  relative  to  the  shoreline;  thus,  for  ANGL 
<90°,  the  breaking  wave  will  generate  a momentum  flux  and  a longshore 
current  to  the  right  (positive  direction) . For  ANGL  >90°  , the  wave- 
generated longshore  current  should  be  directed  to  the  left  (negative) . 

The  longshore  component  of  the  wind  velocity  correlates  with  the  cur- 
rent velocity  at  r equals  0.83.  The  scatter  is  rather  uniform  over  the 
entire  range  of  wind  velocities. 

If  the  longshore  component  of  the  wind  is  considered  as  the  only 
independent  variable,  a simple  regression  equation  for  the  longshore 
current  velocity  can  be  written  as: 

VEL  = 3.42  x WINDL  + 6.32  . (14) 
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This  is  a wholly  empirical  equation,  and  all  parameters  affecting 
the  longshore  current  velocity  are  lumped  into  WINDL.  Equation  (14) 
does  not  indicate  how  much  is  due  to  the  oblique  breaking  of  the  wind- 
generated waves.  However,  the  equation  is  of  predictive  value  for 
longshore  currents  off  Debidue  Island  (and  probably  elsewhere  under 
similar  environmental  conditions) . 

The  correlation  coefficients  between  observed  current  velocity  and 
that  predicted  by  the  three  equations  tested  are  all  reasonably  high. 
This  indicates  proportionality.  However,  the  magnitude  of  the  pre- 
dicted velocity  generally  differs  substantially  from  that  observed 
as  evidenced  by  the  proportionality  factors  in  the  following  three 
regression  equations  (all  velocities  are  in  centimeters  per  second). 


VEL  = 

0.33 

VGAL  ♦ 

5.18, 

r = 

0.57 

(IS) 

VEL  = 

0.45 

VELH  ♦ 

4.6, 

r = 

0.68 

(16) 

VEL  = 

2.0 

VEFD  + 

4.1, 

r = 

0.70 

(17) 

Galvin's  (1963)  and  Longuet-Higgins ' (1970)  equations  are  both 
semiempirical;  their  derivations  are  based  on  the  assumption  that 
the  longshore  current  is  caused  by  the  breaking  of  the  obliquely  in- 
cident wave.  Both  equations  predict  current  velocities  higher  than 
those  observed  at  Debidue  Island.  The  empirical  equation  of  Fox  and 
Davis  (1972)  predicts  current  velocities  which  are  generally  too  low 
by  a factor  of 

To  test  for  the  relative  importance  of  the  independent  variables 
in  explaining  the  variance  in  the  dependent  one,  a stepwise  regression 
analysis  was  used.  The  test  was  performed  both  for  the  pooled  annual 
process  data  and  for  each  seasonal  data  set.  Results  are  summarized 
in  Table  22. 


Table  22.  Percent  of  the  variance  in  VEL. 1 


Variable 

Annual 

1975 

1976 

June 

Sept. 

Jan . 

Mar. 

WINDL 

70 

81 

63 

57 

69 

VGAL 

1 

0 

0 

2 

2 

VELH 

2 

0 

1 

2 

4 

VEFD 

4 

0 

7 

1 

6 

Explained  by  each  of  four  independent  variables  entered 
successively  in  a stepwise  regression  analysis. 


For  all  data  sets,  the  longshore  component  of  the  wind  velocity, 
WINDL,  proved  to  be  the  independent  variable  explaining  most  of  the 
observed  variance  in  current  velocity,  VEL.  The  multiple  correlation 
coefficient  (r2) , for  VEL  versus  WINDL  ranged  from  0.81  for  the  June 
1975  data  set  to  0.57  for  the  January  1976  data.  Expectedly,  data 
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noise  was  at  a maximum  in  January  and  minimum  in  June  because  of  the 
different  weather  conditions  under  which  the  field  observers  had  to 
operate.  It  is  quite  significant  to  note  that  for  the  data  set  which 
was  expected  to  be  the  most  reliable  (June  1975) , WINDL  alone  ex- 
plained 81  percent  of  the  variance.  For  this  data  set,  the  inclusion 
of  breaker  parameters  combined  into  the  predictive  equations  of  Galvin 
(1963) , Longuet-Higgins  (1970) , and  Fox  and  Davis  (1972) , does  not  im- 
prove the  multiple  correlation  coefficient.  For  the  other  data  sets, 
these  variables  add  a few  percentages  of  explained  variance. 

Based  on  these  considerations,  the  best  predictive  equation  for 
longshore  current  velocities  at  Debidue  Island  appears  to  be  that  de- 
rived from  the  data  set  of  62  observations  from  June  1975.  This  equa- 
tion is: 

VEL  = 4.36  x WINDL  ♦ 12.9  . (18) 

In  both  equations  (14)  and  (18),  VEL  is  measured  in  centimeters  per  sec- 
ond and  WINDL  in  miles  per  hour.  For  a longshore  wind  velocity  of  10 
miles  per  hour,  equation  (18)  predicts  a current  velocity  of  56.5  centi- 
meters per  second;  whereas,  equation  (14)  predicts  40.6  centimeters  per 
second,  a difference  of  39  percent.  For  wind  velocities  above  30  miles 
per  hour,  the  difference  between  the  two  predicted  current  velocities 
becomes  less  than  30  percent. 

Using  a standard  multiple  regression  to  test  for  the  amount  of 
variance  explained  by  VGAL,  VELH,  and  VEFD,  if  they  are  entered  as  the 
first  independent  variable  in  three  separate  calculations,  VGAL  can 
explain  32  percent  of  the  variance,  VELH,  46  percent,  and  VEFD,  49 
percent  for  the  annual  data  set. 


d.  Discussion.  A series  of  regression  equations  can  be  developed 
depending  on  which  variables  or  combinations  of  variables  are  included. 
Table  23  summarises  the  two  types  of  equations  derived  in  this  study; 
the  first  (eq.  19)  is  the  regression  equation  derived  by  treating  each 
measured  environmental  parameter  as  an  independent  variable;  the  three 
subsequent  equations  (eqs.  20,  21,  and  22)  are  derived  by  combining  the 
pertinent  breaker  parameters  as  suggested  in  Galvin's  (1963),  Longuet- 
Higgins'  (1970),  and  Fox  and  Davis'  (1972)  predictive  formulas.  Breaker 
parameters  not  accounted  for  in  the  predictive  formulas  are  included  as 
independent  variables.  The  amount  of  variance  in  VEL  explained  by  each 
parameter  combination  is  indicated. 

Based  on  a similar  analysis,  Harrison  (1968)  derived  a regression 
equation  for  longshore  current  velocity  at  Virginia  Beach,  Virginia. 
However,  simultaneous  wind  measurements  were  not  obtained  in  his  study, 
and  the  question  of  wind  stress  or  breaking  wave  dominance  in  current 
generation  was  left  unexplained.  Based  on  earlier  studies  by  Sonu, 
McCIoy,  and  McArthur  (1967),  at  Nags  Head,  North  Carolina,  Harrison 
discounted  the  importance  of  wind  stress.  The  results  at  Debidue  Island 
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Table  23.  Regression  equations  for  longshore  current  velocity.1 
Eg.  No.  Equation  r2 

19  VEL  = 2.78  WINDL  + 44.2  WAVL  + 0.12HGT  - 0.85  PER  ♦ 5.94  0.722 

20  VEL  = 2.73  WINDL  ♦ 0.69  VEFD  + 4.72  0.726 

21  VEL  = 2.83  WINDL  + 0.13  VELH  - 0.94  PER  «•  13.28  0.719 

22  VEL  = 2.91  WINDL  ♦ 0.06  VGAL  » 0.12HGT  - 1.15  0.714 

1 A1 1 measured  environmental  parameters  were  included  in  the  analysis. 
r2  is  the  total  proportion  of  variance  explained. 

clearly  demonstrate  that  wind  stress  can  be  a most  significant  factor 
in  nearshore  current  generation.  The  statistical  approach,  unfortunately, 
does  not  explain  the  actual  mechanism  of  energy  transfer  from  the  wind  to 
the  current  system.  All  current  measurements  were  obtained  by  surface 
floats.  Further  research  is  needed  to  test  the  three-dimensional  struc- 
ture of  this  flow  field  and  to  evaluate  various  genetic  mechanisms. 

V.  SEASONAL  MORPHOLOGIC  CHANGE 


1 . Introduction. 

Mapping  of  the  intertidal  morphology  at  North  Inlet  during  the  quar- 
terly field  studies  (June  1975  to  May  1976),  low  aerial  overflights,  and 
profiling  of  the  channels  and  inlet-associated  beaches  on  a monthly 
basis  provided  the  data  base  for  the  evaluation  of  the  seasonal  changes 
in  morphologic  characteristics.  Sedimentary  subenvironments  were  estab- 
lished according  to  the  strand-plain  classification  of  Fisher  and  Brown 
(1972)  and  the  tidal  delta  models  of  Hayes,  et  al.  (1973).  To  facili- 
tate a systematic  evaluation  of  morphological  development  over  the  year, 
a detailed  description  of  the  observed  change  in  each  subenvironment  was 
made. 

Tidal  deltas  and  associated  characteristic  beach  features  are  formed 
as  a result  of  the  interruption  to  the  continuous  longshore  sediment 
transport  system  by  strong  tidal  currents  flowing  essentially  perpen- 
dicular to  the  regional  shoreline.  Although  additional  research  is 
needed  on  the  origin  and  initial  development  of  the  tidal  deltas  once 
an  inlet  channel  has  been  cut,  the  characteristics  of  morphology  and 
sediment  transport  patterns  on  established  tidal  deltas  are  well  docu- 
mented. Hayes,  et  al.  (1973)  and  Oertel  (1972)  documented  how  the 
typical  ebb  tidal  delta  morphology  of  mesotidal  inlets  results  from  the 
balance  between  a seaward  transportation  of  sediment  by  the  ebb  flow  and 
a landward  transport  due  to  wave  swash  and  flood  currents.  Hayes  (1975) 
considered  the  changes  in  inlet  morphology  with  tidal  range  and  developed 
models  for  microtidal,  mesotidal,  and  macrotidal  estuaries.  Detailed 
studies  of  tidal  delta  bed- form  distributions  and  migration  patterns 
have  been  done  at  a number  of  mesotidal  inlets  (Oertel,  1972;  Booth- 
royd  and  Hubbard,  1975;  Hine,  1975;  Hubbard,  1975;  Finley,  1975,  1976). 
For  some  inlets,  detailed  sediment  circulation  patterns  have  been  pro- 
posed (Masterson,  Machemehl , and  Cavaroc,  1973;  Hine,  1975;  FitzGerald, 
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Nummedal,  and  Kana,  1977).  Past  studies  of  inlet  morphologic  changes 
have  been  based  on  the  assumption  that  MSL  throughout  the  year  is  con- 
stant, which,  as  pointed  out  in  Section  III,  8,  is  not  the  case. 

Future  inlet  studies,  involving  interpretations  of  morphologic  change, 
should  consider  the  fact  that  sea  level  on  the  United  States  east 
coast  varies  some  40  centimeters  during  the  year.  The  first  objective 
of  this  analysis  of  North  Inlet  morphology  was  to  document  how  this 
sea  level  fluctuation  has  affected  particularly  sensitive  inlet  sub- 
environments. 

As  pointed  out  by  Escoffier  (1940)  and  reiterated  by  Inman  and 
Frautschy  (1966)  and  others,  a stable  inlet  appears  to  maintain  a 
balance  between  the  scouring  actions  of  the  tidal  currents  and  the  in- 
filling by  sediment  supplied  through  the  longshore  drift  system.  If 
this  concept  is  correct.it  should  be  documented  that  the  time  his- 
tories of  channel  erosion  and  erosion  on  the  inlet  adjacent  beaches 
are  out  of  phase.  The  extensive  morphological  data  collected  at  North 
Inlet  over  a year,  involving  a wide  range  of  weather  conditions,  make 
such  a test  possible.  This  constituted  the  second  objective  of  the 
morphologic  analysis. 

2.  General  Morphology  and  Sediment  Distribution. 

A depositional  environment  map  (Fig.  22)  was  prepared  from  verti- 
cal aerial  photos,  data  on  historical  changes  in  inlet  morphology 
and  field  mapping.  The  terminology  follows  the  strand-plain  classi- 
fication of  Fisher  and  Brown  (1972). 

Abandoned  flood  tidal  deltas  and  associated  washover  units  were 
identified  on  the  landward  side  of  Debidue  Creek.  The  inlet  has  mi- 
grated to  the  south,  and  Debidue  Spit  has  been  extended  by  about  2 
kilometers  since  1878  (Finley,  1976).  The  abandoned  flood  tidal 
deltas  reflect  short-term  stillstands  in  the  southward  migration  of 
the  inlet.  The  present  inlet  is  backed  by  a flood  tidal  delta  and 
washover  berm  at  the  mouth  of  Sixty  Bass  Creek  and  flanked  by  supra- 
tidal  barrier  beaches  with  small  washover  fans  and  intertidal  and 
shallow  subtidal  bars  and  platforms  forming  a seaward  continuum  of 
North  and  Debidue  Islands;  i.e.,  the  ebb  tidal  delta. 

Sediment  size- and  distribution  at  North  Inlet  are  related  to  a 
number  of  factors  that  vary  on  a regional  scale.  The  inlet  is  lo- 
cated at  the  southern  end  of  an  arcuate  shoreline  that  extends 
approximately  100  kilometers  from  the  North  Carolina  border  to  Winyah 
Bay.  For  most  of  these  100  kilometers,  the  present  beach  borders 
Pleistocene  upland  except  in  the  immediate  vicinity  of  North  Inlet  and 
Murrells  Inlet,  35  kilometers  farther  south.  The  Holocene-Pleistocene 
boundary  at  North  Inlet  is  about  5 kilometers  inland.  Littoral  drift 
along  the  coast  is  derived  from  mid-Pleistocene  beach  deposits.  An 
average  size  of  2.48  phi  (.0.175  millimeter)  was  determined  for  these 
sediments  (Brown,  1975).  The  arcuate  strand  is  broken  by  few  tidal 
inlets.  No  major  rivers  exist  within  this  coastal  segment.  The 
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Figure  22.  Major  depositional  environments  at  North  Inlet. 
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nearest  river  sources  are  located  south  of  North  Inlet.  Winyah  Bay, 
an  estuary  for  the  Sampit,  Pee  Dee,  Black,  and  Waccamaw  Rivers,  is 
connected  to  the  sea  approximately  15  kilometers  to  the  south.  About 
15  kilometers  farther  south  is  the  Santee  River  delta.  Sand-size 
material  reaching  North  Inlet  from  these  sources  is  minimal;  however, 
suspended  fine-grained  organic  and  inorganic  sediments  may  reach  the 
inlet  (Finley,  1976).  The  presence  of  coarser  sediment  south  of  North 
Inlet  (average  size  2.01  phi  or  0.248  millimeter)  indicates  a higher 
percentage  of  fluvially  derived  sand  (Brown,  1975). 

The  sediment  sampling  program  at  North  Inlet  was  designed  to 
evaluate  seasonal  variations  in  sediment  size.  Locations  of  sedi- 
ment sampling  stations  1 to  21  are  shown  in  Figure  23,  with  respect 
to  the  channel  bathymetric  profiles  in  Figure  8.  Most  sample  sta- 
tions were  the  same  as  those  used  by  Finley  (1976) . The  beach  sedi- 
ments were  obtained  midway  between  MSL  and  mean  high  water  after 
removing  the  surface  layer  of  coarse  shell  fragments  and  taking  a 
2-centimeter-deep  sample  immediately  below.  Channel  bed  samples  were 
obtained  with  a Van  Veen  grab  sampler.  Standard  sieving  techniques 
were  applied  in  the  laboratory.  Size  parameters  were  calculated, 
using  Folk's  (1968)  statistics.  The  percentage  of  shell  hash  coarser 
than  2 millimeters  in  each  sample  was  recorded.  The  sample  means 
and  standard  deviations  (see  Table  24)  refer  to  the  sand-sized  frac- 
tion (particles  less  than  2 millimeters). 

No  significant  seasonal  variation  in  grain  size  was  found  at 
any  environment  sampled.  However,  there  is  a distinct  difference 
between  beach  and  channel  sands.  Beach  sands  have  a total  sample 
mean  of  2.375  phi  (0.19  millimeter);  whereas,  channel  sands  have  a 
mean  of  1.214  phi  (0.43  millimeter),  a difference  of  more  than  1 phi 
unit.  More  than  90  percent  of  the  beach  samples  were  well  sorted 
or  very  well  sorted  (S  = 0.50  phi)  in  contrast  to  moderate  to  poor 
sorting  (S  = 0.71  to  2 phi)  for  more  than  85  percent  of  the  channel 
samples.  Shell  fragments  or  coarse  shell  hash  was  very  diagnostic 
of  the  channel  sands,  ranging  from  0.1  to  80  percent  of  the  total 
sample  weight.  This  type  of  variation  is  illustrated  by  the  samples 
collected  at  station  13  (Table  24),  where  a large  amount  of  shell 
hash  is  deposited  in  the  thalweg  of  the  main  ebb  channel.  If  this 
area  was  sampled  continuously  over  a short  period  of  time,  a large 
variation  in  the  amount  and  size  of  shell  fragments  would  be  found 
because  of  changes  in  current  strength  over  the  tidal  cycle.  Sedi- 
ment characteristics  a;  stations  17  and  19  (Fig.  23)  are  not  properly 
defined  by  the  statistics  in  Table  24  as  these  locations  are  at  out- 
crops of  firm,  cohesive  lagoonal  muds.  At  these  thalweg  locations, 
gray  mud  was  often  brought  up  with  the  anchor. 

3.  Morphologic  Response  Characteristics. 

a.  Introduction.  Four  individual  subenvironments  at  North  Inlet 
were  subjected  to  detailed  morphologic  mapping  on  a seasonal  basis: 
inlet  channels,  flood  and  ebb  tidal  deltas,  and  barrier  island  beaches. 
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Table  24.  Quarterly  sediment  data  (expressed  in  phi  units) 


NOTE:--See  Figure  23  for  location  of  stations. 
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To  facilitate  quantitative  comparisons  of  morphology,  a systematic,  con- 
sistent, and  accurate  field  surveying  program  was  emphasized  throughout 
the  study. 


b.  Inlet  Channels.  Temporal  fluctuations  in  cross-sectional  areas 
for  Town  and  Jones  Creeks  (summarized  in  Table  4)  are  graphically  shown 
in  Figure  24  to  facilitate  the  discussion  of  annual  changes.  The  cross- 
sectional  change  (in  square  meters)  is  plotted  relative  to  the  area  of 

the  June  1975  survey  (chosen  as  0).  The  mean  cross-sectional  area  of  * 

Jones  Creek  is  about  half  that  of  Town  Creek.  Consequently,  the  vertical 
scales  in  Figure  24  are  adjusted  such  that  a given  amount  of  deviation 
represents  the  same  relative  change  in  the  two  channel  flow  areas. 

Both  channels  underwent  an  18-percent  variation  in  flow  area  through- 
out the  year,  compared  to  15  percent  for  the  inlet  throat  section.  The 
changes  were  roughly  in-phase  (Fig.  24).  Channel  cross  sections  were 
relatively  stable  through  the  summer,  followed  by  a scouring  trend  in 
early  fall.  Deposition  occurred  in  November  and,  apparently,  continued 
through  mid-winter  to  reach  an  alltime  high  in  January.  For  the  remainder 
of  the  winter  and  the  spring,  two  cycles  of  erosion  and  deposition  oc- 
curred, one  from  January  to  March,  the  other  from  March  to  May.  Both 
creeks  appear  to  respond  to  similar  external  events,  and  both  have  a 
clear  seasonal  pattern  of  channel  fill  and  scour.  An  exact  analysis  of 
flow  area  changes  of  the  inlet  throat  is  impossible  due  to  insufficient 
data,  but  inspection  of  the  data  indicates  a seasonal  behavior  similar 
to  that  of  the  main  creeks. 

c.  Flood  Tidal  Delta.  Both  the  abandoned  flood  tidal  deltas  and  the 
active  one  are  similar  in  that  each  formed  at  the  junction  of  tidal  creeks, 
and  washover  deposits  exist  on  their  landward  sides.  The  active  flood 
tidal  delta  is  formed  at  the  junction  of  Jones  and  Town  Creeks.  The  tri- 
angular planform  geometry  (Fig.  25)  differs  somewhat  from  the  typically 
lobate  flood  tidal  delta  geometry  described  for  New  England  (Hayes,  1975; 

Boothroyd  and  Hubbard,  1975;  Hine,  1975)  and  the  Outer  Banks  of  North 
Carolina  (Nummedal , et  al.,  1977).  The  typical  flood  tidal  delta  is 
generally  located  midchannel  and  is  open  to  ebb  process  modification 
along  its  flank.  As  such,  the  complete  symmetrical  flood  delta  can  only 
form  in  a large  embayment  or  a wide  channel  (Barwis  and  Hubbard,  1976) . 

Of  the  five  components  which  characterize  the  typical  flood  tidal  delta 
(i.e.,  flood  channel,  flood  ramp,  ebb  spit,  spillover  lobe,  and  ebb 
shield),  the  ebb  shield  is  the  least  developed  at  North  Inlet  (Fig.  25). 

This  is  due  to  the  welding  of  the  landward  margin  of  the  delta  to  the 
marsh. 

Between  July  1972  and  March  1973,  a 42-meter  westward  migration  of 
the  landward  edge  of  the  delta  was  measured  (Finley,  1976).  Continued 
monitoring  the  following  2 years  indicated  little  farther  migration,  and 
only  minor  changes  along  the  northern  delta  margin. 

During  the  1975-76  field  season,  the  flood  tidal  delta  changes  were 
monitored  by  monthly  fathometer  transects  (Fig.  8,  profile  D-D'),  quar- 
terly oblique  aerial  photography,  and  a theodolite  and  rod  survey  (Fig.  26). 
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Diagram  of  a flood  tidal  delta  (after  Hayes,  et  al 
1973),  and  oblique  air  photo  of  the  North  Inlet 
flood  tidal  delta,  June  1975,  showing  (1)  ebb  spit 
(2)  flood  ramp,  (3)  ebb  shield,  (4)  flood  channel, 
and  (5)  spillover  lobe. 
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Figure  26.  Topographic  map  of  the  North  Inlet  flood  tidal  delta 


The  bathymetry  recorded  by  profile  track  D-D'  shows  the  wide. 


shallow  delta  platform,  or  flood  ramp,  between  the  two  creeks.  Evi- 
dence of  westward  migration  into  the  aalt -marsh  environment  and  into 
Sixty  Bass  Creek  is  the  washover  berm  deposition  and  the  burial  of 
oyster  reefs  (Figs.  25  and  26).  A steep,  eroding  beach  face  con- 
taining tree  stumps  and  Indian  pottery  fragments  is  in  front  of  the 
southern  washover  deposits.  Dates  of  the  pottery  range  from  500  to 
1500  B.P.,  indicating  seasonal  occupation  by  numerous  Indian  tribes 
(R.J.  Widmer,  Department  of  Archeology,  University  of  South  Carolina, 
personal  communication,  1976). 

Little  seasonal  variation  in  morphology  is  detected  at  the  flood 
tidal  delta.  Its  stability  is  a significant  characteristic  of  the 
North  Inlet  system. 

d.  Ebb  Tidal  Delta.  The  present  ebh  tidal  delta  began  its  growth 
in  the  late  1930's  when  North  Inlet  stabilized  near  its  present  posi- 
tion (Finley,  1976).  By  1964,  additional  spit  growth  at  each  island 
terminus  had  confined  the  inlet  tidal  prism  to  a smaller  cross-sec- 
tional area,  and  the  ebb  delta  had  further  developed.  The  five  major 
ebb  tidal  delta  components  (Hayes,  et  al.,  1973)  are  well  defined  by 
the  present  morphology  (Fig.  27).  The  main  ebb  channel  trends  to  the 
southeast  and  extends  from  the  inlet  throat  to  the  terminal  lobe 
(Fig.  8).  Channel-margin  linear  bars  flank  the  main  ebb  channel. 

The  linear  bar  on  the  north  side  of  the  channel  is  gently  arcuate  and 
located  landward  of  its  counterpart  on  the  south  side  (Fig.  8)  which 
demonstrates  the  downdrift  offset  nature  of  the  inlet  (Hayes,  Gold- 
smith, and  Hobbs,  1970).  Marginal  flood  channels  at  North  Inlet 
differ  greatly  on  either  side  of  the  inlet.  The  northern  channel  is 
generally  indistinguishable  from  the  swash  platform  itself  and  is 
commonly  characterized  by  large-scale,  flood-oriented  bed  forms  along 
the  recurved  spit  margin.  The  southern  channel  is  deeper  and  narrower, 
separating  the  southern  channel-margin  linear  bar  from  North  Island  and 
connecting  the  main  ebb  channel  to  the  swash  platform  interior.  Swash 
bare  are  located  near  the  distal  part  of  the  main  ebb  channel  and  along 
the  seaward  boundary  of  the  swash  platform.  Well- developed  bars  on  the 
south  side  of  the  delta  have  migrated  landward  and  become  attached  to 
North  Island  over  the  past  4 years. 

During  the  1975-76  field  season,  the  ebb  delta  morphology  was 
recorded  on  a quarterly  basis.  Twelve  beach  and  nearshore  profiles, 
extending  generally  out  to  the  7-meter  contour,  or  about  3 kilometers 
offshore,  were  regularly  monitored  (Fig.  28).  Computer  routines 
were  used  to  remove  horizontal  distortions  in  the  offshore  profiles 
and  to  plot  all  corresponding  profiles  at  the  same  scale  (see  App.  I 
for  profiles) . 

In  addition  to  profiling,  six  detailed  theodolite  and  rod  surveys 
of  the  channel -margin  linear  bars  and  swash  bars  were  made.  Seasonal 
evaluations  of  bed  form  and  bar  morphology  were  also  based  on  oblique 
and  vertical  aerial  photos.  Weather  conditions,  equipment  failures, 


and  logistical  problems  prevented  a complete  seasonal  analysis  of  the 


Diagram  of  an  ebb  tidal  delta  (after  Hayes,  et  al., 
1973) , and  oblique  air  photo  of  the  North  Inlet 
tidal  delta,  January  1976,  showing  (1)  main  ebb 
channel,  (2)  channel-margin  linear  bars,  (3)  marginal 
flood  channels,  (4)  swash  bars,  and  (5)  terminal  lobe 


Figure  27 


entire  ebb  tidal  delta.  Sea  conditions  were  often  not  favorable  for 
the  monitoring  of  offshore  profiles,  especially  those  located  over  the 
shoals.  A detailed  bathymetric  map  of  the  North  Inlet  ebb  tidal  delta 
was  constructed  from  the  January  1976  profile  data  (Fig-  29).  The 
evaluation  of  the  last  12  years  of  ebb  delta  development  is  based  on 
the  comparison  between  this  1976  map  and  that  constructed  from  the 
1964  U.S.  Coast  and  Geodetic  Survey,  Chart  No.  8838  (Fig.  29).  The 
general  bathymetric  development  of  the  ebb  delta  since  1964  provides 
the  background  for  the  short-term  morphological  change  detected  in 
the  1975-76  field  season. 

The  following  four  major  bathymetric  changes  have  occurred  since  1964 

(a)  The  7-meter  contour  northeast  of  the  inlet  has  moved 
landward  along  with  the  5-  and  6-meter  contours.  This  repre- 
sents steepening  of  the  northeast  delta  margin  relative  to  the 
1964  configuration. 

(b)  In  the  southeast  map  area,  a topographic  high  out- 
lined by  the  4-meter  contour  in  1976  occupies  an  area  which 
was  gently  sloping  seaward  from  5 to  6 meters  in  1964. 

(c)  At  least  three  distinctive  shoals  (less  than  2 meters 
deep)  have  developed  adjacent  to  the  main  ebb  channel,  one  at 
the  terminal  lobe  and  one  at  either  side  of  the  main  channel 
on  the  swash  platform. 

(d)  The  main  ebb  channel  has  deepened  and  is  now  oriented 
straight  toward  the  southeast. 

Inadequate  long-term  process  data  make  an  attempted  interpretation 
of  these  changes  purely  conjectural.  However,  the  overall  skewness  of 
the  ebb  tidal  delta  reflects  a southward  net  littoral  transport. 

In  the  following  paragraphs,  three  components  of  the  ebb  tidal  delta 
are  discussed  in  detail  with  respect  to  their  short-term  morphological 
changes.  This  discussion  intends  both  to  document  the  rates  and  nature 
of  ebb  delta  changes  as  base-line  data  for  typical  east  coast  mesotidal 
inlets  and  to  provide  an  understanding  for  the  effects  of  annual  changes 
in  sea  level. 

The  northern  channel -mar gin  linear  bar  has  undergone  growth  toward 
the  southwest,  changes  in  surface  elevation,  and  a widening  of  the  spill- 
over channel  on  its  seaward  side  since  observations  began.  Between  July 
and  September  1974,  a 70-centimeter  increase  in  bar  margin  elevation 
occurred,  but  by  March  1975,  there  was  a 70-centimeter  decrease  in  ele- 
vation (Finley,  1976).  This  erosion  was  probably  related  to  a northeast 
storm  in  September  1974.  About  a year  later,  between  January  and  March 
1976,  a period  of  accretion  occurred.  Figure  30  is  a photo  and  contour 
map  of  the  bar  morphology  in  January  1976,  showing  the  broad  and  shallow 
slope  toward  the  swash  platform,  the  steep  channel  margin,  and  ridgelike 
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topographic  highs.  Although  the  swash  platform  is  not  clearly  repre- 
sented by  either  the  contour  map  or  the  photo,  the  open  contours  north 
of  the  bar  imply  a broad  and  shallow  topography.  The  120-centimeter 
contour  defines  well  the  overall  arcuate  geometry.  The  steepness  of  the 
channel  margin  (i.e.,  southern  boundary)  is  indicated  by  the  closeness 
of  the  contours  and  is  most  pronounced  near  the  three  topographic  highs. 
The  "saddle"  at  A (Fig.  30)  is  located  just  across  the  ebb  channel  from 
the  southern  marginal  flood  channel.  This  means  that  the  northwest 
section  of  the  bar  is  landward  of  the  southern  marginal  flood  channel, 
giving  the  inlet  a slight  degree  of  downdrift  offset. 

By  comparing  a topographic  map  of  the  bar  in  March  1976  (Fig.  31) 
to  that  in  January  (Fig.  30),  changes  are  quite  noticeable.  The  bar 
has  accreted  to  the  northeast  while  maintaining  the  same  arcuate  geo- 
metry and  steep  channel  margin.  The  -70-centimeter  contour  is  con- 
tinuous around  the  entire  bar,  and  the  two  topographic  highs  are  sep- 
arated only  at  the  -45-centimeter  contour  level.  During  a 2-month 
period,  in  the  late  winter  of  1976,  the  northern  channel-margin  linear 
bar  underwent  significant  lateral  and  vertical  accretion. 

The  southern  channel-margin  linear  bar  has  grown  almost  in  its  en- 
tirety since  the  obliteration  of  its  predecessor  by  the  February  1973 
northeast  storm.  Following  that  storm,  rapid  swash  bar  accretion  and 
coalescence  formed  a long  linear  bar  between  May  1973  and  July  1974 
(Finley,  1976).  By  January  1975,  the  bar  had  increased  in  width.  One 
year  later  (March  1976,  Fig.  32),  the  bar  was  about  100  meters  wide 
and  had  a length  comparable  to  that  of  its  northern  counterpart.  Un- 
like the  northern  bar,  the  southern  bar  is  topographically  more  iso- 
lated on  the  swash  platform.  This  is  primarily  due  to  major  differences 
in  marginal  flood  channel  morphology  and  topography  on  the  interior  of 
the  swash  platform.  A relatively  deep  and  narrow  marginal  flood  chan- 
nel completely  separates  the  southern  bar  from  North  Island.  In  ad- 
dition, the  low  elevation  of  the  adjacent  swash  platform  enhances  bar 
relief.  The  morphology  of  the  southern  bar  differs  from  that  of  the 
northern.  It  has  a gently  sloping  perimeter  even  at  the  ebb  channel  side, 
a surface  dominated  by  plane  beds  or  ripples,  and  forms  one  continuous 
topographic  high.  Because  of  the  distinctive  light  shade  of  the  dry  sands 
of  the  bar  crest  at  low  tide,  aerial  photography  proved  to  be  an  excellent 
tool  to  document  the  seasonal  changes  in  bar  configuration. 

Oblique  air  photos  were  obtained  in  June  and  October  1975  and  in 
January  and  March  1976  (Fig.  33).  In  June  1975,  the  landward  half 
of  the  bar  had  a distinct  crest  along  the  channel  side;  on  the  sea- 
ward half,  the  CTest  formed  at  the  swash  platform  bar  margin.  By 
October,  only  the  landward  half  of  the  bar  had  a distinct  crest.  In 
January  1976,  a continuous  crest  was  present  along  the  swash  platform 
bar  margin.  In  March,  the  crest  had  again  moved  over  to  the  channel 
side  of  the  bar  at  its  landward  part.  Although  the  March  morphology 
is  not  nearly  as  distinct  as  that  of  the  previous  June,  it  has  the 
same  general  distribution  of  features.  A distinct  "cyclic"  change 
in  bar  crest  location  characterized  the  annual  changes  of  the  mor- 
phology of  the  southern  channel-margin  linear  bar. 


88 


Figure  33..  Seasonal  changes  in  configuration  of  the  southern  channel- 
margin  linear  bar.  Light  areas  are  topographic  highs. 
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The  dominant  processes  controlling  the  morphologies  of  the  north- 
ern and  southern  channel-margin  linear  bars  were  identified.  Strong 
tidal  currents  and  relatively  weak  wave  action  are  responsible  for 
the  overall  morphology  and  the  large  bed  forms  of  the  northern  chan- 
nel-margin linear  bar.  The  southern,  seaward  offset  channel -margin 
bar,  on  the  contrary,  is  dominated  by  wave  swash  causing  a smooth 
compacted  surface  with  few  bed  forms,  similar  to  a beach  face. 

The  North  Island  swash  bar  complex  has  been  monitored  since  1974 
when  it  was  observed  that  the  entire  bar  was  moving  landward  and 
that  its  point  of  beach-face  attachment  was  moving  to  the  north. 

The  large  landward-pointed  swash  bar  had  been  divided,  and  both  parts 
were  migrating  landward  by  September  1974  (Fig.  34, A).  In  March  1975, 
changes  in  the  southern  part  were  apparent  (Fig.  34, B).  Vertical  air 
photos  taken  in  April  1974  and  March  1975  were  used  to  determine  that 
the  southern  part  had  moved  417  meters  or  an  average  of  36.3  meters 
per  month  (Finley,  1976).  The  northernmost  pointed  feature  moved 
much  slower.  During  the  1975-76  field  year,  the  following  detailed 
changes  were  observed: 

(a)  In  October  1975,  an  advanced  stage  of  bar  attachment  to 
the  beach  face  was  evidenced  by  the  presence  of  two  large  inter- 
tidal ridges  only  partly  separated  from  the  beach  face  by  a 
narrow  and  short  runnel  (Fig.  35,  point  B) . Immediately  north 
of  the  point  of  bar  attachment  (Fig.  35,  point  A),  the  beach 
face  is  very  steep. 

(b)  By  January  1976,  the  runnel  had  been  closed  off  by 
accretion  (Fig.  36,  point  C)  forming  an  isolated  saltwater  pond 
(point  B) . Rapid  deposition  had  caused  significant  accretion 
and  seaward  displacement  of  the  berm  crest  off  point  B.  The 
second  intertidal  ridge  (Fig.  37,  point  D)  had  become  morpho- 
logically more  distinct  and  undergone  significant  vertical 
accretion  and  landward  migration  as  illustrated  by  the  topo- 
graphic maps  in  Figures  35  and  36. 

(c)  By  March  1976  (Fig.  37),  the  topography  of  the  site  of 
the  swash  bar  attachment  had  become  similar  to  that  in  October 
1975,  the  primary  difference  being  that  all  features  were  dis- 
placed about  50  meters  seaward  of  their  October  position. 

During  this  6-month  period,  one  episode  of  ridge  welding  had 
been  carried  to  completion. 

e.  Beaches . The  Debidue  Island  shore  in  the  vicinity  of  the  inlet 
consists  of  a recurved  spit  experiencing  seaward  progradation  and  land- 
ward retreat  in  response  to  migrating  ridge-and- runnel  systems.  As 
these  systems  migrate  shoreward,  sediment  is  moved  alongshore  in  dis- 
crete packages,  and  some  is  delivered  to  the  inlet  throat  by  currents 
in  the  marginal  flood  channel.  This  process  has  been  well  documented 
also  at  other  inlets  (Farrell,  1969;  Kumar  and  Sanders,  1974).  Sur- 
veys in  June  1972  and  July  1974  demonstrated  landward  retreat  of  the 
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A.  September  1974 


1 OCT.  1975 


Oblique  air  photo  and  topographic  map  of  swash  bar  complex  off  North  Island, 

1 October  1975.  Late  stage  swash  bar  migration  and  beach  attachment  are  charac 
terized  by  steep  beach  face  (A)  and  confined  runnel  (B) . (Datum  is  MSL. ) 
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beach  accretion  at  C,  ponding  at  B,  and  the  onshore  mig 
ridge  at  D.  (Datum  is  MSL.) 


Debidue  beach  (Finley,  1976).  During  these  2 years,  the  base  of  the 
beach  face  had  retreated  about  60  meters  along  a 300-meter  section. 
Although  the  volume  of  sand  was  not  computed,  a large  amount  of  dep- 
osition probably  occurred  on  the  inlet  shoals.  More  recently,  between 
September  1975  and  March  1976,  a cycle  of  deposition  and  erosion  was 
documented  at  the  recurved  spit.  The  area  monitored  in  this  recent 
survey  was  located  to  the  west  of  the  1972-74  site  at  the  landward 
extension  of  the  inlet  throat  profile  (Fig.  28,  DBI-TS).  The  follow- 
ing is  a discussion  of  variations  in  spit  morphology  at  this  site, 
referencing  the  air  photos  and  topographic  maps  in  Figures  38,  39, 
and  40: 


(a)  Figure  38,  28  September  1975:  Along  the  western  margin 
of  the  mapped  spit  area,  a steep,  irregular  topography  is  illus- 
trated by  the  spacing  and  curvature  of  the  contours.  A reentrant 
at  A defines  the  western  boundary  of  the  lobate  part  of  the  spit. 

A small,  partly  enclosed  runnel  at  B drained  toward  A.  A well- 
defined  berm  crest  existed  only  west  of  B.  Broad  washovers 
characterized  the  rest  of  the  survey  area  and  continued  east  of 
profile  DBI-TS.  Seaward  of  this  broad  overwash  was  a relatively 
steep  beach  face  becoming  gradually  more  gentle  toward  A. 

(b)  Figure  39,  11  January  1976:  The  survey  area  had  under- 
gone extensive  deposition  since  September.  Beach  accretion  had 
displaced  the  berm  crest  toward  the  inlet.  A uniform  beach- 
face  topography  characterized  the  arcuate  margin  of  the  spit. 

The  reentrant  at  A had  changed  into  a narrow  runnel  outlet 
draining  a pond  at  B. 

(c)  Figure  40,  19  March  1976:  Beach-face  erosion  charac- 
terized the  entire  survey  area.  About  25  meters  of  shoreline 
recession  had  occurred  since  January,  leaving  a uniformly 
sloping  beach  face  from  A around  the  recurve  to  profile  DBI-TS. 

The  berm  crest  along  this  margin  had  receded  to  within  10 
meters  of  area  B.  Minor  shoreline  modification  had  occurred 
west  of  area  A.  Little  evidence  of  the  outlet  channel  between 
A and  B was  preserved,  and  B was  a dry,  isolated  depression. 

The  two  ocean-facing  beaches  adjacent  to  North  Inlet  have  very  simi- 
lar intertidal  morphology  and  quite  similar  response  characteristics  on 
a month-by-month  basis.  Beach  morphology  was  monitored  throughout  the 
year  at  nine  profile  locations  at  Debidue  and  North  Islands  (Fig.  28). 
The  individual  profiles  are  included  in  Appendix  I. 

From  the  profiles,  the  amounts  of  erosion  and  deposition  between 
consecutive  surveys  were  computed  (expressed  as  cubic  meters  of  sand 
added  or  removed  per  running  meter  of  beach) . These  results  were 
plotted  as  a function  of  time  using  the  beach  profile  for  June  1975 
as  a reference  datum  (Fig.  41). 

Profile  DBI-25,  which  is  the  one  farthest  away  from  the  inlet, 
is  quite  stable.  Most  other  profiles  recorded  significant  amounts 
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Figure  40.  Oblique  air  photo  and  topographic  map  of  Debidue  Island  recurved 
spit,  19  March  1976.  Letters  indicate  the  same  shoreline 
features  as  in  Figure  35.  (Datum  is  MSL. ) 


Figure  41.  Trends  of  beach  erosion  and  deposition  on  Debidue  and  North  Islands  from 

June  197S  to  May  1976.  Faint  dotted  lines  indicate  missing  data.  Positive 
values  designate  deposition;  negative  values  indicate  erosion. 


of  erosion  and  deposition  during  the  year.  Comparisons  of  the  fluctu- 
ations in  all  profiles  on  Debidue  Island  show  that  variations  are  in 
phase. 

From  October  to  December  1975,  relative  stability  was  reflected 
by  nearly  all  DBI  profiles.  For  the  remainder  of  the  year,  three  dis- 
tinct cycles  of  erosion  and  deposition  were  depicted  by  falling  and 
rising  curves.  One  cycle  occurred  between  late  August  and  late  Octo- 
ber 1975.  Two  cycles  occurred  from  late  December  1975  through  mid- 
May  1976. 

The  same  cycles  are  also  expressed  by  the  beach  profiles  on  North 
Island.  Although  data  were  missing  for  late  September  to  late  Decem- 
ber 1975  at  North  Island,  enough  data  did  exist  for  the  early  fall  to 
indicate  beach  behavior  similar  to  that  on  Debidue  Island.  The  data 
demonstrate  quite  convincingly  that  beach  erosion  or  deposition  occurs 
in  phase  on  the  two  ocean- facing  beaches  adjacent  to  the  inlet. 

4.  Process  Response  Interactions. 

The  morphological  effects  of  seasonal  changes  in  MSL  are  most  evident 
at  the  Debidue  Island  recurved  spit  and  at  the  point  of  swash  bar  attach- 
ment on  North  Island. 

Shoreline  positions  at  the  Dehidue  spit  for  the  three  survey  times 
during  1975-76  are  shown  in  Figure  4£.  The  September  shoreline  is 
landwardmost  and  also  characterized  by  a broad  reentrant,  extensive 
washovers,  and  a steep  beach  face  (Fig.  38).  The  survey,  which  was 
made  at  maximum  MSL  (Fig.  18),  strongly  suggests  a direct  link  be- 
tween sea  level  elevation,  inlet  tidal  prism,  and  throat  cross-sec- 
tional area.  The  annual  changes  in  channel  flow  areas  (Table  4) 
were  consistent  with  this  suggestion. 

By  January,  the  Debidue  spit  accreted  (Figs.  39  and  42)  as  sea 
level  dropped,  and  the  throat  cross  section  was  reduced.  In  March, 
near  the  peak  of  the  spring  high  in  sea  level,  the  entire  Debidue 
spit  was  characterized  by  erosion.  Annual  changes  in  beach-face  mor- 
phology of  the  Debidue  Island  spit  were  directly  controlled  by  MSL 
fluctuations,  rather  than  storm  cycles  or  other  external  events  (Hum- 
phries, 1977a,  1977b). 

A slightly  different  process  response  interaction  occurred  at  the 
swash  bar  attachment  on  North  Island.  The  actual  morphologic  change 
at  this  site  was  complicated  by  the  superposition  of  two  events — 
seasonal  changes  in  sea  level  and  bar  accretion.  The  October  beach 
face  had  a typical  slope  of  0.5°.  A 30-centimeter  reduction  in  MSL 
between  October  and  January  would,  therefore,  correspond  to  a 35-meter 
seaward  progradation  of  the  shoreline.  This  is  in  good  agreement 
with  the  observed  progradation  (Fig.  43).  In  the  same  time  interval, 
significant  vertical  accretion  occurred  on  the  beach  (Figs.  35  and 
36) --a  result  of  the  landward  migration  of  the  swash  bar  conqplex. 
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Figure  42.  Seasonal  change  in  position  of  the  shoreline  at  the  Debidue  Island  recurved  spit 
(based  on  Figs.  21,  22,  and  23).  Monthly  local  MLHf  outlines  the  shoreline. 


This  steepened  the  beach  face  to  about  1.2°  such  that  the  spring  sea 
level  rise  was  accompanied  by,  typically,  only  20  meters  of  shore- 
line transgression  (Fig.  43). 

Perhaps  the  most  significant  response  interaction  was  that  which 
existed  between  inlet  channel  cross  sections  and  the  volume  of  the 
adjacent  beaches.  Figure  44  shows  the  'average'  pattern  of  erosion  and 
deposition  on  the  Debidue  and  North  Island  beaches  (adapted  from  Fig.  41). 
Both  beach  volumes  and  channel  cross  sections  fluctuated  about  a stable 
annual  mean.  The  phase  relationship  between  the  fluctuations  throughout 
the  year  was  quite  striking.  Beach  accretion  (positive  slope)  corres- 
ponded rather  consistently  with  channel  erosion;  beach  erosion  (negative 
slope)  corresponded  to  channel  infilling.  (An  exception  to  this  trend 
appeared  to  exist  for  the  fall  months.  However,  the  North  Island  data 
for  this  time  period  were  rather  sparse,  and  interpretations  of  phase 
relationships  based  on  the  fall  trends  could  be  misleading.)  Although 
more  closely  spaced  data  points  would  have  been  desirable,  the  data 
strongly  suggest  an  out-of-phase  relationship  between  beach  erosion  and 
channel  erosion.  The  data  provided  only  speculation  concerning  the 
mechanisms  of  sediment  exchange  which  tie  inlet  channel  to  beach  mor- 
phology. A feasible  chain  of  events  is  the  following:  Storm-induced 
beach  erosion  would  be  associated  with  the  delivery  of  an  increased 
sediment  load  to  the  inlet  channels,  thus  reducing  their  flow  areas. 
Subsequent  decrease  in  wave  action,  in  turn,  would  permit  channel  scour. 
This  model  is  consistent  with  Escoffier's  (1940)  and  Inman  and  Frautschy's 
(1960)  inlet  stability  concepts. 

VI.  SUMMARY  AND  CONCLUSIONS 

Both  this  report  and  Finley  (1976)  analyzed  the  basic  hydraulic, 
morphologic,  and  littoral  process  characteristics  of  North  Inlet, 

South  Carolina.  Data  were  obtained  over  a 2-year  period  from  July 
1974  to  May  1976.  Water  level  records  from  one  ocean  and  two  bay 
tide  gages,  combined  with  tidal  channel  hydrographic  data,  littoral 
process  measurements,  intertidal  morphological  maps,  bathymetric 
profiles,  and  aerial  photographic  records  (all  obtained  on  a quar- 
terly basis)  provided  the  main  data  base  for  the  study.  Channel 
cross  sections  were  monitored  on  a monthly  basis. 

North  Inlet  is  located  on  the  northern  South  Carolina  coast,  in  a 
transition  zone  between  the  wave-dominated  coast  of  North  Carolina  and 
the  distinctly  tide-dominated  coast  of  southern  South  Carolina  and 
Georgia.  Although  the  mean  ocean  tidal  range  is  only  167  centimeters. 
North  Inlet  is  classified  as  mesotidal  because  tidal  current-generated 
features,  conforming  to  those  of  existing  mesotidal  inlet  models,  domi- 
nate the  system.  The  low  annual  wave  energy  and  the  infrequent  storm 
passage  across  the  South  Carolina  coast  appear  to  permit  the  development 
of  a typical  mesotidal  morphology  for  a tidal  range  slightly  below  the 
proposed  worldwide  mesotidal  interval. 
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1975  1976 

Figure  44.  Average  trends  in  beach  erosion  and  deposition  at  Debidue  and  North  Islands  (top) 
and  channel  cross  sections  (bottom)  (derived  from  data  in  Figs.  24  and  41). 


Within  the  same  seasons  of  the  2 consecutive  years  of  observa- 
tion, the  weather  and  littoral  process  parameters  were  quite  com- 
parable. This  justifies  one  of  the  implicit  assumptions  in  the  study 
methodology:  the  sampling  of  process  parameters  on  a quarterly  basis 
over  2 years  will  provide  a representative  measure  of  the  total  an- 
nual variability  of  the  North  Inlet  environment. 

North  Inlet  is  hydraulically  ebb  dominated.  For  the  throat  section, 
the  peak  ebb  velocity  exceeded  the  peak  flood  by  a factor  of  1.22.  Also, 
the  ebb  duration  was  shorter  than  the  flood  duration  by  an  average  of 
26  minutes.  A model  is  proposed,  which  explains  the  ebb  dominance  as  a 
result  of  the  different  efficiency  of  water  exchange  between  the  ocean 
and  the  bay  at  high  and  low  tide.  Considering  continuity,  the  existence 
of  a much  larger  water  surface  area  in  the  bay  at  high  tide  rather  than 
at  low  tide  results  in  a longer  lag  at  high  water  than  at  low  water,  an 
average  of  46  versus  33  minutes.  This  differential  lag  is  consistent 
with  the  observed  asymmetries  in  the  duration  and  velocity  of  the  flood 
and  ebb. 


Spectral  analysis  of  the  tide  records  further  verify  the  model. 
Overtides,  M4,  M6,  and  Mg  components,  are  present  in  all  the  tide 
records.  These  are  produced  by  propagation  of  the  tidal  wave  in 
shallow  water.  However,  their  amplitudes  are  increased  by  one  to 
two  orders  of  magnitude  as  the  tidal  wave  moved  in  through  the  in- 
let; this  is  due  to  the  presence  of  nonlinear  terms  in  the  equations 
of  motion  for  inlet  flow.  The  development  of  an  asymmetric  bay  tid- 
al curve  relative  to  the  one  in  the  ocean,  is  consistent  with  the 
amplification  of  the  overtides. 

Average  surf  zone  conditions  expressed  in  terms  of  breaker  types, 
breaker  heights,  and  longshore  currents  exhibited  a consistent  season- 
al pattern  from  one  year  to  the  next.  Mean  breaker  heights  ranged 
from  a low  of  54  centimeters  in  March  to  a high  of  70  centimeters  in 
September.  The  longshore  wave  energy  flux  factor  was  computed  from 
the  littoral  process  measurements  and  used  to  estimate  rates  of  long- 
shore sediment  transport.  The  gross  transport  was  800,000  cubic 
meters  per  year  for  both  years  of  observation.  The  net  transport  for 
1974-75  and  1975-76  was  87,000  and  390,000  cubic  meters,  respectively, 
both  years  to  the  south.  This  difference  is  also  reflected  in  the  an- 
nual resultant  wind  vectors  derived  from  anemometer  readings  at  Debi- 
due  Island  during  the  quarterly  field  periods. 

Literature  reviews  suggest  that  the  role  of  wind  stress  in  the 
generation  of  longshore  currents  is  largely  unknown.  From  the  250 
complete  data  sets  of  littoral  environmental  observations  for  Debidue 
Island,  stepwise  regression  analyses  were  performed  to  isolate  the  ex- 
ternal variable  that  would  explain  the  largest  amount  of  variance  in 
the  observed  longshore  current  velocities.  In  all  regression  analyses, 
the  longshore  component  of  the  wind  velocity  was  the  parameter  ex- 
plaining most  of  the  variance.  Other  littoral  parameters,  including 


breaker  angle  and  breaker  height,  used  in  previously  developed  theo- 
retical or  semierapirical  longshore  current  velocity  equations  all 
explained  far  less  of  the  observed  variance.  It  appears  that  the 
longshore  component  of  wind  stress  is  at  leas1-  as  important  as  the 
momentum  of  the  obliquely  breaking  wave  in  the  generation  of  currents 
in  the  surf  zone  at  Debidue  Island,  and  probably  in  similar  coastal 
envi ronments . 

Analysis  of  the  tide  gage  records  demonstrated  a distinct  pat- 
tern of  seasonal  variations  in  MSL.  The  measured  bay  variation  was 
34  centimeters,  the  inferred  ocean  variation  about  45  centimeters. 

The  annual  high  occurred  in  mid- fall,  the  low  in  January,  and  a 
secondary  high  occurred  in  mid-spring.  A combination  of  steric 
effects,  seasonal  wind  setups  and,  perhaps,  river  discharge,  appears 
adequate  to  explain  the  observed  sea  level  fluctuations. 

Repetitive  topographic  mapping  of  the  intertidal  channel-margin 
linear  bars,  the  Debidue  Island  recurved  spit,  and  the  swash  bar  com- 
plex welding  on  to  North  Island  revealed  a distinct  seasonal  pattern 
in  the  morphologic  change.  It  is  suggested  that  this  morphologic  pat- 
tern is,  at  least  in  part,  a response  to  the  seasonal  changes  in  MSL. 
For  example,  the  Debidue  spit  shoreline  occupied  its  landwardmost 
position  during  the  September  1975  survey.  In  January  1976,  the  spit 
shore  had  accreted  several  tens  of  meters.  In  March,  it  had  retreated 
again.  This  suggests  that  the  Debidue  Spit  eroded  during  periods  of 
high  sea  level  and  consequent  large  tidal  prism  in  North  Inlet,  and 
accreted  during  periods  of  low  sea  levels. 

Existing  concepts  of  inlet  stability  suggest  that  there  is  an  ex- 
change of  sediment  between  the  inlet  channels  and  the  adjacent  beaches. 
Thus,  during  periods  of  high  wave  action,  when  the  beaches  will  under- 
go erosion  and  retreat,  sedimentation  and  shoaling  should  be  evident 
in  the  tidal  channels.  During  periods  of  fair  weather,  the  beaches 
will  accrete  and  the  tidal  creeks  scour.  Since  MSL  is  not  a constant, 
the  effects  of  this  sediment  exchange  pattern  could  be  masked  by 
changes  due  to  water  level  fluctuations.  Nevertheless,  in  the  winter 
and  spring  of  1976,  a clear  out-of-phase  relationship  between  beach 
erosion  and  channel  scour  was  detected,  in  good  agreement  with  the 
sediment  exchange  mechanism  stated  above. 
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FATHOMETER  PROFILES  AT  NORTH  INLET 

°fllfs,l'?re  “easured  by  a Bludworth  precision  dept 

rZ*u’  A ”ap  of  channel  profile  locations  is 
onshore  profile  locations  are  shown  in  Figure  28 
computer  reduced  to  the  same  scale. 
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Figure  A- 4.  Jones  Creek,  6 January  to  19  May  1976. 
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Figure  A-6.  Flood  tidal  delta,  6 January  to  20  March  1976 
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Figure  A- 7.  Inlet  throat,  13  June  to  24  November  1975. 
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Stations  NI-6,  June  1975  to  19  March  1976;  DBI-6, 

20  June  1975  and  18  March  1976;  and  DBI-8,  19  June  1975. 
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APPENDIX  B 

TIME-SERIES  DATA  ON  CURRENT  VELOCITIES  AND 
WATER  SURFACE  DIFFERENTIALS 

The  velocities  are  integrated  means  of  vertical  current  profiles  at 
the  indicated  stations  (Figs.  B-l  to  B-15).  Station  locations  are  shown 
in  Figure  9.  All  hydrography  measurements  are  listed  in  Table  5.  The 
water  surface  differentials  are  determined  as  follows:  Town  Creek  and 
Inlet  Throat  differentials  refer  to  the  instantaneous  difference  in  water 
elevation  at  the  ocean  and  the  Town  Creek  gages  (for  gage  location,  see 
Fig. 7);  a positive  difference  indicates  a landward-dipping  slope.  Jones 
Creek  differentials  are  similarly  determined  from  the  ocean  and  Jones 
Creek  gage  records. 

Hi  and  Lo  indicate  times  of  high  and  low  tides  in  the  ocean. 
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CURRENT  VELOCITIES 


Figure  B- 


INLET  THROAT 
7 January  1976 
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. Current  velocity  time  series  and  surface  differential 
at  the  inlet  throat,  7 January  1976. 
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Figure  B-15.  Current  velocity  time  series  at  the  inlet  throat,  2 June 
2 October  1975,  and  17  February  1976,  and  at  Town  and 
Jones  Creeks  24  November  1975. 


APPENDIX  C 

TIDAL  CURRENT  VELOCITY  AND  PHASE  DATA,  1975-76  SEASON 

This  appendix  presents  summary  tables  of  raw  and  reduced  data  per- 
taining to  the  hydrographic  measurements  made  at  North  Inlet  in  the  1975- 
76  field  season  (Tables  C-l  to  C-ll).  Tabulated  frequency  distributions 
and  histograms  are  included  for  all  the  variables.  For  station  locations, 
see  Figure  9.  The  variables  are  defined  as  follows: 

VFMAX:  Peak  flood  velocity  (cm/s) 

VEMAX:  Peak  ebb  velocity  (cm/s) 

DELTF:  Time  between  peak  flood  velocity  and  subsequent 
high  water  in  the  ocean  (hr). 

DELTE:  Time  between  high  water  in  the  ocean  and  sub- 
sequent peak  ebb  velocity  (hr) . 

LAGH:  Time  between  ocean  high  tide  and  subsequent  sta- 
tion slack  (hr). 

LAGL:  Time  between  ocean  low  tide  and  subsequent  sta- 
tion slack  (hr) . 
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Table  C-l.  Measured  hydraulic  parameters  at  North  Inlet. 
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Table  C-l.  Measured  hydraulic  parameters  at  North  Inlet. — 
Continued 
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DATE 
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deltf 

DELTF 
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0.9 

• 

9b 

TCS 

92b7b 

62 

71 

1.7 

3.6 
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Table  C-2.  Frequency  distribution  of  maximum  flood-current 
velocity  (cm/s). 


VfHAX 

FREQUENCY 

CJv  FREQ 

PERCENT 

CUM  PERCENT 
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1.087 

1.087 

25 

1 

2 
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5.635 
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1 
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1 
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16 

6.368 

17.391 

52 

1 

17 

1.087 

18.679 

56 

1 

18 

1.087 

19.565 

56 

6 

22 

6,368 

23.913 

59 

2 

26 

2.176 

26.087 

60 

3 

27 

3.261 

29.368 

61 

2 

29 

2.176 

31.522 

62 

3 

32 

3.261 

36.783 

63 

2 

36 

2.176 

36.957 

66 

2 

36 

2.176 

39.130 

65 

3 

39 

3.261 

62.391 

66 

5 

66 

5.636 

67.826 

67 

2 

66 

2.176 

50.000 

68 

2 

68 

2.176 

52.176 

69 

2 

50 

2.17*. 

56.368 

71 

l 

51 

1.087 

55.635 

72 

2 

53 

2.176 

57.609 

76 

3 

56 

3.261 

60.870 

75 

2 

58 

2.176 

63.063 

76 

2 

60 

2.176 

65.217 

77 

2 

62 

2.176 

67.391 

78 

5 

67 

5.635 

72.826 

80 

2 

69 

2.176 

75.000 

81 

1 

70 

1.087 

76.087 

82 

3 

73 

3.261 

79.368 

66 

1 

76 

1.087 

80.635 

86 

1 

75 

1.087 

81 .522 

90 

1 

76 

1.087 

82.609 

92 

1 

77 

1.087 

83.696 

96 

1 

78 

1.087 

86.783 

99 

1 

79 

1.087 

85.870 

101 

3 

82 

3.261 

89.130 

106 

1 

83 

1.087 

90.217 

105 

1 

86 

1.087 

91.306 

109 

2 

86 

2.176 

93.678 

116 

2 

88 

2.176 

95.652 

117 

1 

89 

1.087 

96.739 

12J 

1 

90 

1.087 

97.626 

125 

2 

92 

2.176 

100.000 
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Table  C-3.  Frequency  distribution  of  maximum  ebb-current 
velocity  (cm/s). 


vemax 

FREQUENCY 

CUM  freQ 

PERCENT 

CUM  percent 

« 

9 

• 

• 

• 

30 

1 

1 

1.099 

1.099 

33 

1 

2 

1.099 

2.198 

36 

1 

3 

1.099 

3.297 

37 

1 

A 

1.099 

A. 396 

38 

1 

5 

1.099 

5.A95 

39 

1 

6 

1.099 

6.593 

A3 

1 

7 

1.099 

7.692 

52 

1 

6 

1.099 

8.791 

53 

1 

9 

1.099 

9.890 

5A 

3 

11 

2.198 

12.008 

57 

1 

12 

1.099 

13.187 

59 

1 

13 

1.099 

1A.2B6 

60 

2 

15 

2.198 

16.A8A 

63 

1 

16 

1.099 

17.582 

63 

3 

19 

3.297 

20.879 

6a 

1 

30 

1.099 

21.978 

65 

1 

21 

1.099 

33.077 

66 

3 

3A 

3.297 

26.37A 

66 

2 

26 

2.198 

28.571 

70 

2 

28 

2.196 

30.760 

71 

3 

31 

3.297 

3A.066 

73 

? 

33 

2.198 

36.36A 

73 

1 

3A 

1.099 

37.363 

TA 

2 

36 

2.198 

39.560 

76 

3 

39 

3.297 

02.657 

7 7 

1 

AO 

1.099 

A3. 956 

>8 

1 

A l 

1.099 

AS. 055 

79 

3 

A3 

2.198 

A7.253 

00 

5 

AB 

S.A95 

52.7A7 

81 

1 

A9 

1.099 

53.8A6 

83 

1 

50 

1.099 

5A.9A5 

83 

1 

51 

1.099 

S6.0AA 

8 A 

3 

53 

2.198 

58.2A2 

85 

3 

56 

3.397 

61.538 

86 

A 

60 

A. 396 

65,93a 

87 

1 

61 

1.099 

67.033 

89 

2 

63 

2.190 

69.331 

90 

1 

6A 

1.099 

70.330 

91 

1 

65 

1.099 

71  ,A?9 

93 

1 

66 

1.099 

72.527 

95 

2 

68 

2.198 

7 A . 725 

100 

1 

69 

1.099 

75.02A 

103 

2 

71 

2.198 

78.022 

106 

1 

72 

1.099 

79.121 

ioa 

A 

76 

A. 396 

■83.516 

109 

1 

77 

1.099 

HA. 615 

113 

1 

78 

1.099 

85. 71 A 

11<* 

2 

80 

2.198 

87.91? 

130 

1 

81 

1.099 

89.011 

131 

1 

82 

1.099 

90.110 

133 

2 

8a 

2.198 

92.308 

136 

1 

85 

1.099 

93.A07 

137 

1 

86 

1.099 

9A.505 

130 

1 

87 

1.099 

95.60A 

133 

3 

89 

2.198 

97.802 

1 aO 

l 

90 

1.099 

90.001 

150 

1 

91 

1.099 

100.000 
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Table  C-4.  Frequency  distribution  of  the  time  difference 

between  peak  flood  velocity  (at  all  hydrography 
stations)  and  high  water  in  the  ocean. 


Otl.Tr 

rxzoufNcr 

cum  raco 

acwcCNi 

CUM  mcqccnt 

e 

8 

e 

« 

. 

0 

1 

1 

1.067 

1.087 

0.2 

1 

2 

1.087 

2.174 

0*3 

1 

3 

1.087 

3*261 

o.<* 

3 

6 

3*261 

6.52? 

OeS 

7 

13 

7.609 

14.130 

0.6 

2 

IS 

2.174 

16.304 

0.7 

3 

18 

3.261 

19.565 

0.9 

4 

22 

4.348 

23.913 

1 

12 

34 

13.043 

36.957 

lei 

5 

39 

5.435 

42.391 

1 .2 

4 

43 

4.348 

46. 739 

lei 

7 

SO 

7.609 

54.346 

1.4 

6 

56 

6.522 

60.270 

l.S 

2 

58 

2.174 

63.043 

1.6 

63 

5.435 

46.476 

1.7 

6 

71 

6.696 

77.174 

1.1 

4 

75 

4.346 

61.52? 

1*9 

4 

79 

4.348 

65.870 

2 

2 

61 

2.174 

86.043 

2.1 

2 

83 

2.174 

90.217 

2.2 

1 

64 

1.067 

91.304 

2.3 

3 

67 

3.261 

94.565 

2.* 

1 

88 

1.067 

95.65? 

2.5 

1 

89 

1.067 

96.739 

2.8 

1 

90 

1.067 

97.826 

2.9 

1 

91 

1.067 

98.913 

1 

1 

1.067 

100.000 

Table  C-5.  Frequency  distribution  of  the  time  difference 
between  high  water  in  the  ocean  and  the  sub- 
sequent peak  ebb  velocity  (at  all  hydrography 
stations) . 


ocltc 

fkEOUCNCf 

com  r«Eo 

nRCCNT 

CUM  M€»CCNT 

• 

16 

• 

• 

• 

0.9 

1 

l 

1.216 

1.220 

1.3 

1 

2 

1.220 

2.439 

1.7 

2 

* 

2.439 

4,678 

1.9 

1 

5 

1.220 

4.096 

2 

1 

« 

1.220 

7.317 

2.1 

1 

T 

1.220 

6.537 

2.2 

3 

10 

3.659 

12.195 

2.3 

5 

IS 

6.096 

16.293 

2.4 

3 

»• 

3.659 

21.951 

2.5 

1 

J* 

1.220 

23.171 

2.6 

3 

22 

3.659 

26.629 

2.7 

3 

2S 

3.659 

30,466 

2.6 

3 

2« 

3.659 

34.146 

2.9 

6 

J* 

7.317 

41.443 

3 

5 

J» 

6.099 

47,561 

3.1 

4 

.1 

4.678 

52.439 

3.2 

7 

so 

8.537 

40.976 

3.3 

3 

S3 

3.659 

44.634 

3.4 

2 

ss 

2.439 

47.073 

3.5 

6 

*1 

7.317 

74.340 

3.6 

6 

*7 

7.317 

81,707 

3.7 

2 

♦» 

2.439 

•4.|46 

3.6 

1 

T» 

1.126 

•9.346 

3.9 

f 

r 2 

2.439 

•7.605 

4 

3 

TS 

3.659 

91.463 

4.1 

3 

7* 

3.659 

95.11? 

4.2 

2 

SO 

2.439 

47.461 

4.5 

1 

• 1 

1.220 

46.760 

4.6 

1 

•2 

1.220 

100.000 
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Table  C-6.  Frequency  distribution  of  the  time  lag  between 
ocean  high  tide  and  subsequent  station  slack. 


LAGM 

FREQUENCY 

CUM  FREQ 

PERCENT 

CUM  PERCENT 

• 

11 

. 

• 

0 

1 

1 

1.126 

1.126 

0.2 

6 

6 

5.618 

6.782 

0.3 

3 

9 

3.371 

10.112 

0.* 

1C 

19 

11.236 

21.368 

o.s 

3 

22 

3.371 

26.719 

0.6 

9 

31 

10.112 

36.831 

0.7 

5 

36 

5.618 

60.669 

0.6 

10 

66 

11.236 

51.685 

0.9 

12 

58 

13.633 

65.169 

1 

13 

71 

16.607 

79.715 

1.1 

6 

77 

6.762 

86.517 

1.2 

3 

60 

3.371 

89.668 

1.3 

2 

82 

2.267 

92.135 

1.6 

3 

85 

3.371 

95.506 

l.S 

1 

86 

1.126 

96.629 

1.7 

1 

87 

1.126 

97.753 

1.6 

1 

68 

1.126 

98.876 

1.9 

1 

89 

1.126 

100.000 

Table  C-7.  Frequency  distribution  of  the  time  lag  between 
ocean  low  tide  and  subsequent  station  slack. 


LAGL 

FRtQUENCY 

CUM  FREQ 

PERCENT 

CUM  PERCENT 

• 

62 

• 

• 

• 

0 

3 

3 

5.172 

5.172 

0.1 

6 

7 

6.897 

12.069 

0.2 

6 

13 

10.365 

22.616 

0.3 

3 

16 

5.172 

27.586 

0.6 

7 

23 

12.069 

39.655 

O.S 

11 

36 

18.966 

58.621 

0.6 

11 

65 

16.966 

77.586 

0.7 

3 

68 

5.172 

82.759 

0.8 

3 

51 

5.172 

87.931 

0.9 

1 

52 

1.726 

89.655 

1 

2 

56 

3.668 

93.103 

1.2 

3 

57 

5.172 

98.276 

1.6 

1 

58 

1.726 

100.000 
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Table  C-8.  Histograms  and  summary  statistics  for  variables 
VFMAX  and  VEMAX. 
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Table  C-9.  Histograms  and  summary  statistics  for  variables 
DELTF  and  DELTE. 
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Table  C-10.  Histograms  and  summary  statistics  for  variables 
LAGH  and  LAGL. 
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Table  C-ll.  Hydraulic  parameters  of  the  North  Inlet  system 
summarized  by  station. 


APPENDIX  D 


TIDAL  CURRENT  VELOCITY  AND  PHASE  DATA,  1974-75  SEASON 

This  appendix  presents  summary  tables  of  raw  and  reduced  data  per- 
taining to  the  hydrographic  measurements  made  at  North  Inlet  in  the 
1974-75  field  season.  Tabulated  frequency  distributions  and  histograms 
are  included  for  all  the  variables  (Tables  D-l  to  D-ll).  For  station 
locations,  see  Figure  9.  The  variables  are  defined  as  follows: 

VFMAX:  Peak  flood  velocity  (cm/s) 

VEMAX:  Peak  ebb  velocity  (cm/s) 

DELTF : Time  hetween  peak  flood  velocity  and  subsequent 
high  water  in  the  ocean  (hr) . 

DELTE:  Time  between  high  water  in  the  ocean  and  subse- 
quent peak  ebb  velocity  (hr) . 

LAGH:  Time  between  ocean  high  tide  and  subsequent  sta- 
tion slack  (hr). 

LAGL:  Time  between  ocean  low  tide  and  subsequent  station 
slack  (hr). 


157 


i 


! 


fBXS  PAGE  IS  BEST  QUALITY  PRA<3IIfijW|, 

IBM!  fltyy  lUwuwu  xv  JOiiQ 

Table  D-l.  Measured  hydraulic  parameters  at  North  Inlet. 


oas 

STATION 

DATE 

VFMAX 

VEMAX 

DELTF 

OELTE 

LAGH 

LAGL 

1 

ITN 

91974 

no 

156 

2.1 

3.1 

0.8 

1.0 

2 

itn 

32175 

46 

91 

2.0 

3.6 

0.8 

• 

3 

itn 

42675 

92 

126 

1.1 

4.4 

1.7 

0.3 

4 

its 

91974 

111 

118 

1.7 

3.6 

0.7 

1.1 

5 

ITS 

32175 

49 

78 

1.9 

3.2 

1.2 

• 

6 

ITS 

42675 

103 

103 

0.7 

4.2 

0.9 

0.3 

? 

TCN1 

82474 

65 

91 

1.7 

3.3 

0.8 

0.8 

a 

TCNl 

91574 

88 

105 

• 

• 

• 

• 

9 

TCN1 

92174 

50 

104 

2.2 

3.3 

0.6 

• 

10 

TCNl 

92774 

62 

• 

2.9 

. 

• 

0.6 

11 

TCNl 

22275 

61 

99 

1.4 

3.9 

0.9 

• 

12 

TCNl 

31975 

42 

81 

l.S 

3.2 

0.4 

• 

13 

TCNl 

12275 

47 

82 

2.2 

3.2 

0.4 

• 

14 

TCNl 

32475 

84 

108 

1.1 

» 

0.7 

0.7 

IS 

TCN2 

02474 

90 

96 

1.6 

3.7 

1.1 

• 

16 

TCN2 

91574 

128 

88 

• 

. 

• 

• 

1? 

TCN2 

92174 

81 

96 

1.4 

3.2 

0.8 

• 

10 

TCN2 

92774 

77 

• 

2.7 

. 

• 

0.6 

19 

TCN2 

22275 

67 

81 

0.7 

2.6 

1.2 

• 

20 

TCN2 

31975 

60 

72 

1.9 

3.9 

0.8 

• 

21 

TCN2 

32275 

60 

69 

1.9 

3.0 

0.6 

t 

2? 

TCN2 

32475 

78 

100 

2.1 

. 

0.8 

1.1 

23 

TCN3 

02474 

72 

87 

1.5 

3.8 

1.2 

• 

24 

TCN3 

91574 

116 

96 

• 

• 

. 

• 

25 

TCN3 

92174 

81 

78 

0.7 

3.2 

0.8 

• * 

26 

TCN3 

92774 

73 

• 

1.8 

• 

• 

0.3 

27 

TCN3 

22275 

70 

76 

0.5 

4.2 

1.4 

• 

20 

TCN3 

31975 

59 

74 

2.0 

3.6 

0.4 

0.3 

29 

TCN3 

32275 

60 

65 

1.9 

3.2 

0.7 

• 

30 

TCN3 

32475 

75 

97 

2.0 

• 

0.9 

1.1 

31 

JCN1 

82474 

72 

84 

0.2 

3.5 

1.4 

1.7 

32 

JCN1 

92174 

77 

89 

0.6 

3.3 

1.0 

• 

33 

JCNl 

92774 

74 

• 

2.1 

• 

• 

0.7 

34 

JCNl 

22275 

68 

84 

0.4 

3.2 

2.0 

• 

35 

JCN1 

31975 

34 

76 

2.2 

3.6 

0.4 

0.4 

36 

JCNl 

32275 

36 

81 

1.1 

2.5 

0.1 

• 

37 

JCNl 

32475 

39 

92 

1.4 

• 

• 

1.8 

30 

JCN2 

02474 

58 

75 

0.7 

3.7 

1.0 

2.1 

39 

JCN2 

92174 

94 

78 

0.9 

3.2 

1.1 

• 

40 

JCN2 

927^4 

90 

• 

1.7 

• 

• 

0.5 

61 

JCN2 

22275 

91 

85 

1.1 

3.6 

1.0 

• 

42 

JCN2 

31975 

44 

65 

2.4 

3.8 
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Table  D-2.  Frequency  distribution  of  maximum  flood-current 
velocity  (cm/s). 


vfma< 

frequency 

CUM  FSEO 

PEHCENT 

CJM  PERCENT 

26 

1 

1 

2,222 

2.222 

34 

1 

2 

2.222 

4,444 

36 

1 

3 

2.222 

6.667 

39 

1 

4 

2.222 

0.009 

42 

1 

5 

2.222 

11.111 

44 

1 

6 

2.222 

13.333 

46 

1 

7 

2.222 

15.556 

47 

2 

9 

4.444 

20.000 

49 

1 

10 

2.222 

22.222 

50 

1 

11 

2.222 

24.444 

56 

1 

12 

2.222 

26.667 

59 

1 

13 

2.222 

20.009 

60 

3 

16 

6.667 

35.556 

61 

l 

IT 

2.222 

37.770 

62 

1 

10 

2.222 

40.000 

65 

1 

19 

2.222 

42.222 

67 

1 

20 

2.222 

44,444 

60 

1 

21 

2.222 

46.667 

70 

2 

23 

4,444 

51.111 

72 

2 

25 

4.444 

55 . 556 

73 

1 

26 

2.222 

57.770 

74 

1 

27 

2.222 

60.000 

75 

1 

20 

2.222 

62.222 

77 

2 

30 

4.444 

66.667 

70 

1 

31 

2.222 

60.009 

01 

2 

33 

4,444 

73.333 

04 

1 

34 

2.222 

75.556 

06 

1 

35 

2.222 

77,77a 

90 

2 

37 

4.444 

02.222 

91 

1 

30 

2.222 

04,444 

92 

1 

39 

2.222 

06.667 

94 

1 

40 

2.222 

00.009 

103 

1 

41 

2.222 

91. Ill 

no 

1 

4? 

2.222 

93.333 

in 

1 

43 

2.222 

95.556 

116 

1 

44 

2.222 

97.770 
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1 

45 

2.222 

100.000 
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Table  D-3.  Frequency  distribution  of  maximum  ebb-current 
velocity  (cm/s) . 

ve«A<  freuuen cy  cu«  ereq  percent  cu«  percent 
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• 

• 
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SR 

1 

1 

2.500 

2.500 

65 

2 

3 

5.000 

7.500 

69 

2 

5 

5.000 

12.500 

72 

1 

6 

2.500 

15.000 

74 

1 

7 

2.500 

17.500 

75 

1 

8 

2.500 

20.000 

76 

2 

10 

5.000 

25.000 

75 

3 

13 

7.500 

32.500 

81 

3 

16 

7.500 

40.000 

82 

1 

17 

2.500 

42,500 

83 

1 

18 

2.500 

45.000 

84 

2 

20 

5.000 

50.000 

85 

l 

21 
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52.500 

87 

1 

22 
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56.000 

88 

1 

23 
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57.500 

89 

1 

24 
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60.000 

91 

2 

26 
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65.000 
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1 

27 
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96 

3 

30 
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97 

1 

31 
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99 

1 
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1 

33 
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1 

34 

2.500 
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1 

35 
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1 

36 
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90.000 
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1 

37 

2.500 

92.500 

118 

1 

38 

2.500 

95.000 

125 

1 

39 

2.500 

97.500 

156 

J 

40 

2.500 

100, 000 

D-4. 

Frequency  distribution  of  the  time  difference 
between  peak  flood  velocity  (at  all  hydrography 
stations)  and  high  water  in  the  ocean. 

OELTE 

epeoufno 

CUv  EREQ  PERCENT 

Ct’M  t PERCENT 

• 

3 

• 

• 

. 

0.2 

l 

1 

2.381 

2.381 

0.4 

1 

2 

2.381 

4.762 

0.5 

1 

3 

2.381 

7.143 

0.6 

1 

4 

2.381 

9.524 

0.7 

4 

8 

9.524 

19.049 

0.9 

1 

9 

2.381 

21.429 

1.1 

5 

14 

11.905 

33.333 

1.4 

4 

18 

9.524 

42.857 

1.5 

2 

20 

4.762 

47.619 

1.6 

1 

21 

2.381 

50.000 

1.7 

3 

24 

7.143 

57.143 

1.5 

1 

25 
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59.524 

1.9 

4 

29 
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69.048 

2 

4 

33 
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78.571 

2.1 

3 

36 

7.143 

85.714 

2.2 

3 

39 

7.143 
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2.4 

1 
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2.381 
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2.7 

1 

41 

2.381 
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2.9 

1 

42 

2.381 

100.000 
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. Frequency  distribution  of  the  time  difference 
between  high  water  in  the  ocean  and  the  sub- 
sequent peak  ebb  velocity  (at  all  hydrography 


stations) . 

oecTe 

EREOUENCY 

cum  race 

PERCENT 

CUM  PERCENT 

e 

13 

e 

e 

e 

2.5 

1 

l 

3.125 

3.125 

2.6 

2 

3 

6.250 

9.375 

3 

1 

6 

3.125 

12.500 

3.1 

1 

5 
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15.625 

3.2 

6 

13 

25.000 
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3.3 

3 

16 
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50.000 

3.5 

1 

17 
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3.6 

5 

22 

15.625 
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2 

26 
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2 
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3.9 

3 

29 
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90.625 

6.2 

2 

31 
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96.875 

6.6 

1 

32 

3.125 

100.000 

D-6. 

Frequency  distribution  of  the  time  lag  between 
ocean  high  tide  and  subsequent  station  slack. 
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CUM  PERCENT 
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1 
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2 
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1 
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D-7. 

Frequency  distribution  of  the 

time  lag  between 

ocean  low  tide  and  subsequent 

station  slack. 
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Table  D-8.  Histograms  and  summary  statistics  for  variables 
VFMAX  and  VEMAX. 

STATISTICS  TABLE  ANO  HISTOGRAM  EO*  VEMAX  STATISTICS  TABLE  ANO  M15IUGWAM  f0R  vt»A' 

NO.  OE  values*  AS  no.  OE  VALUES*  40 

Suw*  31 7T  SUM*  3537 

UNCOBBECTfO  SS*247445  UNCOBRECTED  SS*3255m3 

VABIANCF*  526. 1 000900091  VARIANCE*  327. 789 1 02564 1 

COEE.  VAB1 ANCE *32.488764 94886  COEE.  VARIA5CE*20. 474918  75445 

NO.  missing*  0 *0.  MISSING*  5 

MS  AN*  70.60000000000  “EAN*  88.42500000000 

COBRECTEO  SS*  23148,80000000  CORBfCTEO  SS»  12783. 77500000 

STANOABO  0EVIATI0N*22. 03706905390  STANOABO  OEV I AT  ION* 1 8 . 1 0404600  862 

STANOABO  error*  11.69131313131  STANOABO  ERROR*  8.194727564103 
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Table  D-10.  Histograms  and  summary  statistics  for  variables 
LAGH  and  LAGL. 
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Table  0-11.  Hydraulic  parameters  of  the  North  Inlet  system 
summarized  by  station. 


APPENDIX  E 


HOURLY  CURRENT  VELOCITY  READINGS  AND 
WATER  SURFACE  DIFFERENTIALS 

This  appendix  presents  hourly  velocity  readings  (integrated  means 
through  a vertical  column)  at  all  North  Inlet  hydrography  stations  for  the 
1975-1976  field  season  (Tables  E-l  and  E-2).  The  hourly  water  surface 
differentials  were  obtained  from  simultaneous  tide  records  at  the  ocean, 
Jones  Creek, and  Town  Creek  gages  (Figs.  E-l,  E-2,  and  E-3) . The  numerical 
raw  data  in  this  appendix  correspond  to  the  graphical  presentation  in 
Appendix  B.  For  station  locations  see  Figure  9.  The  variables  are  de- 
fined as  follows : 

DATE:  Date  of  observation 

TIME:  Hour  of  observation  (e.s.t.) 

OCEAN:  Water  surface  elevation  at  ocean  gage  (in 
cm  relative  to  datum). 

TOWN:  Water  surface  elevation  at  Town  Creek  gage  (in 
cm  relative  to  datum) . 

JONES:  Water  surface  elevation  at  Jones  Creek  gage 
(in  cm  relative  to  datum). 

TCN1-FCN:  Current  velocity  at  the  representative 
observation  station. 

DELTOW:  Instantaneous  water  surface  elevation  differ- 
ence between  ocean  gage  and  Town  Creek  gage  (cm).  Positive 
differential  implies  landward  slope. 

DELJON:  Instantaneous  water  surface  elevation  differ- 
ential between  ocean  gage  and  Jones  Creek  gage  (cm). 

VELT : Mean  velocity  of  all  three  Town  Creek  stations 
(cm/s). 

VELJ:  Mean  velocity  of  both  Jones  Creek  stations  (cm/s). 

VELM:  Mean  velocity  of  all  three  inlet  throat  stations 
(cm/s). 


DIFT: 

Absolute 

value 

of 

DELTOW 

(cm) 

DIFJ: 

Absolute 

value 

of 

DELJON 

(cm) 

AVELT: 

Absolute 

value 

of 

VELT 

(cm/s) 

AVELJ : 

Absolute 

value 

of 

VELJ 

(cm/s) 
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Table  E- 

& J 

A • 

1 

Current  velocity  readings  (integrated  means)  at 
each  station  (cm/s),  water  surface  elevations 
(cm),  and  elevation  differentials  (cm). 
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Table  E-l.  Current  velocity  readings  (integrated  means)  at 
each  station  (cm/s),  water  surface  elevations 
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Table  E-2.  Correlation  matrix  for  current  velocity  observations. 
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SCATTERGRAM  OF  DELTOW 


STATISTICAL  PACKAGE  FOR  THE  SOCIAL  SCIENCES  SPSSH -RELEASE  6.02  11/24/76 

statistics: 

CORRELATION  (Rl-  0.77764  R SQUARED  -0.60472  SIGNIFICANCE  - 0.00001 

STD  ERR  OF  EST  -32.85753  INTERCEPT  (A)  -2.10095  SLOPE  (B)  - 2.61776 

PLOTTED  VALUES  -91  EXCLUDED  VALUES  - 0 MISSING  VALUES  - 1 1 7 

Figure  E-l.  Scatter  diagram  of  current  velocity  at  Town  Creek 
(VELT)  versus  the  instantaneous  water  surface 
elevation  differential  between  the  ocean  and  Town 
Creek  gages. 
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SCATTERGRAM  OF  OELJON 

STATISTICAL  PACKAGE  FOR  THE  SOCIAL  SCIENCES  SPSSH-RELEASE  6.02  11/24/76 

statistics: 

CORRELATION  IR)-  0.03751  R SOUAREO  -0.70142  SIGNIFICANCE  - 0.00001 

STD  ERR  OF  EST  -28.18429  INTERCEPT  (A)  -6.37400  SLOPE  IB)  - 2.04415 

PLOTTED  VALUES -89  EXCLUDED  VALUES  - 0 MISSING  VALUES  - 1 1 9 

Figure  E-2.  Scatter  diagram  of  current  velocity  at  Jones  Creek 
(VELJ)  versus  the  instananeous  wat^r  surface 
elevation  differential  between  the  ocean  and  Jones 
Creek  gages. 
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SCATTERGRAM  OF  DELTOW 


STATISTICAL  PACKAGE  FOR  THE  SOCIAL  SCIENCES  SPSSH  - RELEASE  6.02  11/24/76 

statistics: 

CORRELATION  (R)-  0.86393  R SQUARED  - 0.74637  SIGNIFICANCE  - 0.00001 

STD  ERR  OF  EST  -31. 8339-,  INTERCEPT  (A)  - 6.77610  SLOPE  (8)  - 3.84794 

PLOTTED  VALUES -64  EXCLUOED  VALUES  - 0 MISSING  VALUES  - 144 

Figure  E-3.  Scatter  diagram  of  current  velocity  at  the  inlet 

throat  (VELM)  versus  the  instantaneous  water  surface 
elevation  differential  between  the  ocean  and 
Town  Creek  gages. 
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APPENDIX  F 

WATER  SURFACE  ELEVATION  RECORDS 

These  tide  records  were  obtained  by  Stevens  and  Leupold  float  gages 
at  Town  and  Jones  Creeks,  and  by  a Bristol  blubber  gage  at  the  ocean 
(Figs.  F-l  to  F-16).  The  raw  strip-chart  records  have  been  digitized 
and  computer-reduced  to  a common  format.  The  elevations  are  related  to 
a common  datum,  which  is  205  centimeters  below  the  "INLET"  triangulation 
station.  The  available  records  are: 

Ocean  gage : 1975:  10-14  June,  16-28  June,  1-28  July, 

27  Aug. -3  Sept.,  9-23  Sept.,  2-15  Oct.,  18  Oct. -8  Nov., 

7-28  Dec.;  1976:  13-27  Mar.,  25  Apr.  -19  May. 

Town  Creek:  1975:  10-30  June,  28  July-21  Aug., 

30  Aug. -21  Sept.,  24  Oct. -30  Nov.,  30-31  Dec.;  1976: 

1 Jan. -28  Mar.,  31  Mar. -19  May. 

Jones  Creek:  1975:  9 June-24  July,  13-24  Sept., 

25  Oct. -31  Dec.;  1976:  1 Jan-16  Feb.,  18  Feb. -28  Mar., 

31  Mar. -17  May. 


NOTE: — The  Jones  Creek  record  for  the  period  of  1 Nov. -31  Dec.  1975 
is  included  as  Figure  12. 
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Figure  F-l.  Ocean  gage  water  surface  elevation  records,  June  and  July  1975. 


Figure  F-2.  Ocean  gage  water  surface  elevation  records,  August  and  September  1975. 


Figure  F-6.  Town  Creek  gage  water  surface  elevation  records,  June,  July  and  August  1975. 
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ion  records,  December  1975  and  March  1976. 


water  surface  elevation  records,  January  and  February  1976. 
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Figure  F-ll.  Town  Creek  gage  water  surface  elevation  records.  May  1976. 
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Figure  F-14.  Jones  Creek  gage  water  surface  elevation  records,  January  and  February  1976. 


Figure  F-15.  Jones  Creek  gage  water  surface  elevation  records,  March  and  April  1976. 


APPENDIX  G 

SPECTRA  OF  2-WEEK  WATER  SURFACE  ELEVATION  RECORDS 

Spectra  (Figs.  G-l  to  G-5)  were  calculated  by  applying  the  Statistical 
Analysis  System  (SAS)  "SPECTRA"  procedure  to  the  water  surface  elevation 
records  presented  in  Appendix  F. 
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Figure  6-1.  Spectra  of  the  ocean  and  Jones  Creek  gage  water  surface 
records,  9-23  July  197S. 
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Figure  G-3.  Spectra  of  the  ocean  and  Jones  Creek  gage  water  surface 
records,  8-23  December  1975. 
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Figure  G 
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FREQUENCY  t2Thr> 

2.  Spectra  of  the  ocean.  Town  Creek, and  Jones  Creek  gage 
water  surface  records,  26  October-10  November  197S. 
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Figure  G-4.  Spectra  of  the  Town  Creek  and  Jones  Creek  gage  water 
surface  records,  2-17  January  1975. 


APPENDIX  H 


DESCRIPTIVE  STATISTICS  OF  LITTORAL  PROCESS  PARAMETERS 

Measurements  were  made  by  observers  on  the  beach  (Tables  H^l,  H-2,  and 
H-3).  The  techniques  used  were  described  in  some  detail  by  Finley  (1976, 
p.  37).  The  parameters  are  defined  as  follows: 

VEL:  Longshore  current  velocity  (cm/s) 

WIND:  Windspeed,  2 meters  above  ground 

DIR:  Wind  direction  (azimuth) 

WINDL:  Longshore  component  of  windspeed  (in  tenths  of 
miles  per  hour) . Negative  sign  indicates  that  the  wind  was 
blowing  from  right  to  left. 

HGT:  Primary  breaker  height  (cm) 

PER:  Breaker  period  (s) 

ANGL:  Wave  orthogonal  angle  relative  to  the  beach.  Zero 
degrees  is  the  strike  of  the  beach  toward  the  left. 
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Table  H-l.  Summary  table  of  littoral  process  parameters  at 
station  DBI-25. 
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Table  H-l.  Summary  table  of  littoral  process  parameters  at 
station  DBI-2S.—  Continued 
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Table  H-l.  Summary  table  of  littoral  process  paraaeters  at 
station  DBI-2S. --Continued 
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Table  H-l.  Sunnary  table  of  littoral  process  parameters  at 
station  DBI-25.— Continued 
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0 
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Table  H-2.  Frequency  distribution  of  longshore  current  velocities 
at  station  DBI-25.  Negative  sign  indicates  current  to 
the  left  (north),  positive  currents  flow  to  the  right 
(south) . 
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0.403 

61.290 

5*  2 

to 

0.-06 

4,032 

20 

1 

153 

0,403 

61.694 

53  1 

11 
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0.906 

4.452 

25 

4 

162 

1.013 

65.323 

*6  2 

14 
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Table  H-3.  Frequency  distribution  of  wave  heights  at  station  DBI-2S. 


HOT 

rOEOUCNCY 

Com  raco 

OC4CCNT 

cum  Pf«cenr 

30 

9 

9 

3.629 

3.629 

33 

t 

10 

0.403 

4.037 

3* 

4 

14 

1.613 

5.645 

35 

15 

29 

6.049 

11.694 

3% 

7 

31 

0.006 

17.50" 

37 

1 

32 

0.403 

12.901 

31 

2 

34 

0.906 

13.710 

39 

1 

35 

0.403 

14.111 

*0 

20 

55 

9.065 

77.177 

*1 

2 

57 

0.406 

27.994 

5 

62 

2.016 

25.00" 

43 

3 

65 

1.210 

26.71" 

11 

1 

66 

0.403 

26.613 

*9 

6 

72 

2.419 

79.017 

*6 

4 

76 

1.613 

30.645 

*7 

2 

79 

0.006 

31.457 

*9 

5 

43 

2.016 

31.466 

49 

1 

44 

0.403 

31.67] 

SO 

20 

104 

9.065 

41.939 

SI 

1 

105 

0.403 

42.339 

52 

T 

112 

2.023 

45.161 

S3 

2 

114 

0.906 

45.969 

54 

5 

119 

2.016 

47.964 

55 

3 

122 

1.210 

49.194 

56 

6 

126 

2.419 

51.613 

57 

2 

110 

0.806 

52.419 

59 

3 

111 

1.210 

53.629 

60 

17 

150 

6.955 

60.494 

61 

4 

154 

1.613 

67.097 

6? 

2 

156 

0.606 

62.903 

63 

4 

160 

1.613 

64.516 

64 

4 

164 

1.613 

66.129 

65 

6 

170 

2.419 

46.546 

66 

1 

171 

0.403 

66.95? 

67 

1 

172 

0.403 

69.355 

66 

5 

177 

2.016 

71.371 

69 

1 

179 

0.403 

71.774 

70 

13 

191 

5.242 

77.016 

72 

5 

196 

2.016 

79.032 

71 

2 

194 

0.806 

79.639 

74 

2 

200 

0.606 

60.649 

75 

7 

207 

2.823 

93.466 

76 

1 

200 
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APPENDIX  I 

BEACH  PROFILES  AT  DUBIDUE  AND  NORTH  ISLANDS 

Beach  profiles  (Figs.  1-1  to  1-10)  were  measured  from  the  backshore 
to  approximately  2 feet  below  the  waterline.  The  backshore  reference 
height  remained  invariant  over  time.  A map  of  profile  locations  i§ 
given  in  Figure  28. 
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Figure  1-1.  Beach  profiles  at  station  DBI-6,  13  June  1975  to  18  May  1976. 
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Figure  1-9.  Beach  profiles  at  station  NI-14,  9 June  1975  to  18  May  1976 
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